Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,292,402 C→T 5.1% intergenic (+37/+107) 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,4020CT 5.1% 312.7 / 2.3 99intergenic (+37/+107)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base C (48/43);  new base T (2/3);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

GCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTC  >  NC_000913/1292253‑1292550
                                                                                                                                                     |                                                                                                                                                    
gCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGc                                                                                                                                                      <  2:1198434/150‑1 (MQ=255)
          aTTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGcc                                                                                                                                            <  2:1177994/150‑1 (MQ=255)
               gTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGtttt                                                                                                                                       >  1:1042538/1‑150 (MQ=255)
                ttACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTCGATTGCTCAACGCCGtttcg                                                                                                                                      <  2:1637275/150‑3 (MQ=255)
                ttACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTCGATTGCTCAACGCCGtttcg                                                                                                                                      <  1:1413733/150‑3 (MQ=255)
                  aCATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAATAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGttttttt                                                                                                                                    <  2:723904/150‑1 (MQ=255)
                    aTGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAATCCTGGCGTGAAGAAGAGATGGGAATTAAGAAGTAACATCCGTATCGGTGTTATCAACGAAACGGCGTTGAGCAATCGACGCCGTATTTTtat                                                                                                                                  >  1:1104719/1‑150 (MQ=255)
                     tGCAGGCCTTCGTTGAATACGGTATTCTTCATAACAACATTGGCACGGAATTTAAAGACTTGCTTGAAGAAGAGATGGGCATTAAGAAGTAACAAACGTATCGGTTTTATACACTAAACTTAGTTTAGAAATCGACGCCGTTATTTtata                                                                                                                                 >  1:1015577/1‑150 (MQ=255)
                     tGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTCGATTGCTCAACGCCGtttcgtggat                                                                                                                                 <  1:149858/150‑8 (MQ=255)
                     tGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGCCGTTGAGCAATCGACGCCGTTTTTTtata                                                                                                                                 >  1:841333/1‑150 (MQ=255)
                       cAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGc                                                                                                                               >  2:392798/1‑150 (MQ=255)
                        aGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGct                                                                                                                              <  1:696086/150‑1 (MQ=255)
                         ggCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGctt                                                                                                                             >  2:1645677/1‑150 (MQ=255)
                         ggCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGctt                                                                                                                             <  2:678643/150‑1 (MQ=255)
                           ccTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGcttat                                                                                                                           >  1:572027/1‑150 (MQ=255)
                           ccTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCAGTTTTTTTATAGcttat                                                                                                                           >  1:1299472/1‑150 (MQ=255)
                            cTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGcttatt                                                                                                                          >  2:599911/1‑150 (MQ=255)
                                 ttGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGcttattcttat                                                                                                                     >  2:311165/1‑150 (MQ=255)
                                 ttGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGcttattcttat                                                                                                                     <  2:547726/150‑1 (MQ=255)
                                   gAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCttatta                                                                                                                   <  1:1106627/150‑1 (MQ=255)
                                   gAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCttatta                                                                                                                   >  1:270129/1‑150 (MQ=255)
                                     aTACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTaaa                                                                                                                 >  1:810931/1‑150 (MQ=255)
                                      tACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGAATCGGTGTTATCCACGAAACGTCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAAt                                                                                                                <  1:95984/150‑1 (MQ=255)
                                           tATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCt                                                                                                           <  2:1487842/150‑1 (MQ=255)
                                               cGTCATTACACCCTTGGCACGGAATTTAAAGCCAGTCTTGAAGAAGAGATGGCCAATAAGAAGTAACATCCGAATCGGTGTAATCCACGAAACCCCGTTGAGCAAACGACACCGTATTTTTATAGCTTATTCTTATTAAAATGTCTTaca                                                                                                       >  1:1098928/1‑148 (MQ=255)
                                               cGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGAATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGACTTaaa                                                                                                       >  2:302252/1‑150 (MQ=255)
                                                  cATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCg                                                                                                    >  2:559154/1‑150 (MQ=255)
                                                  cATAACACCCTTGGCACGGAATTAAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCg                                                                                                    <  2:1375708/150‑1 (MQ=255)
                                                   aTAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCgg                                                                                                   <  2:1060258/150‑1 (MQ=255)
                                                    tAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGa                                                                                                  <  2:944937/150‑1 (MQ=255)
                                                     aaCACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGAc                                                                                                 >  2:441230/1‑150 (MQ=255)
                                                     aaCACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGAc                                                                                                 >  2:704557/1‑150 (MQ=255)
                                                       cacCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACaa                                                                                               <  2:572027/150‑1 (MQ=255)
                                                       cacCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACaa                                                                                               >  2:1248393/1‑150 (MQ=255)
                                                        acCCTTGGCTCGTAATTTAAAGCCTGGCTTGAAGAAGAGAAGGGCATTAAGAAGTAACATCCGTATCGGTTTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAAt                                                                                              <  2:347084/150‑1 (MQ=255)
                                                        acCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAAt                                                                                              >  1:1136960/1‑150 (MQ=255)
                                                        acCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATa                                                                                             >  2:432455/1‑150 (MQ=255)
                                                          ccTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATaa                                                                                            <  1:666327/150‑1 (MQ=255)
                                                             tGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATaaaaa                                                                                         <  2:1168066/150‑1 (MQ=255)
                                                                cACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTAGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATc                                                                                      <  1:269201/150‑1 (MQ=255)
                                                                  cGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATccc                                                                                    <  1:1697011/150‑1 (MQ=255)
                                                                   gtaaTTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACCTCCGAATCGGTTTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCccg                                                                                   <  1:819466/148‑1 (MQ=255)
                                                                   ggAAGTTCACGCCTGGCTTGAAGACTCGATGCGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAATCCGGACAATAAAAAATCccg                                                                                   <  2:153933/150‑1 (MQ=255)
                                                                          aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTgg                                                                            >  1:241747/1‑150 (MQ=255)
                                                                          aaaGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTgg                                                                            >  1:1421434/1‑150 (MQ=255)
                                                                            aGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCg                                                                          <  1:392798/150‑1 (MQ=255)
                                                                             gcctggcctgAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCgg                                                                         >  2:1186416/1‑150 (MQ=255)
                                                                              ccTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACAGGACAATAAAAAATCCCGCCGCTGGCggg                                                                        >  1:738110/1‑150 (MQ=255)
                                                                               cTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGa                                                                       >  2:896731/1‑150 (MQ=255)
                                                                                 ggCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTATTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtt                                                                     >  1:1039682/1‑150 (MQ=255)
                                                                                   cTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtttt                                                                   >  2:19331/1‑150 (MQ=255)
                                                                                   cTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtttt                                                                   <  1:572913/150‑1 (MQ=255)
                                                                                    ctGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTa                                                                  <  1:1186416/149‑1 (MQ=255)
                                                                                     tGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTaa                                                                 <  1:28409/150‑1 (MQ=255)
                                                                                      gaagaaGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCAGCGGCGGGAtttttta                                                                <  2:1136960/150‑4 (MQ=255)
                                                                                       aagaagCGATCGGCATTAAGAAGTAACATCCGTATCGGTGTTAACCACGAAACCGCGTTGAGCAATCGACGCCGTTTTTTTAAAGCTTATACTTATTAAATTGTCTTAAACCGGACCATAAAAAATCCCGCCGCTGGCGGGATTTTAAGc                                                               <  1:1171628/150‑1 (MQ=255)
                                                                                        agaagaTATGGGCATTAAGTAGTAACATCCGTATCGGTGTTATCCAAGAAACGGCGTTTAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                              <  1:336349/150‑1 (MQ=255)
                                                                                        agaagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCa                                                              >  1:625354/1‑150 (MQ=255)
                                                                                          aagaGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTATTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGTGATTTTAAGCaat                                                            >  1:450812/1‑149 (MQ=255)
                                                                                           agagATGGGCATTAAGAAGTAACATCCGTATTGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGt                                                           >  2:1336837/1‑150 (MQ=255)
                                                                                             agaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                <  1:1422465/143‑1 (MQ=255)
                                                                                             agaTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                >  2:1422465/1‑143 (MQ=255)
                                                                                              caTGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCa                                                        <  1:340917/149‑1 (MQ=255)
                                                                                                tGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAAt                                                      <  1:19331/150‑1 (MQ=255)
                                                                                                 tggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATc                                                     <  1:1459451/149‑1 (MQ=255)
                                                                                                  ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTGTAAGCAAGTGCAATCt                                                    >  2:34385/1‑150 (MQ=255)
                                                                                                        aagaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaaa                                              <  1:355474/150‑1 (MQ=255)
                                                                                                          gaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                            >  2:1677185/1‑150 (MQ=255)
                                                                                                           aagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGat                                           <  2:1680021/150‑1 (MQ=255)
                                                                                                                  aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                    <  1:1645677/150‑1 (MQ=255)
                                                                                                                  aTCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                    <  2:270129/150‑1 (MQ=255)
                                                                                                                    ccGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                  <  2:1421434/150‑1 (MQ=255)
                                                                                                                    acGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTg                                  <  2:1219160/149‑1 (MQ=255)
                                                                                                                       tATCTGTTTTATCCACGAAACTGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATc                               <  1:311165/150‑1 (MQ=255)
                                                                                                                          cGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                            <  2:625354/150‑1 (MQ=255)
                                                                                                                           ggTTTTATCCAAGAAACGGCGTTGAGCTATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGa                           <  2:235243/150‑1 (MQ=255)
                                                                                                                            gtgtTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaa                          >  1:465939/1‑150 (MQ=255)
                                                                                                                               ttATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtc                       >  2:134611/1‑150 (MQ=255)
                                                                                                                                tATCCACGAAACGGCGTTGAGCAATCGACGCCCTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcg                      >  2:1282414/1‑150 (MQ=255)
                                                                                                                                 ttCCAAGAAACGGCGTTGAGAAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGTCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                     <  2:841333/149‑1 (MQ=255)
                                                                                                                                 tataaaaaaaCGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                     >  2:953728/8‑150 (MQ=255)
                                                                                                                                   ccACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                   >  2:1323149/1‑150 (MQ=255)
                                                                                                                                   ccACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                   >  2:1701374/1‑150 (MQ=255)
                                                                                                                                    cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                  >  1:86361/1‑150 (MQ=255)
                                                                                                                                    cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                  >  1:413697/1‑150 (MQ=255)
                                                                                                                                    aaagAAAAGGAGTTGAGCAATAGAAGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                  <  2:1015577/147‑1 (MQ=255)
                                                                                                                                     aCGAAACGGCGTTGATCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                 <  2:1253049/150‑1 (MQ=255)
                                                                                                                                     aCGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAg                 <  2:1634177/150‑1 (MQ=255)
                                                                                                                                     aCGAAACGGAGTTTAGCAATCGACGCCCAATTATAATAGCTAATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCACTGGCGCGATTTAAAGCAAGTGCAAACTACAAAACATTATTGCTTGAACAGGAAATcgacgaa                 >  1:1374542/1‑149 (MQ=255)
                                                                                                                                      cGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                >  1:842333/1‑150 (MQ=255)
                                                                                                                                      cGAAACGGAGTTGAGAAATCGACTCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                <  2:450812/150‑1 (MQ=255)
                                                                                                                                      aaaaaCGGCGTCGATTGCTCAACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAgg                >  1:1274359/3‑150 (MQ=255)
                                                                                                                                         aaCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAt             >  2:1300947/1‑150 (MQ=255)
                                                                                                                                          aCGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAtt            >  1:781834/1‑150 (MQ=255)
                                                                                                                                           cGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                    >  1:1216827/1‑125 (MQ=255)
                                                                                                                                           cGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCt                                    <  2:1216827/125‑1 (MQ=255)
                                                                                                                                             gCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAc         >  2:1556905/1‑150 (MQ=255)
                                                                                                                                               gTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGAAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGTATTTACCt       >  1:1056392/1‑150 (MQ=255)
                                                                                                                                                   aGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCt   >  2:1090670/1‑150 (MQ=255)
                                                                                                                                                    gCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc  >  2:696832/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGGTAATAAATTAGGTTACATGCAGGCCTTCGTTGAATACGGTATTCGTCATAACACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTC  >  NC_000913/1292253‑1292550

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