Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,966,1270TC30.1% 103.4 / 67.4 83E22G (GAA→GGA) flhBflagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base C (9/16);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGG  >  NC_000913/1965988‑1966263
                                                                                                                                           |                                                                                                                                        
gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                                                                                                                          >  1:764427/1‑140 (MQ=255)
gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                                                                                                                          >  1:756266/1‑140 (MQ=255)
gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                                                                                                                          >  1:579002/1‑140 (MQ=255)
gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                                                                                                                          >  1:553276/1‑140 (MQ=255)
gtgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCAACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCtt                                                                                                                                          >  1:321877/1‑140 (MQ=255)
 tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                                                                                                         >  2:61840/1‑140 (MQ=255)
 tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                                                                                                         >  1:477849/1‑140 (MQ=255)
 tgATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTc                                                                                                                                         >  1:3967/1‑140 (MQ=255)
  gATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTcc                                                                                                                                        >  2:563925/1‑140 (MQ=255)
   aTCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTcccc                                                                                                                                       >  2:508320/1‑140 (MQ=255)
   aTCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTccc                                                                                                                                       >  2:514738/1‑140 (MQ=255)
     cAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCcgc                                                                                                                                     <  1:208212/140‑1 (MQ=255)
           gCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCtttt                                                                                                                               >  1:526716/1‑140 (MQ=255)
                 cAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGt                                                                                                                         <  2:77387/140‑1 (MQ=255)
                    cGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCgg                                                                                                                      <  2:1037911/140‑1 (MQ=255)
                     ggAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGt                                                                                                                     <  1:834641/140‑1 (MQ=255)
                      gagaGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTg                                                                                                                    <  2:538937/140‑1 (MQ=255)
                      gagaGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTg                                                                                                                    <  1:799073/140‑1 (MQ=255)
                         agCATGCCCGCCAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCATTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAgg                                                                                                                 >  2:224761/1‑140 (MQ=255)
                         agCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAgg                                                                                                                 >  1:224061/1‑140 (MQ=255)
                           cATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGgtg                                                                                                               >  1:562231/1‑140 (MQ=255)
                             tGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGtgtg                                                                                                             >  1:155584/1‑140 (MQ=255)
                              gCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                                                                                                            >  1:125880/1‑140 (MQ=255)
                              gCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTgg                                                                                                            >  1:42121/1‑140 (MQ=255)
                                 cGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTggggg                                                                                                         <  1:754079/140‑1 (MQ=255)
                                   aCAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCt                                                                                                       >  2:897910/1‑140 (MQ=255)
                                          aCGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGttt                                                                                                >  1:358574/1‑140 (MQ=255)
                                                cAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTgtcgt                                                                                          <  2:562231/140‑1 (MQ=255)
                                                 aGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTgtcgtc                                                                                         >  1:240482/1‑140 (MQ=255)
                                                  gCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGtcgtcg                                                                                        <  2:510285/140‑1 (MQ=255)
                                                  gCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGtcgtcg                                                                                        <  2:200761/140‑1 (MQ=255)
                                                    gACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGct                                                                                      <  2:813370/140‑1 (MQ=255)
                                                        cacCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCg                                                                                  >  2:120883/1‑140 (MQ=255)
                                                        cacCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCg                                                                                  >  2:503498/1‑140 (MQ=255)
                                                        cacCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCg                                                                                  >  1:972229/1‑140 (MQ=255)
                                                                aaCCAGATAACCCTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGacac                                                                          <  2:573190/140‑1 (MQ=255)
                                                                 aCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACg                                                                         >  1:670180/1‑140 (MQ=255)
                                                                 aCCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTACTGTTTTGTCGTCGCTCTCGTCAGACACg                                                                         >  1:769130/1‑140 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     <  2:224061/140‑1 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     <  2:756266/140‑1 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     <  1:514738/140‑1 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     >  1:884664/1‑140 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     <  2:274071/140‑1 (MQ=255)
                                                                     gATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGc                                                                     <  2:634407/140‑1 (MQ=255)
                                                                      aTAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGcc                                                                    <  2:389261/140‑1 (MQ=255)
                                                                      aTAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGcc                                                                    <  1:347287/140‑1 (MQ=255)
                                                                       tAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCa                                                                   <  2:309714/140‑1 (MQ=255)
                                                                          cacTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATc                                                                >  2:459915/1‑140 (MQ=255)
                                                                          cacTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATc                                                                <  1:897910/140‑1 (MQ=255)
                                                                           acTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATcc                                                               >  2:660208/1‑140 (MQ=255)
                                                                            cTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCt                                                              <  2:526716/140‑1 (MQ=255)
                                                                                cGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTCTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGaa                                                          <  2:477849/140‑1 (MQ=255)
                                                                                cGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGaa                                                          <  2:718130/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  2:579002/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  1:115853/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  1:120883/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  2:951081/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  2:843440/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        <  2:321877/140‑1 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        >  2:208233/1‑140 (MQ=255)
                                                                                  cccACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAag                                                        >  2:478311/1‑140 (MQ=255)
                                                                                    caccaGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAagag                                                      >  2:37605/1‑140 (MQ=255)
                                                                                     accaGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGt                                                     >  2:951780/1‑140 (MQ=255)
                                                                                      ccaGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTa                                                    <  1:57996/140‑1 (MQ=255)
                                                                                      ccaGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTa                                                    <  1:61840/140‑1 (MQ=255)
                                                                                      ccaGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTa                                                    <  1:223232/140‑1 (MQ=255)
                                                                                       caGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTAt                                                   <  1:442645/140‑1 (MQ=255)
                                                                                       caGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTAt                                                   <  2:972229/140‑1 (MQ=255)
                                                                                        aGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTAtt                                                  <  2:3967/140‑1 (MQ=255)
                                                                                          cAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTaa                                                >  2:113217/1‑140 (MQ=255)
                                                                                             aaTCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCg                                             <  1:173591/140‑1 (MQ=255)
                                                                                             aaTCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCg                                             <  1:660208/140‑1 (MQ=255)
                                                                                               tCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTg                                           >  1:1024246/1‑140 (MQ=255)
                                                                                                 agcagTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGaa                                         <  1:658581/140‑1 (MQ=255)
                                                                                                 agcagTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGaa                                         <  1:323286/140‑1 (MQ=255)
                                                                                                    agTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGCCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAatga                                      >  2:221015/1‑140 (MQ=255)
                                                                                                         ggTCAGTTCACGGGAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCa                                 <  1:508320/140‑1 (MQ=255)
                                                                                                                  aCGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCaaa                        >  1:1032730/1‑140 (MQ=255)
                                                                                                                    gggAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGc                      >  1:180983/1‑140 (MQ=255)
                                                                                                                     ggAACGCGGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGCCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCg                     <  1:221015/140‑1 (MQ=255)
                                                                                                                         cgcgGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCa                 >  1:819224/1‑140 (MQ=255)
                                                                                                                         cgcgGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTACTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGATTGATGCCAGAGCGCAAAGCGTTCa                 <  2:769130/140‑1 (MQ=255)
                                                                                                                          gcgGGATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCaa                <  1:985480/140‑1 (MQ=255)
                                                                                                                           cgGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAAt               <  2:670180/140‑1 (MQ=255)
                                                                                                                             ggATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATgg             <  2:553276/140‑1 (MQ=255)
                                                                                                                             ggATTTGCCCCTCTCCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATgg             <  2:455507/140‑1 (MQ=255)
                                                                                                                                    ccccTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGt      <  2:764427/140‑1 (MQ=255)
                                                                                                                                     cccTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTa     <  1:37605/140‑1 (MQ=255)
                                                                                                                                     cccTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTa     <  1:563925/140‑1 (MQ=255)
                                                                                                                                     cccTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTa     <  2:205402/140‑1 (MQ=255)
                                                                                                                                        tctcCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAgg  >  2:855521/1‑140 (MQ=255)
                                                                                                                                        tctTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAgg  >  2:419021/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GTGATCAAAATGCAGCCCAGCGGAGAGCATGCCCGACAATCGACGGGCCAGCGACACACCGCCAAACCAGATAACACTAACGCCCACCAGCAAAATCAGCAGTGAGGTCAGTTCACGGGAACGCGGGATTTGCCCCTCTTCCCGCGCTTTTTCTAGTCGGTGAGGTGTGGGGGCTTCTGTTTTGTCGTCGCTCTCGTCAGACACGTCGCCAATCCTGGAAAGAGTATTAAGCGTGAATGATGCCAGAGCGCAAAGCGTTCAATGGTTTGAGTAAGG  >  NC_000913/1965988‑1966263

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