Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,473 C→G 5.2% intergenic (+108/+36) 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,4730CG 5.2% 540.6 / 14.9 154intergenic (+108/+36)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base C (62/84);  new base G (5/3);  total (67/87)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

Only 100 of 156 total aligned reads displayed.

AAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTC  >  NC_000913/1292327‑1292619
                                                                                                                                                  |                                                                                                                                                  
aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                  >  4:1435895/1‑149 (MQ=255)
      tGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGa                                                                                                                                            <  2:687057/149‑1 (MQ=255)
         cTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAttt                                                                                                                                         <  4:277860/149‑1 (MQ=255)
          ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCATGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttt                                                                                                                                        <  3:1430215/149‑1 (MQ=255)
          ttGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtttt                                                                                                                                        >  1:1470285/1‑149 (MQ=255)
            gaagaaGAGATGGGCATTAAGAAGTAATATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                      <  1:1272358/149‑1 (MQ=255)
            gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                                                                                      <  1:1397696/149‑1 (MQ=255)
              agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                                                                                                    >  4:340859/1‑149 (MQ=255)
              agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                                                                                                    >  1:190850/1‑149 (MQ=255)
               gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                                                                                                   >  2:518121/1‑149 (MQ=255)
               gaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                                                                                                   >  1:1574709/1‑149 (MQ=255)
                aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCaa                                                                                                                                  >  3:2447/1‑149 (MQ=255)
                aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATACCGCCGCTGGCGGGATTTTAAGCaa                                                                                                                                  >  3:1359183/1‑149 (MQ=255)
                 agagATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAg                                                                                                                                 <  3:1161003/149‑1 (MQ=255)
                     aTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCa                                                                                                                             >  4:1530772/1‑149 (MQ=255)
                          cATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTa                                                                                                                        <  1:518121/149‑1 (MQ=255)
                           aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                       <  4:16233/149‑1 (MQ=255)
                           aTTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                       >  4:549663/1‑149 (MQ=255)
                             tAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaa                                                                                                                     <  2:1470285/149‑1 (MQ=255)
                               agaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                                                                                                   <  1:77486/149‑1 (MQ=255)
                                 aagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                                 >  4:955627/1‑149 (MQ=255)
                                  agTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGat                                                                                                                <  1:94075/149‑1 (MQ=255)
                                     accttCCTTATCTGTGTTATCCACCAAAGTAGGTTTAATAATGTACCCCCTTTTTTTATAGCTTATTTTTATTAAATTGTCTCAATCCGGATAATAAAAAATCCCGCCGCTGGCGGGATTTAAGCAAGTGCAATCTACAAAAGattatt                                                                                                            <  4:1429246/145‑1 (MQ=255)
                                     aaCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat                                                                                                             >  4:849065/1‑149 (MQ=255)
                                     aaCATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGattat                                                                                                             <  1:728059/149‑1 (MQ=255)
                                        aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                                          <  2:968794/149‑1 (MQ=255)
                                        aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                                          >  1:841874/1‑149 (MQ=255)
                                        aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGc                                                                                                          >  2:1142808/1‑149 (MQ=255)
                                         tCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                                                                                         <  2:319241/149‑1 (MQ=255)
                                          ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                                        <  3:212684/149‑1 (MQ=255)
                                            gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                       >  3:833895/1‑132 (MQ=255)
                                            gTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                                       <  4:833895/132‑1 (MQ=255)
                                             tATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                                     <  4:1636269/149‑1 (MQ=255)
                                               tCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCa                                                                                                   >  1:945835/1‑149 (MQ=255)
                                                          cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                                        >  2:101286/1‑149 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTTTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       <  4:135839/149‑1 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       <  3:849065/149‑1 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       >  2:920892/1‑149 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       >  1:880174/1‑149 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       >  2:1654705/1‑149 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       >  3:1556595/1‑149 (MQ=255)
                                                           aCGAAACGGCGTTGAGCAATCGACGCAGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAATCCGGACTATAAAAAATCCCGCCGCTGGAGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                                       <  2:852329/149‑1 (MQ=255)
                                                            cGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                                                                                      <  1:101286/149‑1 (MQ=255)
                                                             gAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                                                                                     >  3:1352018/1‑149 (MQ=255)
                                                              aaaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAggg                                                                                    <  2:1290177/149‑1 (MQ=255)
                                                               aaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                   <  4:285500/149‑1 (MQ=255)
                                                                 cGGCGTTGAGCAATTGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAtt                                                                                 <  2:380930/149‑1 (MQ=255)
                                                                  ggCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGa                                                                                   >  3:1390489/1‑146 (MQ=255)
                                                                    cGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc                                                                              <  2:1683329/149‑1 (MQ=255)
                                                                    cGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc                                                                              <  1:712689/149‑1 (MQ=255)
                                                                         aGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                                         <  4:603205/149‑1 (MQ=255)
                                                                              tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                                        >  4:1106437/1‑113 (MQ=255)
                                                                              tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                    >  3:1711919/1‑149 (MQ=255)
                                                                              tCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATc                                                                    <  3:646253/149‑1 (MQ=255)
                                                                               cGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATcc                                                                   <  2:913491/149‑1 (MQ=255)
                                                                               cGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATcc                                                                   >  1:1054510/1‑149 (MQ=255)
                                                                               cGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATcc                                                                   <  1:1752430/149‑1 (MQ=255)
                                                                                 aCGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCAt                                                                 <  2:1574709/149‑1 (MQ=255)
                                                                                   gCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                                               <  4:400354/149‑1 (MQ=255)
                                                                                          ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                                        <  2:119101/149‑1 (MQ=255)
                                                                                          ttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                                        <  3:219924/149‑1 (MQ=255)
                                                                                            ttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGa                                                      <  1:1286869/149‑1 (MQ=255)
                                                                                               taGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTa                                                   >  4:108326/1‑149 (MQ=255)
                                                                                                 gCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACa                                                 <  1:1142808/149‑1 (MQ=255)
                                                                                                   ttattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGc                                               >  3:71844/1‑149 (MQ=255)
                                                                                                    tattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCTGCCGCTGGCTGTATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc                                                                              >  4:913265/1‑117 (MQ=255)
                                                                                                    tattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc                                                                              <  3:913265/117‑1 (MQ=255)
                                                                                                    tattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                                              <  1:159837/149‑1 (MQ=255)
                                                                                                    tattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCt                                              <  4:1359183/149‑1 (MQ=255)
                                                                                                     attcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTg                                             <  3:1504185/149‑1 (MQ=255)
                                                                                                         ttattaAGTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCCATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTa                                                                               <  1:261584/111‑1 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                  <  3:898563/44‑1 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTg                                                                                                                                                  >  4:898563/1‑44 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTa                                                                               >  2:261584/1‑111 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                                          <  2:1224984/132‑1 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                                          >  1:1224984/1‑132 (MQ=255)
                                                                                                         ttattaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGt                                         >  2:1538900/1‑149 (MQ=255)
                                                                                                            ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                      >  4:46018/1‑149 (MQ=255)
                                                                                                            ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGACTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACg                                      >  3:698617/1‑149 (MQ=255)
                                                                                                                  tGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTg                                <  3:73907/149‑1 (MQ=255)
                                                                                                                   gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTgg                               <  1:39024/149‑1 (MQ=255)
                                                                                                                    tCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                              >  2:1152531/1‑149 (MQ=255)
                                                                                                                    tCTTAAACCGAACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGc                              <  2:99982/149‑1 (MQ=255)
                                                                                                                     cTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGcc                             <  1:1654705/149‑1 (MQ=255)
                                                                                                                      ttAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCa                            >  2:1325002/1‑149 (MQ=255)
                                                                                                                        aaaCCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                                     <  2:132793/138‑1 (MQ=255)
                                                                                                                            cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGCAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGACAGTCCAg                      <  4:1567528/149‑1 (MQ=255)
                                                                                                                            cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                      >  1:1489340/1‑149 (MQ=255)
                                                                                                                            cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATACATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCac                      <  2:1489340/149‑2 (MQ=255)
                                                                                                                            cGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTCGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATAGATTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg                      >  4:1498091/1‑149 (MQ=255)
                                                                                                                             ggACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAgg                     >  3:395413/1‑149 (MQ=255)
                                                                                                                                  aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                <  3:955627/149‑1 (MQ=255)
                                                                                                                                  aTAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTa                <  3:46018/149‑1 (MQ=255)
                                                                                                                                   tAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAg               >  2:742728/1‑149 (MQ=255)
                                                                                                                                    ctaaaaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTAc                                       <  1:1585151/122‑1 (MQ=255)
                                                                                                                                       aaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc           <  2:1014536/149‑1 (MQ=255)
                                                                                                                                       aaaTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTc           <  3:549663/149‑1 (MQ=255)
                                                                                                                                        aaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg          >  2:1618344/1‑149 (MQ=255)
                                                                                                                                        aaTCCCGCCAGCGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCg          >  1:1121174/1‑149 (MQ=255)
                                                                                                                                             cgccgcTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGtt     >  1:1117542/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTTTTAGTTTCGCCGTTTTC  >  NC_000913/1292327‑1292619

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