Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,070,116 A→T 7.8% intergenic (‑96/+25) plsC ← / ← PP_1845 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,070,1160AT 7.8% 528.9 / 23.9 166intergenic (‑96/+25)plsC/PP_18451‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase/enoyl‑CoA hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base A (76/76);  new base T (8/5);  total (85/81)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 169 total aligned reads displayed.

CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2069968‑2070264
                                                                                                                                                     |                                                                                                                                                    
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCaa                                                                                                                                                     <  2:2844697/150‑1 (MQ=255)
ccAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCaa                                                                                                                                                     <  1:1515514/150‑1 (MQ=255)
     ggCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCcg                                                                                                                                                <  1:2948727/150‑1 (MQ=255)
        aaGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGa                                                                                                                                             <  1:3671855/150‑1 (MQ=255)
        aaGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGa                                                                                                                                             >  1:418150/1‑150 (MQ=255)
          gAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaa                                                                                                                                           >  1:259382/1‑150 (MQ=255)
           aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  2:2509212/150‑1 (MQ=255)
            aaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                          <  1:2360447/150‑2 (MQ=255)
             aaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAgg                                                                                                                                        <  2:4897843/150‑1 (MQ=255)
                 cAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGccc                                                                                                                                    <  1:3034495/150‑1 (MQ=255)
                  aGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCg                                                                                                                                   <  2:4314713/150‑1 (MQ=255)
                            ggAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                         <  1:4396822/150‑1 (MQ=255)
                                 aTACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGACAAGGGCCCGTGCGGCTAGCCAata                                                                                                                    >  2:2107175/1‑150 (MQ=255)
                                 aTACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAata                                                                                                                    <  2:4981496/150‑1 (MQ=255)
                                 aTACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGGGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAata                                                                                                                    >  2:2107220/1‑150 (MQ=255)
                                  tACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGACAGGGGCCCGTGCGGCTAGCCAatag                                                                                                                   >  2:3224031/1‑149 (MQ=255)
                                   aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtg                                                                                                                  >  1:5435255/1‑150 (MQ=255)
                                   aCGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtg                                                                                                                  >  1:5435361/1‑150 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  2:3698183/150‑1 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  1:3126173/150‑1 (MQ=255)
                                      cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  1:3160505/150‑1 (MQ=255)
                                       aaTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctgc                                                                                                              >  1:885876/1‑150 (MQ=255)
                                        aTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                             >  2:644129/1‑150 (MQ=255)
                                           aaTAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGAt                                                                                                          >  1:3144051/1‑150 (MQ=255)
                                              aaaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGTACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGTGGGCAAGCCCGCGACAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCa                                                                                                       >  2:4858721/1‑150 (MQ=255)
                                              aaaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgc                                                                                                                 <  2:150558/140‑1 (MQ=255)
                                              aaaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgc                                                                                                                 >  1:150558/1‑140 (MQ=255)
                                                 gCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCt                                                                                                    >  1:4823593/1‑150 (MQ=255)
                                                     ggTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGAt                                                                                                <  1:2555075/150‑1 (MQ=255)
                                                        aCGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGc                                                                                             <  1:4875017/150‑1 (MQ=255)
                                                         cGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                            <  2:4417480/150‑1 (MQ=255)
                                                          gACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               <  1:2657063/146‑1 (MQ=255)
                                                          gACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               >  2:2657063/1‑146 (MQ=255)
                                                          gACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCt                                                                                           >  2:622028/1‑150 (MQ=255)
                                                           aCTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCtt                                                                                          <  1:3254353/150‑1 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                        >  2:1907650/1‑150 (MQ=255)
                                                             tCACACCAGACGGTCGCGGGGCAGTGCACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                        >  2:4166529/1‑150 (MQ=255)
                                                              cacaCCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTc                                                                                       <  1:622028/150‑1 (MQ=255)
                                                               acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTcc                                                                                      >  1:1785367/1‑150 (MQ=255)
                                                                 acCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAg                                                                                    <  2:5178954/150‑1 (MQ=255)
                                                                  ccAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  2:1439159/95‑1 (MQ=255)
                                                                  ccAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgc                                                                                   <  1:1085466/150‑1 (MQ=255)
                                                                    aGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgcgc                                                                                 <  2:1554705/150‑1 (MQ=255)
                                                                      aCGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCt                                                                               <  1:1654169/150‑1 (MQ=255)
                                                                      aCGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCt                                                                               <  2:1577328/150‑1 (MQ=255)
                                                                         gTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCca                                                                            <  2:242999/150‑1 (MQ=255)
                                                                         gTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCca                                                                            <  2:202720/150‑1 (MQ=255)
                                                                                ggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGccc                                                                     >  1:3498036/1‑150 (MQ=255)
                                                                                 gCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCt                                                                    >  2:4617738/1‑150 (MQ=255)
                                                                                  cAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCtt                                                                   <  1:4792871/150‑1 (MQ=255)
                                                                                     tgtACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGATAAGGGCCCGTTCGGCTAGCCAATTTTCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                <  2:4450218/150‑1 (MQ=255)
                                                                                     tgtACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                >  1:3651555/1‑150 (MQ=255)
                                                                                       tACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgcgc                                                              >  2:1308189/1‑150 (MQ=255)
                                                                                          gggTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                         <  2:1613045/88‑1 (MQ=255)
                                                                                            gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                            >  2:2866198/1‑115 (MQ=255)
                                                                                            gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                            <  1:2866198/115‑1 (MQ=255)
                                                                                            gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                            >  2:2866273/1‑115 (MQ=255)
                                                                                            gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCttt                                                                  >  2:5183056/1‑141 (MQ=255)
                                                                                            gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACttt                                                         >  1:2333633/1‑150 (MQ=255)
                                                                                             tgtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTc                                                        >  1:4285874/1‑150 (MQ=255)
                                                                                             tgtgCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTc                                                        >  1:4663524/1‑150 (MQ=255)
                                                                                                    tggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                                        <  1:4636406/63‑1 (MQ=25)
                                                                                                    tggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTAGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACATTAAGCTTgg                                                 >  1:1384157/1‑150 (MQ=255)
                                                                                                     ggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGccctga                                                <  2:5409480/150‑6 (MQ=255)
                                                                                                      gctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgtg                                               >  2:843072/1‑150 (MQ=255)
                                                                                                      gctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGAGCCTCCAGGAGCCCTTTGCGGACTTTCAGCTTGgtg                                               <  2:4524202/150‑1 (MQ=255)
                                                                                                          gggCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCt                                           >  2:5162265/1‑150 (MQ=255)
                                                                                                            gCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtg                                         >  2:3076485/1‑150 (MQ=255)
                                                                                                              aaaCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTtgtg                                       >  1:1829100/1‑150 (MQ=255)
                                                                                                               aaCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTgc                                      <  1:3759313/150‑1 (MQ=255)
                                                                                                                  cTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGcc                                   >  1:1632216/1‑150 (MQ=255)
                                                                                                                  cTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGcc                                   >  2:3409859/1‑150 (MQ=255)
                                                                                                                  cTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAACCCCGCGAAAAGGGCCCGTGCGGATAGCCAATAAGCTGCAGATCCACCTCGATCGCCATGTCCAGCGCCGCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTGGTGCGcc                                   >  2:4905804/1‑150 (MQ=255)
                                                                                                                   ttGGTGCTGGCCTTTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                                  <  1:1438072/150‑1 (MQ=255)
                                                                                                                   ttGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                                  <  2:225290/150‑1 (MQ=255)
                                                                                                                   ttGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCa                                  <  2:1914192/150‑1 (MQ=255)
                                                                                                                    tGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAg                                 >  1:1606440/1‑150 (MQ=255)
                                                                                                                           ggCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                          <  1:5102634/150‑1 (MQ=255)
                                                                                                                           ggCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTc                          <  1:4820180/150‑1 (MQ=255)
                                                                                                                             ccTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCACCAGCATGTTCAg                        <  2:2567040/150‑1 (MQ=255)
                                                                                                                              cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                       >  1:47202/1‑150 (MQ=255)
                                                                                                                              cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCGTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGc                       >  1:1139266/1‑150 (MQ=255)
                                                                                                                               tGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGct                      <  1:5364640/150‑1 (MQ=255)
                                                                                                                                   tcGGCCCTTTCGCGGGCAAGCCCGCGAAATGGGCCCGTGCTGCTAGCGAATATTCTGCTGATGCAGCTCGATCGCATTGTCGAGCGCCTACAGGAGGCCTTTTCGCACTTTCAGATTGGTGTTCTTGTGCGCCAGGATGTTCTGcttctt                  <  2:4562667/149‑1 (MQ=255)
                                                                                                                                   aCGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt                  >  1:5193595/1‑150 (MQ=255)
                                                                                                                                     ggCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                <  1:1616397/150‑1 (MQ=255)
                                                                                                                                      gcccttttCGCGGGCTTGCCCGCGAAAGGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                >  2:1483453/6‑150 (MQ=255)
                                                                                                                                      gCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAg               <  2:4663524/150‑1 (MQ=255)
                                                                                                                                       cccTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGc              >  2:5053472/1‑150 (MQ=255)
                                                                                                                                       cccTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCGTGTGCGCCAGCATGTTCAGCTTCTTCAGc              >  2:4401274/1‑150 (MQ=255)
                                                                                                                                        ccttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCa                <  1:4706243/144‑1 (MQ=255)
                                                                                                                                        ccttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCt             >  2:4854920/4‑150 (MQ=255)
                                                                                                                                        ccttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCt             >  1:2238579/4‑150 (MQ=255)
                                                                                                                                        ccttttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGc              >  1:1366818/4‑150 (MQ=255)
                                                                                                                                          ttttCGCGGGCTTGCCCGCGAAAGGGCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTccc          >  2:2998576/2‑150 (MQ=255)
                                                                                                                                           tttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGcttctt                  >  1:2025185/1‑143 (MQ=255)
                                                                                                                                           tttCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTcc           >  2:2641571/1‑150 (MQ=255)
                                                                                                                                               gcgGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgcgc                                                                                 >  1:1555748/1‑76 (MQ=255)
                                                                                                                                                 gggCCAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCGTTGTCCAGCGCCTCCAGCAGCCCTTGGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCTTGTTCAGCTTCTTCAGCTCCCGCGCt     >  1:5288407/1‑150 (MQ=255)
                                                                                                                                                  ggCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTg    >  2:4147160/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2069968‑2070264

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