Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,645,945 A→T 5.6% L32Q (CTG→CAG)  PP_3212 ← Rieske 2Fe‑2S family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,645,9450AT 5.6% 348.8 / 8.8 108L32Q (CTG→CAG) PP_3212Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (50/52);  new base T (3/3);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 113 total aligned reads displayed.

AGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACC  >  NC_002947/3645796‑3646093
                                                                                                                                                     |                                                                                                                                                    
aGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCca                                                                                                                                                      <  1:1485435/150‑1 (MQ=255)
aGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATGGCCCTCAGCACTGACGAAGCCTTTGGACAGGGTGCCGGTGCGGTGCCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAAAACGAGCGCTTCACCca                                                                                                                                                      >  1:1977036/1‑150 (MQ=255)
  tCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCcagc                                                                                                                                                    >  1:1830692/1‑150 (MQ=255)
   cGTATGCCCAGCCGTGATCCCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGACTTCCCAACGGTGTACCATCGCTACGCTGCCACAACACCAGCGCTTCACCcagca                                                                                                                                                   <  2:26784/150‑1 (MQ=255)
     tATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTgg                                                                                                                                                 <  2:1088650/150‑1 (MQ=255)
     tATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAgg                                                                                                                                                 <  1:95854/150‑1 (MQ=255)
     tATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAgg                                                                                                                                                 <  1:177579/150‑1 (MQ=255)
      aTTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGAACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAggg                                                                                                                                                >  2:2082533/1‑150 (MQ=255)
          ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGaa                                                                                                                                            >  1:403538/1‑150 (MQ=255)
          ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGaa                                                                                                                                            >  2:584612/1‑150 (MQ=255)
                 tGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTg                                                                                                                                     <  2:1850396/150‑1 (MQ=255)
                 tGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTg                                                                                                                                     <  2:837683/150‑1 (MQ=255)
                   aTACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                   <  2:182498/150‑1 (MQ=255)
                   aTACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                   >  1:434136/1‑150 (MQ=255)
                    tACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTgggg                                                                                                                                  <  2:403538/150‑1 (MQ=255)
                      cccGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGcc                                                                                                                                <  2:2149169/150‑1 (MQ=255)
                              cGATATTGCCCTCAGCACTGACGAAGCACCTGGACCATCTGGCTGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGa                                                                                                                        <  2:804241/150‑1 (MQ=255)
                                   ttGCCCTCAGCACTGACGAAGCCTTTGGACTGTTTGGCGGTGCGGCGACAGCAGCGGTCCCGCATTCCCACCCGGGTACCCTCGGGACCCTTCCACAACACCAGCGCTTGACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   <  2:131146/150‑1 (MQ=255)
                                   ttGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   >  2:2025647/1‑150 (MQ=255)
                                   ttGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   <  1:354880/150‑1 (MQ=255)
                                         tCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATa                                                                                                             >  1:849609/1‑150 (MQ=255)
                                           aGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAgg                                                                                                           <  2:857878/150‑1 (MQ=255)
                                           aGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAgg                                                                                                           <  1:1427977/150‑1 (MQ=255)
                                             cACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGt                                                                                                                                               <  1:359362/112‑1 (MQ=255)
                                             cACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGt                                                                                                                                               >  2:359362/1‑112 (MQ=255)
                                                tGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                   <  2:255570/121‑1 (MQ=255)
                                                tGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                   >  1:255570/1‑121 (MQ=255)
                                                tGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCa                                                                                                      >  1:1510768/1‑150 (MQ=255)
                                                  aCGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGATGAATGGCTTGGGGCCGTCTTCTAGGGCTGTCATAGGGAGCacg                                                                                                    >  2:768203/1‑148 (MQ=255)
                                                    gAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGc                                                                                                  >  1:1263068/1‑150 (MQ=255)
                                                    gAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGc                                                                                                  >  2:105388/1‑150 (MQ=255)
                                                         cTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAg                                                                                                                       <  1:2181315/124‑1 (MQ=255)
                                                         cTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAg                                                                                                                       >  2:2181315/1‑124 (MQ=255)
                                                          tttGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCa                                                                                            <  1:1760548/150‑1 (MQ=255)
                                                          tttGGACAGTTTGGCGCTGCGGCGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGCGACGCTTCCACAACACCAGCGTTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCa                                                                                            <  1:10780/150‑1 (MQ=255)
                                                              gACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGaaa                                                                                        >  1:1573767/1‑150 (MQ=255)
                                                                   tttGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgcg                                                                                   >  1:2153182/1‑150 (MQ=255)
                                                                      ggCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCaa                                                                                <  2:1689491/150‑1 (MQ=255)
                                                                           tgcggtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAg                                                                           <  1:2025647/150‑1 (MQ=255)
                                                                           tgcggtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAg                                                                           <  2:499594/150‑1 (MQ=255)
                                                                           tgcggtgGCAGCAGCGGTCCCGCATTGCCACCGCGGCCGCGGCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAg                                                                           <  1:126592/150‑1 (MQ=255)
                                                                             cggtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgt                                                                         <  1:584612/150‑1 (MQ=255)
                                                                               gtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   <  2:1474757/106‑1 (MQ=255)
                                                                               gtgGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   >  1:1474757/1‑106 (MQ=255)
                                                                                  gcagcaGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                                    <  2:1573767/150‑1 (MQ=255)
                                                                                   cagcagCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg                                                                   >  1:1800419/1‑150 (MQ=255)
                                                                                    agcagcGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTgg                                                                  >  2:1661830/1‑150 (MQ=255)
                                                                                     gcagcGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGt                                                                 >  2:732548/1‑150 (MQ=255)
                                                                                        gcGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAAc                                                              <  2:1830692/150‑1 (MQ=255)
                                                                                            tCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAAc                                                          >  2:811811/1‑150 (MQ=255)
                                                                                             cccGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAg                                                                                                     <  2:924594/106‑1 (MQ=255)
                                                                                              ccGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAAAAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACAt                                                        >  1:263726/1‑150 (MQ=255)
                                                                                               cGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACg                                                                                      >  2:2399199/1‑119 (MQ=255)
                                                                                               cGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACg                                                                                      <  1:2399199/119‑1 (MQ=255)
                                                                                                gCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATaa                                                      >  2:116874/1‑150 (MQ=255)
                                                                                                 cATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGt                                                                 >  1:1357877/1‑138 (MQ=255)
                                                                                                 cATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGt                                                                 <  2:1357877/138‑1 (MQ=255)
                                                                                                      ccACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGa                                                <  1:435708/150‑1 (MQ=255)
                                                                                                      ccACCGGGGTACCATCGGGACGAAACCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGa                                                >  1:933449/1‑150 (MQ=255)
                                                                                                        aCCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCa                                                   <  1:30849/145‑1 (MQ=255)
                                                                                                        aCCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCa                                                   >  2:30849/1‑145 (MQ=255)
                                                                                                              gTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTcac                                        >  1:763278/1‑150 (MQ=255)
                                                                                                               tACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTcaca                                       <  1:811811/150‑1 (MQ=255)
                                                                                                                   aTCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGt                                   >  1:1743659/1‑150 (MQ=255)
                                                                                                                        gACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCa                                                                                            >  2:1405304/1‑88 (MQ=255)
                                                                                                                        gACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCa                                                                                            <  1:1405304/88‑1 (MQ=255)
                                                                                                                           gCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAAcc                                                             <  2:580520/116‑1 (MQ=255)
                                                                                                                           gCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAAcc                                                             >  1:580520/1‑116 (MQ=255)
                                                                                                                            cTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGt                          <  2:347669/150‑1 (MQ=255)
                                                                                                                             ttCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   >  1:60577/1‑60 (MQ=255)
                                                                                                                             ttCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                   <  2:60577/60‑1 (MQ=255)
                                                                                                                             ttCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGtt                         >  2:1747017/1‑150 (MQ=255)
                                                                                                                                cacaACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATg                      <  2:1977036/150‑1 (MQ=255)
                                                                                                                                cacaACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGGCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAAACAACATAATCAGGATTCCTCACAACGTTGAGCTTGTTAtt                      >  2:1084844/1‑149 (MQ=255)
                                                                                                                                 acaacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAgg                                                                                                    <  2:942288/71‑1 (MQ=255)
                                                                                                                                 acaacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAgg                                                                                                    >  1:942288/1‑71 (MQ=255)
                                                                                                                                 acaacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATgg                     >  1:1347709/1‑150 (MQ=255)
                                                                                                                                  caacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTcaca                                       <  2:102920/131‑1 (MQ=255)
                                                                                                                                  caacaCCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTcaca                                       >  1:102920/1‑131 (MQ=255)
                                                                                                                                    acacCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgtgt                  <  1:1159289/150‑1 (MQ=255)
                                                                                                                                    acacCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgtgt                  >  2:1690804/1‑150 (MQ=255)
                                                                                                                                     cacCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGTCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTTAGCGGATTATGGTGtt                 >  1:498271/1‑150 (MQ=255)
                                                                                                                                      acCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCa                                                                                                      <  2:423678/64‑1 (MQ=255)
                                                                                                                                      acCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCa                                                                                                      >  1:423678/1‑64 (MQ=255)
                                                                                                                                        cAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAg              <  2:1263068/150‑1 (MQ=255)
                                                                                                                                         accgcTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTTCTCCCTATCTCTCATTGGGTGTAGGTCATACCAAACATCACGCAATACATGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGc             <  1:542137/148‑1 (MQ=255)
                                                                                                                                         aGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgtgt                  >  2:223715/1‑145 (MQ=255)
                                                                                                                                         aGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgtgt                  <  1:223715/145‑1 (MQ=255)
                                                                                                                                            cattCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgt                    >  1:690709/3‑140 (MQ=255)
                                                                                                                                            cattCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgt                    <  2:690709/138‑1 (MQ=255)
                                                                                                                                             cTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGGATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCt         <  2:933449/150‑1 (MQ=255)
                                                                                                                                             cTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCt         <  2:1743659/150‑1 (MQ=255)
                                                                                                                                             attCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                                    >  1:1287823/2‑91 (MQ=255)
                                                                                                                                             attCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                                    <  2:1287823/90‑1 (MQ=255)
                                                                                                                                             attCACCCTGCTGGGAGAATGGCTTGGGGCCGTCTTCGCAAACACAGATTGTGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCGCACATCGTTTAGCTTGTTATGGGGTGCGGCGGCt         >  2:944343/2‑150 (MQ=255)
                                                                                                                                              ttCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgt                                                                         <  1:1596361/85‑1 (MQ=255)
                                                                                                                                              ttCACCCAGCAGGGTCAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGa                                <  1:1934960/126‑1 (MQ=255)
                                                                                                                                                 aCCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCg     >  1:455850/1‑150 (MQ=255)
                                                                                                                                                  cccAGCAGGGCGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGa    >  2:314640/1‑150 (MQ=255)
                                                                                                                                                    cagcagGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGAcc  >  1:1788202/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGTCGTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGCAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACC  >  NC_002947/3645796‑3646093

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