Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,329,444 T→A 6.5% L10Q (CTG→CAG)  PP_4690 → Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,329,4440TA 6.5% 391.4 / 12.4 123L10Q (CTG→CAG) PP_4690Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base T (64/51);  new base A (3/5);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 123 total aligned reads displayed.

TTGTGGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGAT  >  NC_002947/5329298‑5329586
                                                                                                                                                  |                                                                                                                                              
ttGTGGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGc                                                                                                                                             <  1:242874/150‑1 (MQ=255)
    ggCGAGCGGCGGAACTGCGCACGAGAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGGCAAGCCATCACGTGTTTCCCGGGTCATCCTTGCCATGCACTTCCTTTGTACATCTATCGGCCTGGCCGAg                                                                                                                                         <  2:739523/150‑1 (MQ=255)
     gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                        <  1:449824/150‑1 (MQ=255)
     gCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                        <  2:1801257/150‑1 (MQ=255)
      ggAGCTGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATCCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAggg                                                                                                                                       <  1:2042560/149‑1 (MQ=255)
        aGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                     <  2:2224027/150‑3 (MQ=255)
        aGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                     <  2:1195594/150‑3 (MQ=255)
                aaCTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCcgc                                                                                                                             <  2:499767/150‑1 (MQ=255)
                 aCTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                     <  1:755358/141‑3 (MQ=255)
                   tGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGcc                                                                                                                          >  2:1285425/1‑150 (MQ=255)
                   tGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATATGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGcc                                                                                                                          >  1:2299685/1‑150 (MQ=255)
                      gcACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTc                                                                                                                       >  2:1412657/1‑150 (MQ=255)
                            tAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTa                                                                                                                 <  2:1853326/150‑1 (MQ=255)
                            tAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGAAAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTa                                                                                                                 >  2:107044/1‑150 (MQ=255)
                                 aGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATGCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACgg                                                                                                            >  1:544015/1‑150 (MQ=255)
                                  gCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGc                                                                                                           >  1:1274930/1‑150 (MQ=255)
                                         gcAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAAc                                                                                                    >  2:2391355/1‑150 (MQ=255)
                                          cAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACt                                                                                                   >  2:1881218/1‑150 (MQ=255)
                                           agagAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTg                                                                                                  >  1:905876/1‑150 (MQ=255)
                                           agagAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTg                                                                                                  >  2:1044836/1‑150 (MQ=255)
                                             agaATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGtt                                                                                                >  2:240315/1‑150 (MQ=255)
                                               aaTGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCg                                                                                              >  1:537800/1‑150 (MQ=255)
                                               aaTGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCg                                                                                              >  1:2080292/1‑150 (MQ=255)
                                                   cGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgt                                                                                          <  1:1044836/150‑1 (MQ=255)
                                                     aTCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgc                                                                                        >  1:568222/1‑150 (MQ=255)
                                                     aTCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgc                                                                                        <  1:240315/150‑1 (MQ=255)
                                                      tCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgcg                                                                                       >  2:1161659/1‑150 (MQ=255)
                                                      tCAATTAGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgcg                                                                                       >  1:2020025/1‑150 (MQ=255)
                                                         aTTTGATAATTATTCGCATTGGATCGCCATGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCg                                                                                    >  1:932971/1‑150 (MQ=255)
                                                             gATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCg                                                                                                                   <  1:1575672/115‑1 (MQ=255)
                                                             gATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCg                                                                                                                   >  2:1575672/1‑115 (MQ=255)
                                                              aTAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAg                                                                               >  2:1433358/1‑150 (MQ=255)
                                                              aTAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTAACGGGTCCTACTTGCCATGAAAATCCTTTGTACATCCCTCGGGCTGGCCTAGTGCCAAATCCGGGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAg                                                                               >  2:1478475/1‑150 (MQ=255)
                                                                aaTTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGTACTGTTCGGCGTGCGTCGCCAAGGc                                                                             <  2:1954764/150‑1 (MQ=255)
                                                                     ttCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCa                                                                                                                      >  2:1192420/1‑104 (MQ=255)
                                                                     ttCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGgt                                                                        <  1:484784/150‑1 (MQ=255)
                                                                     ttCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGgt                                                                        >  2:2386321/1‑150 (MQ=255)
                                                                        gCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgta                                                                     >  1:235456/1‑150 (MQ=255)
                                                                         cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                    <  1:1285425/150‑1 (MQ=255)
                                                                         cATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                    <  1:1412657/150‑1 (MQ=255)
                                                                            tGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctg                                                                 <  1:1431795/150‑1 (MQ=255)
                                                                              gAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCCCGCTGGTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgta                                                               >  1:1933772/1‑150 (MQ=255)
                                                                               aGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgtac                                                              <  1:226417/150‑1 (MQ=255)
                                                                                   tCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCa                                                          >  1:650660/1‑150 (MQ=255)
                                                                                        ccATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGc                                                     >  1:1345489/1‑150 (MQ=255)
                                                                                         cATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGGGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgta                                                               >  2:583717/1‑139 (MQ=255)
                                                                                          aTCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGcc                                                                                                          >  1:906925/1‑95 (MQ=255)
                                                                                          aTCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGcc                                                                                                          <  2:906925/95‑1 (MQ=255)
                                                                                           tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAgg                                                                                                                                        <  2:463761/64‑1 (MQ=255)
                                                                                           tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGc                                                  <  1:1826675/150‑1 (MQ=255)
                                                                                           tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGc                                                  <  2:2043410/150‑1 (MQ=255)
                                                                                           tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                     <  1:798495/67‑3 (MQ=255)
                                                                                           tCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCAGGCCGAgggat                                                                                                                                     >  2:798495/1‑65 (MQ=255)
                                                                                            cACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGcc                                                 <  1:1433358/150‑1 (MQ=255)
                                                                                               gtgtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCaaa                                                                                                                                  >  1:125126/1‑66 (MQ=255)
                                                                                               gtgtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCaaa                                                                                                                                  <  2:125126/66‑1 (MQ=255)
                                                                                               gtgtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACGGCTGc                                                  >  2:177594/1‑146 (MQ=255)
                                                                                               gtgtTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGc                                                  <  1:177594/146‑1 (MQ=255)
                                                                                                   ttCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCAGCGCCTTCAGCGTAGACGGCCTCGAACTTTTCGGCGTGCGTCTCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCg                                          <  2:2080292/150‑1 (MQ=255)
                                                                                                    tCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCgg                                         <  1:107044/150‑1 (MQ=255)
                                                                                                    tCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCgg                                         <  1:983329/150‑1 (MQ=255)
                                                                                                      aCGGGTCATCCTTGCCATGCACTTCCTTTGTCCCTACCTCGGCCTTGCGCATGTCCAAAGCCGCGCCTTCAAAGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGa                                       <  1:1142503/150‑1 (MQ=255)
                                                                                                      aCGGGTCATCCTTGCAATGCACTTCCTTTGTCACTCCCTCGGCCTGGCCGAGGGCCACAGCCGCGCCTTCAGCGTAGACGGCTTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGa                                       <  2:1933772/150‑1 (MQ=255)
                                                                                                        gggTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacc                                                                   <  1:1261954/120‑1 (MQ=255)
                                                                                                              tCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGc                               >  2:556619/1‑150 (MQ=255)
                                                                                                              tCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGc                               >  1:322012/1‑150 (MQ=255)
                                                                                                                cTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTg                             >  1:221700/1‑150 (MQ=255)
                                                                                                                 ttGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGa                            >  1:1846690/1‑150 (MQ=255)
                                                                                                                    ccATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACt                         >  2:957053/1‑150 (MQ=255)
                                                                                                                    ccATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACt                         >  2:353748/1‑150 (MQ=255)
                                                                                                                      aTGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTgg                       >  1:190419/1‑150 (MQ=255)
                                                                                                                       tGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTggg                      <  2:235456/150‑1 (MQ=255)
                                                                                                                        gCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGc                     >  2:1109765/1‑150 (MQ=255)
                                                                                                                         cACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCa                    <  1:1903810/150‑1 (MQ=255)
                                                                                                                          aCTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtac                                                                    <  1:2044830/101‑1 (MQ=255)
                                                                                                                           cTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCaga                  <  1:1881218/150‑1 (MQ=255)
                                                                                                                              ccTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGc                                                                                                                    >  1:907206/1‑49 (MQ=255)
                                                                                                                              ccTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGc                                                                                                                    <  2:907206/49‑1 (MQ=255)
                                                                                                                               cTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAc              >  2:787669/1‑150 (MQ=255)
                                                                                                                                tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAAc                                                                                                    <  1:1377976/63‑1 (MQ=255)
                                                                                                                                tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAAc                                                                                                    >  2:1377976/1‑63 (MQ=255)
                                                                                                                                tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgtacc                                                             >  2:764298/1‑102 (MQ=255)
                                                                                                                                tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGtgtacctgtacc                                                             <  1:764298/102‑1 (MQ=255)
                                                                                                                                tttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCAGCGCCTTAAGAGTAAACAGCCTCGAACTGTTCGGCTGGCGTCGCCCAAGCCAGGTGTACCTGTACCGCAACCGCTCCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGCcg             >  1:46822/1‑150 (MQ=255)
                                                                                                                                 ttGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAcgc            <  2:2226669/150‑1 (MQ=255)
                                                                                                                                  tGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGAcgcg           >  2:1009263/1‑150 (MQ=255)
                                                                                                                                   gTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGt          <  2:1345489/150‑1 (MQ=255)
                                                                                                                                    tACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCaa                                                         >  2:1862007/1‑102 (MQ=255)
                                                                                                                                    tACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCgg                                         >  2:2257774/1‑118 (MQ=255)
                                                                                                                                    tACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCgg                                         <  1:2257774/118‑1 (MQ=255)
                                                                                                                                     aCATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAg                                                                          >  1:1157822/1‑84 (MQ=255)
                                                                                                                                      cATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTc       <  2:1846690/150‑1 (MQ=255)
                                                                                                                                        tCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCa                                                                           >  2:853085/1‑80 (MQ=255)
                                                                                                                                         cccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTgg    >  2:483994/1‑150 (MQ=255)
                                                                                                                                          ccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTCGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTCCCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGa   <  2:650660/150‑1 (MQ=255)
                                                                                                                                          ccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGa   >  2:358280/1‑150 (MQ=255)
                                                                                                                                          ccTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGa   <  2:221700/150‑1 (MQ=255)
                                                                                                                                           cTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATc                                   <  2:936186/117‑1 (MQ=255)
                                                                                                                                           cTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGAt  >  1:1087920/1‑150 (MQ=255)
                                                                                                                                           cTCGGCCAGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGgcgtgc                                                                                        >  2:2001943/1‑64 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TTGTGGCGAGCGGCGGAACTGCGCACCATAGGCAGCCGGGCGCAGAGAATGCGATCAATTTGATAATTATTCGCATTGGAGCGTCAAGCCATCACGTGTTTCACGGGTCATCCTTGCCATGCACTTCCTTTGTACATCCCTCGGCCTGGCCGAGGGCCAAAGCCGCGCCTTCAGCGTAGACGGCCTCGAACTGTTCGGCGTGCGTCGCCAAGGCCAGGTGTACCTGTACCGCAACCGCTGCCCGCACCGGGACATCCCGCTGAACTGGGCAGAGGACGCGTTTCTGGAT  >  NC_002947/5329298‑5329586

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