Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,068,194 A→T 6.1% 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,068,1940AT 6.1% 230.9 / 6.4 82intergenic (‑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 (39/37);  new base T (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.91e-01

CAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068047‑2068342
                                                                                                                                                    |                                                                                                                                                    
cAACGGAAAGGAAAAGCAGCCCGAGCAGGAACATACGCAATGAAGAAATCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAg                                                                                                                                                    >  2:1310231/1‑150 (MQ=255)
          aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCTTGCCCGCGaaag                                                                                                                                          <  2:633289/150‑2 (MQ=255)
          aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  1:340903/150‑1 (MQ=255)
          aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  2:2150715/150‑1 (MQ=255)
          aaaaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACCGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGaaaa                                                                                                                                          <  1:528050/150‑1 (MQ=255)
            aaTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTTCGCGGGCTTGCCCGCGaaagg                                                                                                                                         <  1:1508438/150‑3 (MQ=255)
             aTGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGGCGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAggg                                                                                                                                       <  1:354437/150‑1 (MQ=255)
              tacAGCCACAGCAAGAACATACGCAATGAATAAATCATATTACGACCCACCCCAAACGGTCGCGGGGCAGTGTACGGGTGTGCACTGTCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGAAAGCCCGCGAAAAGGGc                                                                                                                                      <  2:850520/148‑1 (MQ=255)
                   ccccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTg                                                                                                                                 >  1:1656484/1‑150 (MQ=255)
                     ccAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCg                                                                                                                               >  1:1957115/1‑150 (MQ=255)
                          agtaaCATACGCAATTAATAAAGCATGGTACGACTCACACCAGACGGTGGCGGGGCAGTGTACGGGTGTGCAACGGCAAGCGCAACCCGACGTGCGGCACAGGTTGAGCCCTTTTGCTTGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAg                                                                                                                          <  1:1797552/147‑1 (MQ=255)
                          agGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAg                                                                                                                          >  2:2257477/1‑150 (MQ=255)
                           ggcaCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGACCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                         <  1:140788/147‑1 (MQ=255)
                            gAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCCTGGTGCAGGCCTGTTACGCCCGTTTCGCGGGCCAGCCCGCCAAAAGGGCCCGTGCGGCTAGcc                                                                                                                        >  2:140788/1‑150 (MQ=255)
                              aCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                         <  1:4380/147‑1 (MQ=255)
                              aCATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGc                                                                                                                         >  2:4380/1‑147 (MQ=255)
                                     cAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  1:1004015/150‑1 (MQ=255)
                                       aTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                             >  2:520764/1‑150 (MQ=255)
                                       aTGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATgctgct                                                                                                             >  1:701273/1‑150 (MQ=255)
                                             aaaGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCCGGCAAGCCCGCGAAAATGGCCCGTGCGCCTAGCCAATAGGCTGCTGATCCa                                                                                                       >  2:231968/1‑150 (MQ=255)
                                                      tACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGcc                                                                                                                        >  2:2192531/1‑124 (MQ=255)
                                                      tACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGcc                                                                                                                        <  1:2192531/124‑1 (MQ=255)
                                                      tACGACTCACACCAGACGGTCGCGGGCCAGTGTACGGGTGTGCGCGGGCTGGGGCAAAGCTTGGTGCAGGACTGTTATGGCCCTTTCGCGGGGAAGCCCGCCAAATGGAAACGTGCGGCTAGCCAATATGCTGCTGATCAAGCTCGATCg                                                                                              >  2:813706/1‑150 (MQ=255)
                                                        cGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGcc                                                                                            <  1:622327/150‑1 (MQ=255)
                                                            tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                        >  1:905698/1‑150 (MQ=255)
                                                            tCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGt                                                                                        >  1:1077311/1‑150 (MQ=255)
                                                             cacaCCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTc                                                                                       >  1:538316/1‑150 (MQ=255)
                                                              acacCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTcc                                                                                      >  2:939666/1‑150 (MQ=255)
                                                               ccccAGACGGTCGCGGGGCGGTGTACGGGTGTGCACTGGCTGGGCCAAAACTTGTTGAAGGCCTGTTACGGCCCTCCCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTACa                                                                                     <  2:624094/148‑1 (MQ=255)
                                                                acCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAg                                                                                    >  2:1763466/1‑150 (MQ=255)
                                                                  cAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctgc                                                                                                              <  2:1883486/122‑1 (MQ=255)
                                                                  cAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAggctgc                                                                                                              >  1:1883486/1‑122 (MQ=255)
                                                                  cAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgcg                                                                                  >  1:1739102/1‑150 (MQ=255)
                                                                  cAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGAAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAgcg                                                                                  >  2:1950594/1‑150 (MQ=255)
                                                                   aGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtg                                                                                                                  >  2:1644624/1‑117 (MQ=255)
                                                                   aGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtg                                                                                                                  <  1:1644624/117‑1 (MQ=255)
                                                                        gTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCca                                                                            <  2:1739102/150‑1 (MQ=255)
                                                                        gTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCca                                                                            <  1:219066/150‑1 (MQ=255)
                                                                         tcgaggggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCcag                                                                           <  2:2209235/146‑1 (MQ=255)
                                                                         tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGc                                                                                                                                                   <  1:1027052/78‑1 (MQ=255)
                                                                         tCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGc                                                                                                                                                   >  2:1027052/1‑78 (MQ=255)
                                                                             ggggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGa                                                                                                                                             >  1:1170543/1‑80 (MQ=255)
                                                                             ggggCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGa                                                                                                                                             <  2:1170543/80‑1 (MQ=255)
                                                                                    tgtACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                <  1:2257477/150‑1 (MQ=255)
                                                                                    tgtACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                <  1:63709/150‑1 (MQ=255)
                                                                                    tggaCGGGTGTGCACTGGCTGGGGCAACCCTTGCTGCAGTTCTGTTACGGCCCTTTCGCGGGCAAGCCCGAGGGAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgc                                                                <  1:520764/147‑1 (MQ=255)
                                                                                      tACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTgcgc                                                              >  2:821772/1‑150 (MQ=255)
                                                                                         gggTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATACTGCTGATCCAGCTCGATCGcc                                                                                            >  2:789752/1‑117 (MQ=255)
                                                                                         gggTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATACTGCTGATCCAGCTCGATCGcc                                                                                            <  1:789752/117‑1 (MQ=255)
                                                                                           gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               <  1:1388942/96‑1 (MQ=255)
                                                                                           gtgtGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATAtgctg                                                                                                               >  2:1388942/1‑96 (MQ=255)
                                                                                                 aCTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGACGCTCCAGAGCTGTCGCCGTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCtt                                                   >  1:1687484/1‑150 (MQ=255)
                                                                                                  ctggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  <  2:120848/150‑1 (MQ=255)
                                                                                                  ctggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  <  2:701273/150‑1 (MQ=255)
                                                                                                    ggctggGGCAAACCTTGGTGCAGGCCTGTTCCGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCCCCAGCCCTTCGCGCACTTTCTTCTTGgt                                                <  2:514098/150‑1 (MQ=255)
                                                                                                    ggctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgt                                                <  1:821772/150‑1 (MQ=255)
                                                                                                     gctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCa                                                                                     >  2:947041/1‑112 (MQ=255)
                                                                                                     gctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCa                                                                                     <  1:947041/112‑1 (MQ=255)
                                                                                                      ctggGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGgtgt                                              <  2:905698/150‑1 (MQ=255)
                                                                                                         gggCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCt                                           <  1:925833/150‑1 (MQ=255)
                                                                                                          ggCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGGGACTGGAGAGGCGCTTTGCGGACTTTCATGTTGGTgttgtt                                          >  1:1072528/1‑150 (MQ=255)
                                                                                                           gaaaaCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTCCTtg                                         <  1:1950594/148‑1 (MQ=255)
                                                                                                            cAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  <  2:587011/140‑1 (MQ=255)
                                                                                                            cAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTg                                                  >  1:587011/1‑140 (MQ=255)
                                                                                                               aCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCagcagc                                                                       <  1:1328805/116‑1 (MQ=255)
                                                                                                               aCCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCagcagc                                                                       >  2:1328805/1‑116 (MQ=255)
                                                                                                                 cTTGGTGCACCCCTCTTCAGGCCATTTAGCGGCTAAGCGCACTGAAAGAGGCCGTGCGGCTAGCCAATATGCTGCTGAACCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGcc                                   >  2:1092588/1‑150 (MQ=255)
                                                                                                                    ggTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGc                                >  1:26860/1‑150 (MQ=255)
                                                                                                                     gTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               <  1:795610/86‑1 (MQ=255)
                                                                                                                     gTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               >  2:795610/1‑86 (MQ=255)
                                                                                                                      tGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCAt                              >  1:216435/1‑150 (MQ=255)
                                                                                                                            ccTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAACAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAg                        <  2:1077311/150‑1 (MQ=255)
                                                                                                                             cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCtt                                          >  1:1763477/1‑131 (MQ=255)
                                                                                                                             cTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCtt                                          <  2:1763477/131‑1 (MQ=255)
                                                                                                                             cTGTTACGGCCCTTTCGCGGGCAAGCCCGAGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATGCAGCTCGATCGCCTTGTCCAGCGCCTACAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGCTCGCCCGCGTGTTCatc                       >  2:390788/1‑148 (MQ=255)
                                                                                                                              tGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTGTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGct                      >  1:1642215/1‑150 (MQ=255)
                                                                                                                                      cccTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGc              >  1:1865971/1‑150 (MQ=255)
                                                                                                                                         ttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               >  2:1818968/2‑66 (MQ=25)
                                                                                                                                         ttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATc                                                                                               <  1:1818968/65‑1 (MQ=25)
                                                                                                                                         ttttCGCGGGCTTGCCCGCGAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTcc           >  2:1511495/2‑150 (MQ=255)
                                                                                                                                                gggCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCAAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCt     >  2:190597/1‑150 (MQ=255)
                                                                                                                                                   cAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCACTTTGCTGCTTATCCAGCTCGATCGCCTTGCCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTGAGCTTGGTGTTCTTGTGCGCCAGCCTGTTCAGCTTCTTCAGCTCCCGCGCTGCg  >  1:187078/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAACGGCAAGAAAATGCAGCCCCAGCAGGAACATACGCAATGAATAAAGCATGGTACGACTCACACCAGACGGTCGCGGGGCAGTGTACGGGTGTGCACTGGCTGGGGCAAACCTTGGTGCAGGCCTGTTACGGCCCTTTCGCGGGCAAGCCCGCGAAAAGGGCCCGTGCGGCTAGCCAATATGCTGCTGATCCAGCTCGATCGCCTTGTCCAGCGCCTCCAGCAGCCCTTTGCGCACTTTCAGCTTGGTGTTCTTGTGCGCCAGCATGTTCAGCTTCTTCAGCTCCCGCGCTGCG  >  NC_002947/2068047‑2068342

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