Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,474 T→C 5.7% intergenic (+109/+35) galU → / ← hns glucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,292,4740TC 5.7% 1016.2 / 36.4 322intergenic (+109/+35)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base T (126/168);  new base C (14/4);  total (148/174)
Fisher's exact test for biased strand distribution p-value = 5.82e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 325 total aligned reads displayed.

TTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTC  >  NC_000913/1292325‑1292622
                                                                                                                                                     |                                                                                                                                                    
ttAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCt                                                                                                                                                      >  1:3293463/1‑150 (MQ=255)
 tAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                     >  1:925872/1‑150 (MQ=255)
 tAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                     >  2:759586/1‑150 (MQ=255)
  aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTgg                                                                                                                                                    >  1:1600372/1‑150 (MQ=255)
   aaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGc                                                                                                                                                   >  2:5351162/1‑150 (MQ=255)
    aGCCTGGCTTGAAGAAGAGATGGGAATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGAAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCg                                                                                                                                                  <  1:1560432/150‑1 (MQ=255)
    aGCCGGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCg                                                                                                                                                  <  2:2122627/150‑1 (MQ=255)
     gCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGCCAATAAAAAATCCCGCCGCTGGCgg                                                                                                                                                 >  2:3005669/1‑150 (MQ=255)
     gCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCgg                                                                                                                                                 <  2:4279913/150‑1 (MQ=255)
       cTGGCTTGCAGCAGCGCGGGGCATGAAGAAGTAACATCCGTATCGGTGTTATCCACGCACCGGCGGTGCGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGa                                                                                                                                               <  2:4940457/150‑1 (MQ=255)
        tGGCTTGAAGAAGAGAGGGGCATTAAGAAGTAACATCCGTATCGGTGGTATCCACGAAACGGCGTTGAGCAAACGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAt                                                                                                                                              <  2:4029712/150‑1 (MQ=255)
        gggCTTGAAGATGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAt                                                                                                                                              <  2:2856094/149‑1 (MQ=255)
         ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtt                                                                                                                                             <  1:4395419/150‑1 (MQ=255)
          gCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAttt                                                                                                                                            <  2:2973112/150‑1 (MQ=255)
            ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAttttt                                                                                                                                          <  2:1947706/150‑2 (MQ=255)
            ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTa                                                                                                                                          >  1:4262118/1‑150 (MQ=255)
             tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttttt                                                                                                                                         <  2:2783292/150‑3 (MQ=255)
              gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                                                                                        <  1:244335/150‑1 (MQ=255)
                    gagaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTg                                                                                                                                  >  1:3927701/1‑150 (MQ=255)
                     agaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGc                                                                                                                                 >  1:1193039/1‑150 (MQ=255)
                       aTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGGTTATTGTTATTAAATTGTCTTAAACCGGACAATAAAAAATACCGCCGCTGGCGGGATTTTTAGCAAGTGCaa                                                                                                                               >  1:1225385/1‑150 (MQ=255)
                          ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAAGTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                                            >  2:1251858/1‑150 (MQ=255)
                              ttAAGCCGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTGATAGCTTATTCTTATTACATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCGGGCGGGATTTTAAGCAAGTGCAATCTACaa                                                                                                                        <  2:4785431/150‑1 (MQ=255)
                              ttAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaa                                                                                                                        >  1:673992/1‑150 (MQ=255)
                              ttAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaa                                                                                                                        >  1:728687/1‑150 (MQ=255)
                              ttAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaa                                                                                                                        <  1:1630949/150‑1 (MQ=255)
                               tAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAAACGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaa                                                                                                                       >  1:1778888/1‑150 (MQ=255)
                                aagaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                      <  1:3428/150‑1 (MQ=255)
                                aagaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                      >  1:2908051/1‑150 (MQ=255)
                                    agTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt                                                                                                                  <  1:3633105/150‑1 (MQ=255)
                                          aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                            >  1:2235609/1‑150 (MQ=255)
                                          aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCGGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                            >  1:4430215/1‑150 (MQ=255)
                                            ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCCCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                                                                                          >  1:4567610/1‑150 (MQ=255)
                                              gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                                        <  2:3955405/150‑1 (MQ=255)
                                              gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                                        <  1:5115297/150‑1 (MQ=255)
                                               tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                       <  1:2172733/150‑1 (MQ=255)
                                               tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                       <  1:5460841/150‑1 (MQ=255)
                                                aTCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCa                                                                                                      <  2:1014797/150‑1 (MQ=255)
                                                 tCGGTGTTATCCACGAATCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAg                                                                                                     <  1:137779/150‑1 (MQ=255)
                                                  cGGAGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                                    <  2:2183420/150‑1 (MQ=255)
                                                     tgtTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaaa                                                                                                 <  1:1279410/150‑1 (MQ=255)
                                                          tCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTTCTTGATCAGGAAAtcgtc                                                                                            >  1:5269955/1‑150 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATCGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCGGGCGGGATTTTAAGCAAGGGCAATCTACAAAAGATTATTGCTTGATCAGGAACTCGTCGa                                                                                          <  2:2389460/150‑1 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGGTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATACCGCCGCTGGCGGGATTTTAAGCAAGTGCAAGCAAAAAAAGATTATTTCTTGTTCAGGAAAAAGTCGa                                                                                          >  1:2258486/1‑150 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTTTCTTAAACCGGACAATAAAAAATACCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAATATTATTTCTTGATCAGGAAATCGTCGa                                                                                          >  1:2662193/1‑150 (MQ=255)
                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAAACGGACAATAAAAAATCCCGCCGCGGGCGGGATTTTAAGCAAGTGCAATATACAAAAGATTATTGCTTGATCAGGAAAAAGTCGa                                                                                          >  1:4696344/1‑150 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                         <  1:4874563/150‑1 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATAAGGAAATCGTCGAg                                                                                         >  1:3371352/1‑150 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCCTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCACGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                         <  2:3124233/150‑1 (MQ=255)
                                                             aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCGTATTATTATTAAATTGTATTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAGCGTCGAg                                                                                         >  1:10378/1‑150 (MQ=255)
                                                              cGAAACGGCGTTGAGCAATCGACGCCGTTTTTGTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATGTTAAGCACGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                                                                                        >  2:335810/1‑150 (MQ=255)
                                                                   cGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                                   >  2:3177395/1‑150 (MQ=255)
                                                                     gCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTATTATTAAATTGTCTTAAAACGGACAAAAAAAAAACACGCCCCGGGCGGGATTTTAATCAATTGCAATCTACAAAAGATTATTTCTTGATAAGGAAAACGTCGAGGGATTTaa                                                                                 >  1:2763090/1‑149 (MQ=255)
                                                                        ttGAGCAATCGACGCCGTTTTTTTTTAGCTTATTCTTATTAAAATGTCTTAAAACGGACAAAAAAAAAAACCGCCGCTGGCGGGGTTTTAAGCAAGTGCAATCAACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCtt                                                                              >  1:632947/1‑150 (MQ=255)
                                                                         tGAGCACTCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTg                                                                             <  1:4323222/150‑1 (MQ=255)
                                                                          gAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                                            <  1:3663304/150‑1 (MQ=255)
                                                                               aTCGACGCAGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                       <  2:3474012/150‑1 (MQ=255)
                                                                                         tttttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                                             <  1:1493343/150‑1 (MQ=255)
                                                                                             tttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGa                                                         <  2:2925184/150‑1 (MQ=255)
                                                                                             tttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGa                                                         <  1:1971705/150‑1 (MQ=255)
                                                                                                  aGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACa                                                    <  2:5137898/150‑1 (MQ=255)
                                                                                                     ttattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                                                 <  2:2235609/150‑1 (MQ=255)
                                                                                                     ttattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                                                 >  1:1071153/1‑150 (MQ=255)
                                                                                                        ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGa                                              <  2:728687/150‑1 (MQ=255)
                                                                                                        ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGa                                              <  2:5212052/150‑1 (MQ=255)
                                                                                                         tcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAg                                             >  1:4033275/1‑150 (MQ=255)
                                                                                                          cttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGt                                            <  1:2064821/150‑1 (MQ=255)
                                                                                                          cttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGt                                            <  1:2823727/150‑1 (MQ=255)
                                                                                                           ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTa                                           >  2:938028/1‑150 (MQ=255)
                                                                                                           ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTa                                           >  2:4307586/1‑150 (MQ=255)
                                                                                                            tattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTAc                                          >  2:4493853/1‑150 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                                        >  1:3109901/1‑150 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                                        >  2:4139993/1‑150 (MQ=255)
                                                                                                              ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                                        >  1:3376145/1‑150 (MQ=255)
                                                                                                                aaaTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGcc                                      >  2:167171/1‑150 (MQ=255)
                                                                                                                 aaTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGTTTACAGCTGGAGTACGGCCt                                     >  1:5170095/1‑150 (MQ=255)
                                                                                                                 aaTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCt                                     >  1:3577381/1‑150 (MQ=255)
                                                                                                                  aTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt                                    <  1:5033871/150‑1 (MQ=255)
                                                                                                                    tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                                  <  1:412255/150‑1 (MQ=255)
                                                                                                                    tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                                  <  2:2593785/150‑1 (MQ=255)
                                                                                                                     gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                                 <  2:1519220/150‑1 (MQ=255)
                                                                                                                          aaaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc                            <  2:5164716/150‑1 (MQ=255)
                                                                                                                             ccGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                         <  2:200906/150‑1 (MQ=255)
                                                                                                                               ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGt                       <  1:1732150/150‑1 (MQ=255)
                                                                                                                                 aCAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGttt                     <  2:3781889/150‑1 (MQ=255)
                                                                                                                                 aCAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGttt                     <  1:2594781/150‑1 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                   >  1:5483603/1‑150 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                   >  1:2739485/1‑150 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                   >  1:3049357/1‑150 (MQ=255)
                                                                                                                                   aaTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                   <  2:1604130/150‑1 (MQ=255)
                                                                                                                                      aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGtt                >  1:2541125/1‑150 (MQ=255)
                                                                                                                                      aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAACTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGtt                <  2:831095/150‑1 (MQ=255)
                                                                                                                                       taaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGttt               >  1:2571244/2‑150 (MQ=255)
                                                                                                                                       aaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGttt               <  2:4622961/150‑1 (MQ=255)
                                                                                                                                       aaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATAAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTGTTTGTTTACAGCTGGAGTACGGACTGGGCCAGGCCAGGTTTTAGttt               >  1:4387847/1‑150 (MQ=255)
                                                                                                                                        aaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc              >  1:3694468/1‑150 (MQ=255)
                                                                                                                                        aaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc              >  2:1726200/1‑150 (MQ=255)
                                                                                                                                           aTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGcc           <  1:974516/150‑1 (MQ=255)
                                                                                                                                           aTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGcc           >  1:4607756/1‑150 (MQ=255)
                                                                                                                                                  cgcTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCg    <  1:4651713/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTC  >  NC_000913/1292325‑1292622

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