Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,552,961 T→G 12.0% intergenic (+576/+85) PP_4882 → / ← PP_5717 iron ABC transporter permease/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,552,9610TG12.0% 276.5 / 30.9 110intergenic (+576/+85)PP_4882/PP_5717iron ABC transporter permease/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/46);  new base G (6/7);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-02

Only 100 of 118 total aligned reads displayed.

TTCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAG  >  NC_002947/5552812‑5553110
                                                                                                                                                        |                                                                                                                                                     
ttCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCt                                                                                                                                                       >  2:1656560/1‑150 (MQ=255)
   ggTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCg                                                                                                                                                    <  1:2282950/150‑1 (MQ=255)
   ggTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCg                                                                                                                                                    >  1:726882/1‑150 (MQ=255)
       cTCATTTTCGCGCTCTCATAAGTCACTGAGCATCAGTTGGTACGTGTGA‑‑GTTACATTTGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCaaa                                                                                                                                                <  2:63411/150‑1 (MQ=18)
        tCATTTTCGCGCTCTCATAAGTCACTGGGTTGCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                               >  1:63411/1‑150 (MQ=255)
        tCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAg                                                                                                                                               >  2:1383307/1‑150 (MQ=255)
          aTTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCt                                                                                                                                                       <  1:1262451/140‑1 (MQ=255)
           ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGccc                                                                                                                                            <  2:1220774/150‑1 (MQ=255)
           ttttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                            <  1:1597279/150‑1 (MQ=255)
            tttCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                           <  1:1405719/150‑1 (MQ=255)
              tCGCGCTCTCATAAGTCACTGCGCATCAGTCGAAACGTCTGGACCGTGCGCCGTGGACCAAGAAGCAGGCTTCTACTAAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCGGa                                                                                                                                         <  1:2341432/150‑1 (MQ=12)
                    tctcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAACAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACggg                                                                                                                                   >  1:319680/1‑150 (MQ=255)
                     ctcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                  <  1:1493357/150‑1 (MQ=255)
                     ctcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                  <  1:903756/150‑1 (MQ=255)
                      tcatcaGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGccc                                                                                                                                            <  2:1675178/139‑1 (MQ=18)
                      tcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                 <  1:2399106/150‑1 (MQ=255)
                      tcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgggg                                                                                                                                 >  2:1541956/1‑150 (MQ=255)
                      tcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                  >  2:592377/1‑149 (MQ=255)
                      tcATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTggg                                                                                                                                  <  1:592377/149‑1 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                     >  2:856588/1‑128 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTc                                                                                                                                                     <  1:856588/128‑1 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGtt                                                                                                                               >  2:1913484/1‑150 (MQ=255)
                        aTAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGtt                                                                                                                               >  1:2007824/1‑150 (MQ=255)
                              cACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTc                                                                                                                             >  1:1726936/1‑146 (MQ=255)
                              cACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTc                                                                                                                             <  2:1726936/146‑1 (MQ=255)
                                cTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtgg                                                                                                                       <  1:483871/150‑1 (MQ=255)
                                cTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCtgg                                                                                                                       <  2:1540662/150‑1 (MQ=255)
                                cTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGCCCg                                                                                                                                           >  1:992448/1‑130 (MQ=255)
                                 tGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                           <  2:2347259/129‑1 (MQ=255)
                                 tGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCg                                                                                                                                           >  1:2347259/1‑129 (MQ=255)
                                    gcATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTg                                                                                                                                    >  2:1339329/1‑133 (MQ=255)
                                    gcATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                   >  1:392252/1‑150 (MQ=255)
                                    gcATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                   >  1:391782/1‑150 (MQ=255)
                                    gcATCAGTTGGTACGTCTGG‑‑GGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGgc                                                                                                                   >  2:643894/1‑150 (MQ=255)
                                               tACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAgtgt                                                                                                        <  1:187935/150‑1 (MQ=255)
                                                aCGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGtt                                                                                                       <  2:391782/150‑1 (MQ=255)
                                                 cGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGttt                                                                                                      <  1:316537/150‑1 (MQ=255)
                                                  gTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                     >  2:1290922/1‑150 (MQ=255)
                                                      gg‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGGATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGaaa                                                                                                 >  1:1377880/1‑150 (MQ=255)
                                                              cAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGa                                                                                                                                         <  1:2407479/104‑1 (MQ=255)
                                                               aGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                          <  1:1383307/150‑1 (MQ=255)
                                                                  cGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                       <  1:1656560/150‑1 (MQ=255)
                                                                          cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACa                                                                                              >  2:1184911/1‑135 (MQ=255)
                                                                          cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACa                                                                                              <  1:1184911/135‑1 (MQ=255)
                                                                          cAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTg                                                                               >  2:25483/1‑150 (MQ=255)
                                                                          cAAGAAGCAGGCTTACGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                          >  2:629516/1‑139 (MQ=255)
                                                                          cAAGAAGCAGGCTTACGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATgg                                                                                          <  1:629516/139‑1 (MQ=255)
                                                                           aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgg                                                                                                                                   <  2:2388972/97‑1 (MQ=255)
                                                                           aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTgg                                                                                                                                   >  1:2388972/1‑97 (MQ=255)
                                                                           aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                       <  2:1776618/141‑1 (MQ=255)
                                                                           aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGtt                                                                                       >  1:1776618/1‑141 (MQ=255)
                                                                           aagaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCGAAGCAAAGcc                                                                                                                                             >  1:210128/1‑87 (MQ=17)
                                                                            agaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                             <  2:2007824/150‑1 (MQ=255)
                                                                             gaagCGGGCTTGTGCCGAGCAAGGTCTTGAGGAATTCGGACACGGAGGGCTAGACCGGCAGGTGCGGGGTTTGGTTCGCATAGTCCGGAAACTGGGGTTTGGGCTGGTGGCTGCTTGAGTGTTTAGAACACAAGGGGTTGCATTCTGttt                                                                            <  2:319680/150‑1 (MQ=255)
                                                                             gaagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                            <  1:25483/150‑1 (MQ=255)
                                                                              aagCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTg                                                                           <  1:1913484/150‑1 (MQ=255)
                                                                               agCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGa                                                                          <  1:97326/150‑1 (MQ=255)
                                                                                  aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                              <  2:1373724/111‑1 (MQ=255)
                                                                                  aGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCtt                                                                                                              >  1:1373724/1‑111 (MQ=255)
                                                                                    gCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGtt                                                                             >  1:1427448/1‑142 (MQ=255)
                                                                                      ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCGCGGAAACGGGGGTGTCAGCTGTTGGCTGCTTGAGTGTTTAGAAAACAATGGATTGCATTCTGTTTGACCAGGGc                                                                   >  1:1328723/1‑150 (MQ=17)
                                                                                      ttATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGATTGAGTGTTTAGAAAAAAATGGGTTGAATTCTGTTTGagatcgga                                                                   >  2:1260837/1‑143 (MQ=17)
                                                                                         tGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTa                                                                                                     <  2:1963950/113‑1 (MQ=255)
                                                                                             gAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACAGGGCGGGCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTTTGTTTGACCAGGGCTTACCTc                                                            >  1:151270/1‑150 (MQ=17)
                                                                                                cAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                         >  1:1403289/1‑150 (MQ=255)
                                                                                                cAAGGTCTTGAGGAATTCCGAACCAGAAGCACCGACGCGCAGGCGGGTAGCCCTTTTCGCAAAGCCCGGAAGCTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCCGGGCTTACCTCAAg                                                         >  1:304507/1‑150 (MQ=18)
                                                                                                  aGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTa                                                       >  2:220965/1‑150 (MQ=255)
                                                                                                    gTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                     <  2:2116334/150‑1 (MQ=255)
                                                                                                       ttGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTa                                                  >  2:959044/1‑150 (MQ=255)
                                                                                                        tGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTAt                                                 <  1:643894/150‑1 (MQ=255)
                                                                                                         gAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAAt                                                                                            >  2:234858/1‑106 (MQ=255)
                                                                                                           ggAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCa                                              >  1:188265/1‑150 (MQ=255)
                                                                                                           ggAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATAGGTTGAATTCAAGTTGACCAGTGCTTACCTCAAGTACTCTATCCa                                              >  1:1964593/1‑150 (MQ=18)
                                                                                                             aaTTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCaaa                                            <  1:1541956/150‑1 (MQ=255)
                                                                                                               ttCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAGCTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGa                                          <  2:304507/150‑1 (MQ=255)
                                                                                                                  cGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaaa                                       <  1:1290922/150‑1 (MQ=255)
                                                                                                                     aCCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCaaaaaaaaaaaa                                    <  2:1403289/150‑1 (MQ=255)
                                                                                                                       cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaa                                        <  2:763649/144‑1 (MQ=255)
                                                                                                                       cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGaaa                                        >  1:763649/1‑144 (MQ=255)
                                                                                                                       cAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAg                                  >  1:202684/1‑150 (MQ=255)
                                                                                                                         gAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGaa                                <  1:2171158/150‑1 (MQ=255)
                                                                                                                               aaGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTATTTGACTAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCa                          >  2:443926/1‑150 (MQ=255)
                                                                                                                                  aCCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                         <  2:850924/116‑1 (MQ=255)
                                                                                                                                  aCCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAg                                                         >  1:850924/1‑116 (MQ=255)
                                                                                                                                      gCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAAt                   >  1:811371/1‑150 (MQ=255)
                                                                                                                                      gCAGGTCCGGGCTTTGCATCGAAAAAACAGGAAAAAGGGGTTTCAGCTAGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCAAAAAAAAAAAAAAAAAGAACAATCAGTGCAAt                   <  2:1475981/150‑1 (MQ=21)
                                                                                                                                           tCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGt                                                                              >  2:952126/1‑86 (MQ=255)
                                                                                                                                             cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                            <  2:2285265/86‑1 (MQ=18)
                                                                                                                                             cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                            >  1:1179296/1‑86 (MQ=18)
                                                                                                                                             cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                            <  2:1179296/86‑1 (MQ=18)
                                                                                                                                             cGGGCTTTGCGAAGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGttt                                                                            >  1:2285265/1‑86 (MQ=18)
                                                                                                                                             cGGGCTTTGCGAAGCAAAGCCCAGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTa                                                                >  1:1732836/1‑98 (MQ=21)
                                                                                                                                                   tttGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTg       >  1:1048080/1‑150 (MQ=255)
                                                                                                                                                       cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGAc                                                                         >  2:304512/1‑80 (MQ=255)
                                                                                                                                                       cTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACGCGATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACa   >  1:2331034/1‑150 (MQ=255)
                                                                                                                                                        ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                     <  2:439111/99‑1 (MQ=255)
                                                                                                                                                        ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                     >  1:439111/1‑99 (MQ=255)
                                                                                                                                                        ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                     >  2:2282674/1‑99 (MQ=255)
                                                                                                                                                        ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTAct                                                     <  1:2282674/99‑1 (MQ=255)
                                                                                                                                                        ttCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAg  >  1:440111/1‑150 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
TTCGGTACTCATTTTCGCGCTCTCATAAGTCACTGCGCATCAGTTGGTACGTCTGG‑‑TGTACAGCCGGTTCACCAAGAAGCAGGCTTATGCTGAACAAGGTCTTGAGGAATTCCGAACCAGAAGCCAAGACCGGCAGGTCCGGGCTTTGCTTCGCAAAGCCCGGAAACTGGGGTTTCAGCTGGTGGCTGCTTGAGTGTTTAGAAAACAATGGGTTGCATTCTGTTTGACCAGGGCTTACCTCAAGTACTCTATCCAAAGAAAAAAAAGAACAATCAGTGCAATCAGCAAAATTTGTACAG  >  NC_002947/5552812‑5553110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: