Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,212,555 C→A 7.8% intergenic (‑552/+148) iraM ← / ← ycgX RpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,212,5550CA 7.8% 381.5 / 8.5 114intergenic (‑552/+148)iraM/ycgXRpoS stabilzer during Mg starvation, anti‑RssB factor/DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base C (28/74);  new base A (10/2);  total (38/76)
Fisher's exact test for biased strand distribution p-value = 2.61e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAG  >  NC_000913/1212406‑1212702
                                                                                                                                                     |                                                                                                                                                   
gCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATAAAATTTAAAATATc                                                                                                                                                     >  1:1947502/1‑150 (MQ=255)
gCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGATGTTTTTATAAAAAAGTTAACCGTTATAAAATAACGTATAAAAAGTCAAGTGATCACATTTCCAATATc                                                                                                                                                     >  1:954358/1‑150 (MQ=255)
     aaTACATATAACCTGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGtt                                                                                                                                                <  2:2142840/150‑1 (MQ=255)
       tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGa                                                                                                                                              <  2:1992513/150‑1 (MQ=255)
       tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGa                                                                                                                                              <  2:2013218/150‑1 (MQ=255)
       tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGa                                                                                                                                              <  1:1703083/150‑1 (MQ=255)
       tACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGa                                                                                                                                              <  2:61964/150‑1 (MQ=255)
        aCATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGAGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGAt                                                                                                                                             >  2:1219686/1‑150 (MQ=255)
         cATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATa                                                                                                                                            <  2:1331856/150‑1 (MQ=255)
         cATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATa                                                                                                                                            <  1:1430113/150‑1 (MQ=255)
           tataACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGGTTTTATAAATATTTTAAACGTTATAAAAAAAAGTATCAAAAGTCAAGTGATGACAATTCAAATATAAAATTGTTAGt                                                                                                                                          >  1:815866/1‑150 (MQ=255)
            ataACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTa                                                                                                                                         >  2:2073174/1‑150 (MQ=255)
              aaCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTAtt                                                                                                                                       >  1:1230937/1‑150 (MQ=255)
              aaCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATAAAAAGTAAAGTTATAACATTTAAAATATCAATTTTATAGTAtt                                                                                                                                       >  1:1480359/1‑150 (MQ=255)
               aCCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTa                                                                                                                                      <  1:2829360/150‑1 (MQ=255)
                ccGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTTTAAATATGTTAAACGTTATAAAATAAAGTATCAAAAGTCAAGTGATCACATTTAAAATATAAAGTTGAGAGTATTAg                                                                                                                                     >  1:1052134/1‑150 (MQ=255)
                 cGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGt                                                                                                                                    >  2:2626748/1‑150 (MQ=255)
                 cGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATAAAAAGTCAAGTGATAACATTTAAAATATAAATTTGTTAGTATTAGt                                                                                                                                    >  1:492637/1‑150 (MQ=255)
                     tgatgaCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAAACGTTATAAAATAACGTATCAAAAGTCAAGTTTTAACAATTAAAATATAAATTTGATAGTTTTAGtgggg                                                                                                                                >  1:1720882/1‑146 (MQ=255)
                          aCTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGatt                                                                                                                           <  1:1242455/150‑1 (MQ=255)
                              aaCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttattt                                                                                                                       <  1:63257/150‑1 (MQ=255)
                               aCCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGACAGTATTAGTCTGGTGAttattta                                                                                                                      <  2:425250/150‑1 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTTAACGTTATAAAATAAAGTATAAAAAGTAAAGTGATCACATTTAAAATATCAAGTTGGTATTATTAGTCTGGTGTttatttat                                                                                                                     >  1:2895569/1‑150 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                                                                                                     <  2:1640678/150‑1 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                                                                                                     <  2:2063454/150‑1 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGAttatttat                                                                                                                     <  2:669934/150‑1 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATAATTTAAACCGTAAAAAAAAAAGTATAAAAAATAAAGTGATAACAATTCAAAAATAAAGTTTATATTATTAGTATGGTGAttatttat                                                                                                                     >  1:160850/1‑150 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGTTTTTATAAATATGTTTAACGTTATAAAATAACGTATCAAAAGTCAAGTGATAAAAATTAAAAAATAAAGTTGATAGTATTAGTATGGTGAttatttat                                                                                                                     >  1:160896/1‑150 (MQ=255)
                                ccAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTGGTTTTTATAAAAATGTTAAACGTTATAAAATAAAGTATAAAAAGTCAAGTGATCACAATTAAAATATAAAGTTGATAGTATTAGTATGGTGAttatttat                                                                                                                     >  1:161054/1‑150 (MQ=255)
                                 cAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATg                                                                                                                    <  2:1342547/150‑1 (MQ=255)
                                 cAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTAAAAATAATTTAAACCTTAAAAAATAAAGTATCAAAATTCAAGTGATAACAATTAAAATAACAAGTTGATAGTATTAGTCTGGTGATTATTTTTg                                                                                                                    >  2:252179/1‑150 (MQ=255)
                                  aaaaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATgg                                                                                                                   <  2:2456574/150‑1 (MQ=255)
                                    aaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGt                                                                                                                 <  1:68121/150‑1 (MQ=255)
                                    aaTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACAGTTCAAATCTCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGt                                                                                                                 <  1:2289472/150‑1 (MQ=255)
                                     aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATAAGTTAAACGTTATAAAATAAAGTATAAAAAATAAAATGATAACAATTAAAATATAAAGTTTATAGTTTTATTCTTGTGATTATTTAtgtgtg                                                                                                                >  1:1790564/1‑150 (MQ=255)
                                     aTACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAAAATTTTAAAACTAATAAAAAAAAGTATAAAAATAAAAGTGATCACAATTCAAAAAAAAAATTTATATTATTAGTGTGGTTTTTTTTTATTGGTg                                                                                                                >  1:1796119/1‑150 (MQ=255)
                                        cATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACa                                                                                                             <  2:2390794/150‑1 (MQ=255)
                                          tgtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAAt                                                                                                           <  1:1948337/150‑1 (MQ=255)
                                           gtgCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATTAATATTTTAAACATTATAAAAAAAAGTATAAAAAGTAAACAGATAAAAATTAAAAAATAAAATTTATAGTATTAGTGTGGTGGTTATTTATTGGTGaaaaaa                                                                                                          >  1:2591515/1‑148 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        <  2:156958/150‑1 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        <  2:5380/150‑1 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        <  1:200284/150‑1 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        <  1:816965/150‑1 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        >  2:660800/1‑150 (MQ=255)
                                             gCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATaaa                                                                                                        <  1:200291/150‑1 (MQ=255)
                                                 tAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGa                                                                                                    <  1:2680000/150‑1 (MQ=255)
                                                 tAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGa                                                                                                    <  2:1343037/150‑1 (MQ=255)
                                                 tAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGa                                                                                                    <  1:803659/150‑1 (MQ=255)
                                                  aaGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAATTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGAc                                                                                                   >  2:1537033/1‑150 (MQ=255)
                                                  aaGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGAc                                                                                                   >  2:402939/1‑150 (MQ=255)
                                                   aGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACa                                                                                                  <  2:714248/150‑1 (MQ=255)
                                                       aaacaaCGGCCGTGCCACACGTCCGTGGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCCAAAGTCCCGGGATCCCATTTCAAATATCAAGTTGATCGTATTAGTCTGGGGATGATTTATGGGTGACAATAAAAAGCCAGGAt                                                                                              >  1:1201676/5‑150 (MQ=255)
                                                               gaggtgtGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAatcatc                                                                                      <  1:254038/147‑1 (MQ=255)
                                                                 cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCATGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa                                                                                    <  1:636251/150‑1 (MQ=255)
                                                                 cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa                                                                                    <  2:1858753/150‑1 (MQ=255)
                                                                 cGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCa                                                                                    <  2:213480/150‑1 (MQ=255)
                                                                   tgtgGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATa                                                                                  <  2:258592/150‑1 (MQ=255)
                                                                    gccacACGTCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATag                                                                                 >  2:1579156/5‑150 (MQ=255)
                                                                      ggCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATagag                                                                               <  1:110807/150‑1 (MQ=255)
                                                                        cACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCAAATTTCAAATATCCATTTGATAGTATTAGTCTGGGGATGATTTATGGGTGACAATAAAAAGACAGTATTATTAATCAATAGAGAt                                                                             >  2:7334/1‑150 (MQ=255)
                                                                         aCGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATa                                                                            <  1:2073174/150‑1 (MQ=255)
                                                                               gttgttTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTctct                                                                      <  1:113758/150‑1 (MQ=255)
                                                                               gttgttTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTctct                                                                      <  2:1666/150‑1 (MQ=255)
                                                                               gttgttTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTctct                                                                      >  2:2394097/1‑150 (MQ=255)
                                                                                ttgttTTTATTAATATGTTAACCGTTATAAAATAACGTATCATAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTTGTGTTTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTg                                                                     <  2:1482267/150‑1 (MQ=255)
                                                                                    ttttATAAATATGTTAACCGTTATAAAATAACGTATCTAAAGTCAAGTTATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACt                                                                 <  2:495498/150‑1 (MQ=255)
                                                                                    ttttATAAATATGTTAACCGTTATAAAATAACGTATCAATAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACt                                                                 <  2:2119008/150‑1 (MQ=255)
                                                                                     tttATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACtt                                                                >  1:1064650/1‑150 (MQ=255)
                                                                                      ttATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACttt                                                               <  1:2090300/150‑1 (MQ=255)
                                                                                         taAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTAt                                                            <  1:1338913/150‑1 (MQ=255)
                                                                                            atatGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTc                                                         <  1:1670070/150‑1 (MQ=255)
                                                                                            atatGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATAAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTc                                                         <  1:178916/150‑1 (MQ=255)
                                                                                                   aaCCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCTAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATg                                                  <  2:1190/150‑1 (MQ=255)
                                                                                                    aCCGTTCTAAAATAACGTCTCAAAAGTCAAGTGATCACATGTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGc                                                 <  1:1784557/150‑1 (MQ=255)
                                                                                                      cGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTa                                               >  2:1056832/1‑150 (MQ=255)
                                                                                                       gTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTaa                                              <  2:264344/150‑1 (MQ=255)
                                                                                                         tataAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATg                                            <  2:2883303/150‑1 (MQ=255)
                                                                                                          ataAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGc                                           <  1:660927/150‑1 (MQ=255)
                                                                                                              aaTAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTa                                       <  2:2278279/150‑1 (MQ=255)
                                                                                                                    gTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAAt                                 <  2:492637/150‑1 (MQ=255)
                                                                                                                    gTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAAt                                 <  1:2626748/150‑1 (MQ=255)
                                                                                                                      aTCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTa                               >  1:2925473/1‑150 (MQ=255)
                                                                                                                        cAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATg                             <  2:1780679/150‑1 (MQ=255)
                                                                                                                         aaaaGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGa                            <  2:917941/150‑1 (MQ=255)
                                                                                                                         aaaaGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGa                            >  1:1115899/1‑150 (MQ=255)
                                                                                                                         aaaaGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGa                            >  1:1115923/1‑150 (MQ=255)
                                                                                                                          aaaGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGaa                           <  1:512967/150‑1 (MQ=255)
                                                                                                                                   tGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCtt                  >  2:2285893/1‑150 (MQ=255)
                                                                                                                                     aTCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTaa                >  1:894733/1‑150 (MQ=255)
                                                                                                                                      tCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAg               <  1:401539/150‑1 (MQ=255)
                                                                                                                                          aTTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATc           <  1:812588/150‑1 (MQ=255)
                                                                                                                                          aTTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATc           <  1:1537033/150‑1 (MQ=255)
                                                                                                                                             tCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCaa        >  1:2450661/1‑150 (MQ=255)
                                                                                                                                                aaTATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACtt     <  1:1199417/150‑1 (MQ=255)
                                                                                                                                                  tatCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGACTTATGAAGCATTTCTTAAGTATCCAACTTTa   <  1:208698/150‑1 (MQ=255)
                                                                                                                                                   atCAATTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAg  <  2:1146563/150‑1 (MQ=255)
                                                                                                                                                   atCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAg  <  1:357850/150‑1 (MQ=255)
                                                                                                                                                   atCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAg  <  2:301163/150‑1 (MQ=255)
                                                                                                                                                   atCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAg  <  1:1813526/150‑1 (MQ=255)
                                                                                                                                                   atCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAg  <  1:328120/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCTAAAATACATATAACCGGATGATGACTAAACCAAAATACATGTGCGTTAAGTATTGAAACGGACGTGTGGCACGGCCGTTGTTTTTATAAATATGTTAACCGTTATAAAATAACGTATCAAAAGTCAAGTGATCACATTTCAAATATCAAGTTGATAGTATTAGTCTGGTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAG  >  NC_000913/1212406‑1212702

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