Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,442,017 G→C 12.4% G479A (GGG→GCG)  thiC → thiamin (pyrimidine moiety) biosynthesis protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,442,0170G→C12.4% 87.4 ‑3.0 ‑0.1241G479A (GGG→GCG) thiCthiamin (pyrimidine moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  new base (20/10):  ref base (97/114):  total (117/124)
Fisher's exact test strand distribution p-value = 4.96e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.59e-01

Only 100 of 256 total aligned reads displayed.

AATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAG  >  NC_007779/3441917‑3442117
                                                                                                    |                                                                                                     
   atGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAt                                                                                                    <  1:10931457/101‑1 (MQ=255)
   atGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAt                                                                                                    >  1:6243427/1‑101 (MQ=255)
    tGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGAtt                                                                                                   <  1:520305/101‑1 (MQ=255)
     gACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTg                                                                                                  <  1:6429762/101‑1 (MQ=255)
      aCCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTgg                                                                                                 <  1:4503677/101‑1 (MQ=255)
       ccGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGt                                                                                                >  1:2376257/1‑101 (MQ=255)
        cGAGGAGTTAGAGCACTGCCACGAAGCGCCTTTTTACACTCTGGGGCCGCTAACTGCCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                               <  1:13462738/101‑1 (MQ=255)
        cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                               >  1:1669206/1‑101 (MQ=255)
        cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                               <  1:5101800/101‑1 (MQ=255)
        cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                               >  1:14095601/1‑101 (MQ=255)
        cGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTg                                                                                               >  1:10130077/1‑101 (MQ=255)
              gTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCga                                                                                         >  1:467872/1‑101 (MQ=255)
                tAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatg                                                                                       <  1:13805272/101‑1 (MQ=255)
                 agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                                                                      >  1:11993120/1‑101 (MQ=255)
                 agagCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCgatga                                                                                      >  1:13253561/1‑101 (MQ=255)
                     tggtgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgc                                                                                  >  1:8341683/10‑100 (MQ=255)
                     cACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTgg                                                                                 >  1:10226117/1‑101 (MQ=255)
                     cACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTgg                                                                                 >  1:1309990/1‑101 (MQ=255)
                        tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                                               <  1:5492815/95‑5 (MQ=255)
                        tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                                               <  1:11374226/95‑5 (MQ=255)
                        tgtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgca                                                                               >  1:11756615/7‑97 (MQ=255)
                         gtgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcag                                                                              >  1:573771/6‑96 (MQ=255)
                          tgaggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcagg                                                                             >  1:11026797/5‑95 (MQ=255)
                          gcACGAAGAGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACAAATTCACGTCGGGGATTGGTGCGGCGGTGTATGGATGGt                                                                            >  1:4012016/2‑101 (MQ=255)
                          ccACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGt                                                                            >  1:1517994/1‑101 (MQ=255)
                             ggAAGCTCCGTTCTACACCCTCGGCCCGCTGGCAACCGACATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctg                                                                          >  1:7270025/2‑92 (MQ=255)
                             cGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTg                                                                         >  1:8909127/1‑101 (MQ=255)
                              gAAGCTCCGTTCTACGCCCTCGGCCCGCTGGCAGCCGTCATCGCGCCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTgcgcaggctgg                                                                         <  1:8505637/101‑11 (MQ=255)
                              gAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTgg                                                                        <  1:7011556/101‑1 (MQ=255)
                                aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                                                      <  1:4235961/101‑1 (MQ=255)
                                aGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCt                                                                      >  1:2984189/1‑101 (MQ=255)
                                 gcgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGTGATGATTGGCTGGTTTGGCTg                                                                     <  1:7649784/101‑1 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                    >  1:6749588/1‑101 (MQ=255)
                                  cgcCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgc                                                                    <  1:7117127/101‑1 (MQ=255)
                                    ccGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTgcgc                                                                  <  1:2604202/101‑1 (MQ=255)
                                     cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                 >  1:8278608/1‑101 (MQ=255)
                                     cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                 >  1:9025261/1‑101 (MQ=255)
                                     cGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGcgcg                                                                 <  1:10712738/101‑1 (MQ=255)
                                      gTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGa                                                                >  1:505797/1‑101 (MQ=255)
                                       ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                               <  1:9631677/101‑1 (MQ=255)
                                       ttttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGAt                                                               <  1:5173937/101‑1 (MQ=255)
                                        tttACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATg                                                              <  1:4583599/101‑1 (MQ=255)
                                          tACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGct                                                            <  1:79210/101‑1 (MQ=255)
                                           acacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGctc                                                           <  1:10032789/101‑1 (MQ=255)
                                            cacTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGGTGCGCGATGctct                                                          >  1:3744810/1‑101 (MQ=255)
                                             acTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTg                                                         <  1:10364610/101‑1 (MQ=255)
                                              ctctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGt                                                        <  1:5581123/101‑1 (MQ=255)
                                              atctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGt                                                        <  1:3686834/100‑1 (MQ=255)
                                               tctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGtt                                                       <  1:3270501/101‑1 (MQ=255)
                                                ctGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTa                                                      <  1:13399857/101‑1 (MQ=255)
                                                  ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                    <  1:7721419/101‑1 (MQ=255)
                                                  ggggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACg                                                    >  1:13308615/1‑101 (MQ=255)
                                                    ggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTa                                                  <  1:1239999/101‑1 (MQ=255)
                                                    ggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTa                                                  >  1:5694349/1‑101 (MQ=255)
                                                    ggCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTa                                                  <  1:10572756/101‑1 (MQ=255)
                                                      ccGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAAc                                                >  1:8986606/1‑101 (MQ=255)
                                                       cGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACg                                               <  1:9736532/101‑1 (MQ=255)
                                                         cTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGcc                                             <  1:760251/101‑1 (MQ=255)
                                                           aaCTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaa                                           <  1:3415353/101‑1 (MQ=255)
                                                             cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                         <  1:4804565/101‑1 (MQ=255)
                                                             cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCaaaa                                         <  1:5378619/101‑1 (MQ=255)
                                                             cTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGAGTGGCTGGTTTGGGTGCGCGATGCTCTGTTACGTAACGCCaaaa                                         <  1:13033198/101‑1 (MQ=255)
                                                              tACCGATATTGCGCCGGGCTATGTCCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                        <  1:8721583/101‑1 (MQ=255)
                                                              tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                        <  1:7199490/101‑1 (MQ=255)
                                                              tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                        <  1:7395214/101‑1 (MQ=255)
                                                              tACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAag                                        <  1:2806844/101‑1 (MQ=255)
                                                                ccGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAagag                                      >  1:5298895/1‑101 (MQ=255)
                                                                 cGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGc                                     >  1:9184795/1‑101 (MQ=255)
                                                                  gacatcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCa                                    >  1:12888273/7‑101 (MQ=255)
                                                                    tatTTCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTTCGGCGATGATTGGCTGGTTTTGCTGCGCAATGGCCCGCTACGCAATGCCAAAAGAGCaac                                  >  1:13649561/1‑99 (MQ=255)
                                                                    tatTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATc                                  >  1:8082767/1‑101 (MQ=255)
                                                                     agcgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCt                                 <  1:2866559/98‑1 (MQ=255)
                                                                       tGCGCCGGGCTATGACCACTTCACGTCGGCGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTgg                               >  1:8142239/1‑101 (MQ=255)
                                                                       cgcgcCGGGCTATGACCACATCACGTCGGCTATCGGCGCGGCAATCATTGCGCAGGCTGGTACCGCCATGCTCTGCTACGTCACGCCGAAGGAGCACCTgg                               >  1:5645003/2‑101 (MQ=255)
                                                                          gcCGGGCTATGACCACTTCACGTCGGGGATTGGAGCGGCGTTGTTTGGCTGGTTTGGCTGCGCGACGCTCTGTTTCGTAACGCCAAAAGAGCATCTgggcc                            >  1:1385426/1‑99 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGGTGCTCTGTTACGCAACGACAAAAGAGCATCTgggagtg                          >  1:9298916/1‑97 (MQ=255)
                                                                            cGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTg                          >  1:11791363/1‑101 (MQ=255)
                                                                             gggCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGc                         <  1:5817937/101‑1 (MQ=255)
                                                                               gCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGccc                       >  1:12103995/1‑101 (MQ=255)
                                                                                cTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCa                      <  1:13576131/101‑1 (MQ=255)
                                                                                 tATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCaa                     >  1:12167644/1‑101 (MQ=255)
                                                                                   tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATa                   <  1:1049386/101‑1 (MQ=255)
                                                                                   tGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCCCGATGCCCTGTTACGTAACGCCAAAAAAGAAACGGGGTCTGCCCAATa                   >  1:6347068/1‑101 (MQ=255)
                                                                                    gACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATaa                  >  1:2204692/1‑101 (MQ=255)
                                                                                      ccACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaag                <  1:13273840/101‑1 (MQ=255)
                                                                                       cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga               >  1:7730998/1‑101 (MQ=255)
                                                                                       cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga               >  1:11290704/1‑101 (MQ=255)
                                                                                       cACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAaaga               >  1:8627172/1‑101 (MQ=255)
                                                                                          ttCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGTCCAATAAagaaga            <  1:1947835/101‑1 (MQ=255)
                                                                                            cACGTCGGGGATTGGTGCGGCGATGATTGGCTGGGTTGGCTGCGGGAGGCTCTGTTACGGAAAGCCAAAAGAGCATCCGGGTCTGCCCAATGAAGAAGATg          >  1:319961/1‑101 (MQ=255)
                                                                                             aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt         <  1:4174154/101‑1 (MQ=255)
                                                                                             aCGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGt         <  1:13998729/101‑1 (MQ=255)
                                                                                               gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa       <  1:6704658/101‑1 (MQ=255)
                                                                                               gTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTa       <  1:352040/101‑1 (MQ=255)
                                                                                                tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa      >  1:12202770/1‑101 (MQ=255)
                                                                                                tCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTaa      <  1:4894091/101‑1 (MQ=255)
                                                                                                 cGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGTCCAATAAAGAAGATGTTAAg     <  1:1785350/101‑1 (MQ=255)
                                                                                                  ggggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGc    <  1:9592434/101‑1 (MQ=255)
                                                                                                   gggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCa   >  1:12786883/1‑101 (MQ=255)
                                                                                                    ggATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAg  <  1:8534782/101‑1 (MQ=255)
                                                                                                    |                                                                                                     
AATATGACCGAGGAGTTAGAGCACTGCCACGAAGCGCCGTTTTACACTCTGGGGCCGCTAACTACCGATATTGCGCCGGGCTATGACCACTTCACGTCGGGGATTGGTGCGGCGATGATTGGCTGGTTTGGCTGCGCGATGCTCTGTTACGTAACGCCAAAAGAGCATCTGGGTCTGCCCAATAAAGAAGATGTTAAGCAG  >  NC_007779/3441917‑3442117

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.