Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,474 T→C 6.6% 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 6.6% 441.6 / 15.7 150intergenic (+109/+35)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base T (72/68);  new base C (5/5);  total (77/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-02

Only 100 of 151 total aligned reads displayed.

TAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTC  >  NC_000913/1292326‑1292622
                                                                                                                                                    |                                                                                                                                                    
 aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTgg                                                                                                                                                    >  1:428948/1‑150 (MQ=255)
  aaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGc                                                                                                                                                   <  1:401491/150‑1 (MQ=255)
    gCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGCCGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCgg                                                                                                                                                 >  2:1285691/1‑150 (MQ=255)
     ccTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCggg                                                                                                                                                >  1:430351/1‑150 (MQ=255)
        ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtt                                                                                                                                             <  1:953756/150‑1 (MQ=255)
         gCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtt                                                                                                                                             <  2:23418/149‑1 (MQ=255)
         gCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAttt                                                                                                                                            <  2:1346166/150‑1 (MQ=255)
          cTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtttt                                                                                                                                           <  2:1256471/150‑1 (MQ=255)
           ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTa                                                                                                                                          <  2:1477680/150‑1 (MQ=255)
            tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                         >  1:653787/1‑150 (MQ=255)
            tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                         >  1:475779/1‑150 (MQ=255)
            tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttttt                                                                                                                                         <  2:6741/150‑3 (MQ=255)
            tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttttt                                                                                                                                         <  2:1187308/150‑3 (MQ=255)
             gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCGGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                                                                                        >  1:518325/1‑150 (MQ=255)
             gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                                                                                        >  2:1427267/1‑150 (MQ=255)
              aagaagAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                                                                                                       <  1:1281926/150‑1 (MQ=255)
              aagaagAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                                                                                                       <  2:1341798/150‑1 (MQ=255)
               agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                                                                                                      <  1:691058/150‑1 (MQ=255)
                 aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAg                                                                                                                                    <  1:675445/150‑1 (MQ=255)
                 aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAg                                                                                                                                    <  2:757518/150‑1 (MQ=255)
                  agagATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                                                                                                   >  1:525696/1‑150 (MQ=255)
                         ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                                            >  2:672298/1‑150 (MQ=255)
                         ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                                            >  1:1473415/1‑150 (MQ=255)
                         ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                                            <  1:984283/150‑1 (MQ=255)
                          gCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTa                                                                                                                           <  1:1003112/150‑1 (MQ=255)
                          gCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTa                                                                                                                           >  2:902623/1‑150 (MQ=255)
                            aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACa                                                                                                                         <  1:1554941/150‑1 (MQ=255)
                               aagaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                      >  2:897704/1‑150 (MQ=255)
                                   agTAACATCAGTATCTGTGTTATCCACGAATCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTATATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt                                                                                                                  <  1:1137367/150‑1 (MQ=255)
                                        cATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCACGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTga                                                                                                             >  1:1620815/1‑149 (MQ=255)
                                         aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                            <  1:1250865/150‑1 (MQ=255)
                                           ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                                                                                          <  2:390405/150‑1 (MQ=255)
                                           ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                                                                                          <  2:428948/150‑1 (MQ=255)
                                             gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                                        <  1:82470/150‑1 (MQ=255)
                                              tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                       <  2:1398315/150‑1 (MQ=255)
                                              tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                                                                                                       <  2:210280/150‑1 (MQ=255)
                                                   gtgtTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaa                                                                                                  >  2:1326586/1‑150 (MQ=255)
                                                   gtgtTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaa                                                                                                  >  2:1247290/1‑150 (MQ=255)
                                                        aTCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                                    <  1:1410507/143‑1 (MQ=255)
                                                        aTCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                                                                                             >  2:837705/1‑150 (MQ=255)
                                                        aTCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                                                                                             >  1:1697149/1‑150 (MQ=255)
                                                           cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                          >  2:796767/1‑150 (MQ=255)
                                                            aaaaaaCGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                         >  1:60755/4‑150 (MQ=255)
                                                            aaaaaaCGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATAGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                         >  2:1114982/4‑150 (MQ=255)
                                                             cGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                                                                                        <  2:1319626/150‑1 (MQ=255)
                                                               aaaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                      <  1:1321771/150‑1 (MQ=255)
                                                               aaaCGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGCCAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                      >  2:1223208/1‑150 (MQ=255)
                                                                  cGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAttt                                                                                   >  2:1401903/1‑150 (MQ=255)
                                                                     cGTTGATCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAcc                                                                                <  1:999646/150‑1 (MQ=255)
                                                                      gTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCt                                                                               >  1:1014561/1‑150 (MQ=255)
                                                                           tcAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc                                                                          <  1:1241575/149‑1 (MQ=255)
                                                                              aTCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                       >  1:883534/1‑150 (MQ=255)
                                                                                cGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCa                                                                     <  1:25621/150‑1 (MQ=255)
                                                                                   cGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                                                  <  2:653787/150‑1 (MQ=255)
                                                                                         ttttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCttttt                                                            <  2:538224/150‑1 (MQ=255)
                                                                                           ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTg                                                          <  1:1601971/150‑1 (MQ=255)
                                                                                           ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTg                                                          <  2:141596/150‑1 (MQ=255)
                                                                                           ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTg                                                          <  1:541751/150‑1 (MQ=255)
                                                                                             ttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                                        <  1:1278439/150‑1 (MQ=255)
                                                                                              tataGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAtt                                                       >  2:650741/1‑150 (MQ=255)
                                                                                              tataGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAtt                                                       >  2:824618/1‑150 (MQ=255)
                                                                                              tataGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAtt                                                       <  2:1061848/150‑1 (MQ=255)
                                                                                              tataGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAtt                                                       <  1:1116179/150‑1 (MQ=255)
                                                                                                taGCTTATTCTTATTAAATAGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTAc                                                     <  1:1114982/150‑1 (MQ=255)
                                                                                                      attcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTgg                                               <  2:883534/150‑1 (MQ=255)
                                                                                                      attcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTgg                                               <  1:837705/150‑1 (MQ=255)
                                                                                                       ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGa                                              <  1:952549/150‑1 (MQ=255)
                                                                                                         cttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGTTTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGt                                            >  2:152234/1‑150 (MQ=255)
                                                                                                           tattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTAc                                          <  1:91466/150‑1 (MQ=255)
                                                                                                              taAATTTTCTTAAACCGGACAATAAAAAATCCCGGCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTTATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGc                                       <  2:884048/150‑1 (MQ=255)
                                                                                                              taAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGc                                       <  1:1478728/150‑1 (MQ=255)
                                                                                                                  ttGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTg                                   <  2:1661132/150‑1 (MQ=255)
                                                                                                                   tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                                  >  1:632116/1‑150 (MQ=255)
                                                                                                                    ttCTTAAACCTGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACTAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                                 <  1:752882/149‑1 (MQ=255)
                                                                                                                    gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                                 <  1:408294/150‑1 (MQ=255)
                                                                                                                      cTTAAACCGGACAATAGAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCa                               >  1:65580/1‑150 (MQ=255)
                                                                                                                       ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAg                              >  1:1390355/1‑150 (MQ=255)
                                                                                                                        tAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTGATCCATTGCTTTTTTGATTGCAGCTGGAGTACGGCCTTGGGCAGt                             >  1:308280/1‑150 (MQ=255)
                                                                                                                        tAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGt                             >  1:965804/1‑150 (MQ=255)
                                                                                                                        tAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGt                             >  1:179811/1‑150 (MQ=255)
                                                                                                                        tAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGt                             >  1:350627/1‑150 (MQ=255)
                                                                                                                         aaaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTTCGGCCTTGGCCAGTc                            >  2:128582/1‑150 (MQ=255)
                                                                                                                         aaaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTc                            <  2:430351/150‑1 (MQ=255)
                                                                                                                            ccGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                         <  2:1115588/150‑1 (MQ=255)
                                                                                                                             agtaCAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGTTTTACCTTGCTCATCCTTTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAgg                        <  1:337687/147‑1 (MQ=255)
                                                                                                                              ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGt                       >  1:1361291/1‑150 (MQ=255)
                                                                                                                                 cAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGtttt                    >  2:966274/1‑150 (MQ=255)
                                                                                                                                  aaTAAAAACTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGTAAATCTTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACTGCTCGAGTACGGCCTTGGCCAGTCCAGGTTTTa                   <  1:163585/150‑1 (MQ=255)
                                                                                                                                   aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAg                  >  2:620931/1‑150 (MQ=255)
                                                                                                                                    tAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGAGCAATCTACAAAAGATTATTGCATGATCAGGAATTCGTCGAGTGATTTACCTTGCACATCCATTGCTTTTTTTAATACCGCTGGAGTACGGCCTTGGCAAGTACAGTTTTTatt                 >  2:1083810/1‑148 (MQ=255)
                                                                                                                                    cttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGt                 >  1:959856/4‑150 (MQ=255)
                                                                                                                                     ttaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGtt                >  1:704368/3‑150 (MQ=255)
                                                                                                                                     aaaaaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGtt                >  2:475076/1‑150 (MQ=255)
                                                                                                                                       aaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc              >  2:1713006/1‑150 (MQ=255)
                                                                                                                                        ataTCCGGCCGGTGGCGGGATTTTAAGCATGTGGAAGCTACAAAAGATTATTTCTTGATCAGGTATTCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg             <  2:262003/148‑1 (MQ=255)
                                                                                                                                        aaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg             >  2:582917/1‑150 (MQ=255)
                                                                                                                                         aaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGc            >  2:614364/1‑150 (MQ=255)
                                                                                                                                            cccGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGt         >  2:1013569/1‑150 (MQ=255)
                                                                                                                                             ccgccgCTGGCGGGATTTTTAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATAGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGGTTTAGTTTCGCCGtt        >  1:114529/1‑150 (MQ=255)
                                                                                                                                                   gtGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTc  <  2:518325/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTCGTC  >  NC_000913/1292326‑1292622

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