Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,256,369 A→G 5.8% L33P (CTG→CCG)  plsB ← glycerol‑3‑phosphate O‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,256,3690AG 5.8% 421.4 / 11.3 140L33P (CTG→CCG) plsBglycerol‑3‑phosphate O‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (66/64);  new base G (4/4);  total (71/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 143 total aligned reads displayed.

CATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAGGATAGC  >  NC_000913/4256220‑4256516
                                                                                                                                                     |                                                                                                                                                   
caTAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTACCGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCa                                                                                                                                                     >  1:653007/1‑150 (MQ=255)
caTAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCa                                                                                                                                                     >  2:1311116/1‑150 (MQ=255)
   aGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGtt                                                                                                                                                  >  1:1446608/1‑150 (MQ=255)
    gcgcGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGATc                                                                                                                                                 >  2:421067/1‑150 (MQ=255)
     cgcgGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTcc                                                                                                                                                >  2:1181155/1‑150 (MQ=255)
     cgcgGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTcc                                                                                                                                                <  1:1443181/150‑1 (MQ=255)
       cgGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGCCACTGGGCGGGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCgg                                                                                                                                              >  2:1150026/1‑150 (MQ=255)
         gCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCgggg                                                                                                                                            >  1:1662457/1‑150 (MQ=255)
          cAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCACCGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGc                                                                                                                                           >  1:320410/1‑150 (MQ=255)
          cAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGc                                                                                                                                           >  1:769697/1‑150 (MQ=255)
           aGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCa                                                                                                                                          <  1:1275030/150‑1 (MQ=255)
           aGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCa                                                                                                                                          <  1:158404/150‑1 (MQ=255)
           aGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCa                                                                                                                                          <  2:569054/150‑1 (MQ=255)
            gTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCCTGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAg                                                                                                                                         <  1:1190920/150‑1 (MQ=255)
              aGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGa                                                                                                                                       <  1:479262/150‑1 (MQ=255)
              aGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGa                                                                                                                                       >  2:579113/1‑150 (MQ=255)
               gCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCGGAACTGGGGCTGGAt                                                                                                                                      <  2:816917/150‑1 (MQ=255)
               gCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGAt                                                                                                                                      >  1:578531/1‑150 (MQ=255)
                 gTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCTGCCCCGGAACTGGGGCTGGAtac                                                                                                                                    <  2:934084/150‑3 (MQ=255)
                    ccGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCTGCCCCGGAACTGGGGCTGGAtacctc                                                                                                                                 <  1:724179/150‑6 (MQ=255)
                     cGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCg                                                                                                                                >  2:44253/1‑150 (MQ=255)
                         gATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAAt                                                                                                                            >  2:710/1‑150 (MQ=255)
                         gATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAAt                                                                                                                            <  1:1268158/150‑1 (MQ=255)
                               cAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACtt                                                                                                                      >  2:419592/1‑150 (MQ=255)
                                aGCGGCGCTAACGGGTCAGGCAATTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTg                                                                                                                     <  2:219741/150‑1 (MQ=255)
                                   ggCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCtt                                                                                                                  <  2:800927/150‑1 (MQ=255)
                                       ctAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCtttttt                                                                                                              >  1:1100420/1‑149 (MQ=255)
                                       ctAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTa                                                                                                              >  1:368192/1‑150 (MQ=255)
                                       ctAACGGGTCAGGCAAGTCATGAGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTa                                                                                                              <  1:1057391/150‑1 (MQ=255)
                                       ctAACGGGACAGGCAAGTCATGTGCCAGACACTGGGCCCGCAACGTCAGCAAAACTGCTTTCGAGTTGTACGGGTTCACGTGCATTATTGGATGAGAGGTTTCCTGCCCCAGTTCCGGGGCCGGATCTGCGGGAATAGACTTGcagtgta                                                                                                              >  2:328613/1‑144 (MQ=255)
                                       ctAACGCGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGTAATAGACTTGCTTTTTa                                                                                                              <  1:895392/150‑1 (MQ=255)
                                        tAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTAc                                                                                                             <  2:1311891/150‑1 (MQ=255)
                                        tAACGGGGCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTAc                                                                                                             <  1:1139844/150‑1 (MQ=255)
                                        aaaCGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCGTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTAc                                                                                                             <  1:838573/149‑1 (MQ=255)
                                          aCGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCa                                                                                                           >  1:61567/1‑150 (MQ=255)
                                           cGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAg                                                                                                          >  1:1083189/1‑150 (MQ=255)
                                                  ggCAAGTCATGTGCCAGACACTGGGCGCGCAACGTAAGCTAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCGCAGTTCCGGGGCAGGTTCTGCCGGAATAGACTGGCTTTTTAGCAGGATGGtt                                                                                                   >  1:1132980/1‑150 (MQ=255)
                                                   gCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTa                                                                                                  >  1:983559/1‑150 (MQ=255)
                                                   gCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTa                                                                                                  >  1:870584/1‑150 (MQ=255)
                                                    cAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTaa                                                                                                 >  2:91375/1‑150 (MQ=255)
                                                     aaGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAAt                                                                                                >  2:833299/1‑150 (MQ=255)
                                                      aGTCATGGGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATg                                                                                               <  2:1553104/150‑1 (MQ=255)
                                                       gTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGGGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATgg                                                                                              >  1:784393/1‑150 (MQ=255)
                                                       gTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATgg                                                                                              >  1:546835/1‑150 (MQ=255)
                                                        tCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGt                                                                                             <  2:1446608/150‑1 (MQ=255)
                                                         cATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTa                                                                                            <  1:857556/150‑1 (MQ=255)
                                                          aTGTTCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACTAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTaa                                                                                           <  1:1563032/150‑1 (MQ=255)
                                                              gCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTc                                                                                       >  1:327803/1‑150 (MQ=255)
                                                                cAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGATGTATCCAGCCCCTGTTCCGGGGCCGGATCTGCCGGAATAGAATTGCTTTTTACCAGGATGCTTAATTGTAAATTCAg                                                                                     >  2:293189/1‑150 (MQ=255)
                                                                 aGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTTGACGAGAGGTATCCATCCCCAGTTCCGGGGCAGGATCTGCCGTAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGt                                                                                    <  1:857283/150‑1 (MQ=255)
                                                                    cacTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAAt                                                                                 <  2:962224/150‑1 (MQ=255)
                                                                     acTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAAtt                                                                                >  2:1474025/1‑150 (MQ=255)
                                                                      cTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAAttt                                                                               <  2:600960/150‑1 (MQ=255)
                                                                       tGGGCGCGCAACGGCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGATAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTg                                                                              <  1:1469396/150‑1 (MQ=255)
                                                                        gggCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTgt                                                                             <  1:1356885/150‑1 (MQ=255)
                                                                          gcgcgcAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTgtag                                                                           <  2:870584/150‑1 (MQ=255)
                                                                             cgcAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTaa                                                                        >  2:1526989/1‑150 (MQ=255)
                                                                               cAACGTCAGCAATTCTGCGTTCGCGTTGTACGGTAAACCGTACATAATTGGACGAGAGGGTTACATCCACTGATACGGGGCAGTATCTGACTCAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAAt                                                                      <  1:293189/150‑1 (MQ=255)
                                                                               cAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAAt                                                                      <  1:1382723/150‑1 (MQ=255)
                                                                                     cAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTgg                                                                <  1:579113/150‑1 (MQ=255)
                                                                                     cAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTgg                                                                <  1:1001125/150‑1 (MQ=255)
                                                                                     cAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTAAGTAATATGTAGTAAATTAGTgg                                                                >  1:1686639/1‑150 (MQ=255)
                                                                                     cAGCAAATCAGCTTTCGAGTTTTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTgg                                                                <  2:868986/150‑1 (MQ=255)
                                                                                       gCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGcc                                                              <  1:421067/150‑1 (MQ=255)
                                                                                         aaaTCTGCTTTCGAGTTGTTCGGTAAAAAGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGGAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAg                                                            <  1:914258/150‑1 (MQ=255)
                                                                                         aaaTCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAg                                                            <  1:1195177/150‑1 (MQ=255)
                                                                                            tCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTGAATTCGTGGCCAGCCg                                                         <  1:1691691/150‑1 (MQ=255)
                                                                                              tGCTTTCGAGTTGTACTGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGa                                                       <  1:446343/150‑1 (MQ=255)
                                                                                               gCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGAc                                                      <  2:1080583/150‑1 (MQ=255)
                                                                                                cTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACa                                                     <  2:769697/150‑1 (MQ=255)
                                                                                                cTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACa                                                     >  1:156965/1‑150 (MQ=255)
                                                                                                cTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACa                                                     >  2:1439007/1‑150 (MQ=255)
                                                                                                cTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACa                                                     >  1:1475479/1‑150 (MQ=255)
                                                                                                    cGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATaaa                                                 >  1:1305807/1‑150 (MQ=255)
                                                                                                      aGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACg                                               >  2:1463335/1‑150 (MQ=255)
                                                                                                              ggTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAg                                       >  1:250792/1‑150 (MQ=255)
                                                                                                                      gTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGt                               >  2:525883/1‑150 (MQ=255)
                                                                                                                      gTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGt                               >  1:1347150/1‑150 (MQ=255)
                                                                                                                       tACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGtt                              <  2:320410/150‑1 (MQ=255)
                                                                                                                          aTAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTtaat                           <  2:983559/150‑1 (MQ=255)
                                                                                                                              ttGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATg                       <  1:1319295/150‑1 (MQ=255)
                                                                                                                                ggACGAGAGGTATCCATCCCCAGTTCCGGGGCAGGAGCTGCCGGAATAGACTTGCTTGTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCa                     <  2:1323387/150‑1 (MQ=255)
                                                                                                                                ggACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCa                     >  1:388767/1‑150 (MQ=255)
                                                                                                                                ggACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCa                     <  2:1128765/150‑1 (MQ=255)
                                                                                                                                  aCGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCaaa                   <  1:54400/150‑1 (MQ=255)
                                                                                                                                  aCGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCaaa                   <  2:1675302/150‑1 (MQ=255)
                                                                                                                                   cGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGTACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAAt                  >  2:19783/1‑150 (MQ=255)
                                                                                                                                    gagaGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATg                 <  1:1181155/150‑1 (MQ=255)
                                                                                                                                     agagGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATgc                >  1:952302/1‑150 (MQ=255)
                                                                                                                                        ggTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGcgcg             <  1:1662614/150‑1 (MQ=255)
                                                                                                                                         gTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCgg            >  2:558852/1‑150 (MQ=255)
                                                                                                                                         agaTCCTGCCCCGGAACTGGGGCTGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCgg            >  1:23766/3‑150 (MQ=255)
                                                                                                                                          gaTCCTGCCCCGGAACTGGGGCTGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGc           >  2:1369851/2‑150 (MQ=255)
                                                                                                                                           aTCCTGCCCCGGAACTGGGGCTGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCa          >  1:857879/1‑150 (MQ=255)
                                                                                                                                           aTCCTGCCCCGGAACTGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCa          >  2:18375/1‑150 (MQ=255)
                                                                                                                                             ccAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAg        >  2:624059/1‑150 (MQ=255)
                                                                                                                                             ccAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAg        >  2:1332848/1‑150 (MQ=255)
                                                                                                                                               aGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAGGa      >  1:688352/1‑150 (MQ=255)
                                                                                                                                                   ccAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATTCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAGGATAgt  <  1:1065224/150‑2 (MQ=255)
                                                                                                                                                   ccAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAGGATAGc  <  2:543364/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CATAGCGCGGCAGTAGCGTGCCGTCGATTTCCAGCGGCTCTAACGGGTCAGGCAAGTCATGTGCCAGACACTGGGCGCGCAACGTCAGCAAATCTGCTTTCGAGTTGTACGGTAAAACGTACATAATTGGACGAGAGGTATCCAGCCCCAGTTCCGGGGCAGGATCTGCCGGAATAGACTTGCTTTTTACCAGGATGCTTAATGGTAAATTCAGTAATTTGTAGTAAATTCGTGGCCAGCCGGACATAAACGATGTAAAGCCTCTGGTTAATAATGCAAATGCGCGGCAAGGATAGC  >  NC_000913/4256220‑4256516

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