Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,411,384 C→T 10.1% A189T (GCC→ACC)  ttcA ← tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,411,3840CT10.1% 2182.1 / 166.6 705A189T (GCC→ACC) ttcAtRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base C (337/297);  new base T (25/46);  total (362/343)
Fisher's exact test for biased strand distribution p-value = 5.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 707 total aligned reads displayed.

GCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTT  >  NC_000913/1411237‑1411533
                                                                                                                                                   |                                                                                                                                                     
   cTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCg                                                                                                                                                  >  1:679953/1‑150 (MQ=255)
     ggATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGa                                                                                                                                                <  1:10018954/150‑1 (MQ=255)
      gATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGAc                                                                                                                                               >  1:1367839/1‑150 (MQ=255)
          aCGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAAt                                                                                                                                           <  1:7697301/150‑1 (MQ=255)
            gTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATaa                                                                                                                                         <  2:6219218/150‑1 (MQ=255)
                aTCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGat                                                                                                                                     <  1:11206393/150‑1 (MQ=255)
                 tCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGata                                                                                                                                    >  2:11239734/1‑150 (MQ=255)
                 tCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGata                                                                                                                                    >  1:11717398/1‑150 (MQ=255)
                      gTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGttt                                                                                                                               >  1:7017329/1‑150 (MQ=255)
                      gTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGttt                                                                                                                               <  1:8891953/150‑1 (MQ=255)
                       tCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTg                                                                                                                              <  2:6511736/150‑1 (MQ=255)
                           gCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcat                                                                                                                          >  1:10874757/1‑150 (MQ=255)
                           gCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcat                                                                                                                          <  2:9521458/150‑1 (MQ=255)
                            cAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcatc                                                                                                                         >  1:366046/1‑150 (MQ=255)
                            cAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcatc                                                                                                                         >  2:4319536/1‑150 (MQ=255)
                            cAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTATTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGATTGCcatc                                                                                                                         >  2:2998857/1‑150 (MQ=255)
                             aaCATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCcatca                                                                                                                        <  1:6013193/150‑1 (MQ=255)
                             aaCATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcatca                                                                                                                        <  1:1837111/150‑1 (MQ=255)
                             aaCATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCcatca                                                                                                                        <  1:6552816/150‑1 (MQ=255)
                                aTGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCg                                                                                                                     <  1:4244044/150‑1 (MQ=255)
                                  gTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCt                                                                                                                   <  1:11824383/150‑1 (MQ=255)
                                              cTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATCATCGGGAACGCTTTTGCAGCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCgg                                                                                                       >  2:10049054/1‑150 (MQ=255)
                                              cTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCgg                                                                                                       >  2:1578885/1‑150 (MQ=255)
                                                    ttGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCAt                                                                                                 >  2:9767248/1‑150 (MQ=255)
                                                    ttGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCAt                                                                                                 <  1:972971/150‑1 (MQ=255)
                                                    ttGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCAt                                                                                                 >  1:2779704/1‑150 (MQ=255)
                                                      gCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATAc                                                                                               >  2:11493173/1‑150 (MQ=255)
                                                       cAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATAcc                                                                                              >  2:11089483/1‑150 (MQ=255)
                                                          gTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCttt                                                                                           <  2:8125055/150‑1 (MQ=255)
                                                           ttAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTc                                                                                          <  1:4404883/150‑1 (MQ=255)
                                                           ttAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTc                                                                                          >  2:6647525/1‑150 (MQ=255)
                                                           ttAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTc                                                                                          >  2:2690468/1‑150 (MQ=255)
                                                              ggCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATc                                                                                       <  1:6202471/150‑1 (MQ=255)
                                                               gCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCt                                                                                      >  1:3025435/1‑150 (MQ=255)
                                                                 tgtgAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTa                                                                                    <  1:5110141/150‑1 (MQ=255)
                                                                  gtgAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTc                                                                                          >  1:4949501/1‑143 (MQ=255)
                                                                   tgAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTAcc                                                                                  >  2:873488/1‑150 (MQ=255)
                                                                      aCCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTAccgcc                                                                               <  1:9065088/150‑1 (MQ=255)
                                                                      aCCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTAccgcc                                                                               <  2:4097841/150‑1 (MQ=255)
                                                                          cacaGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAg                                                                           >  1:3346112/1‑150 (MQ=255)
                                                                           acaGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGa                                                                          >  2:8391026/1‑150 (MQ=255)
                                                                               gTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAAc                                                                        <  2:3367403/148‑1 (MQ=255)
                                                                                    aCGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTa                                                                 >  2:7021287/1‑150 (MQ=255)
                                                                                     cGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCTCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTaa                                                                <  1:825674/150‑1 (MQ=255)
                                                                                     cGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTaa                                                                <  2:7155559/150‑1 (MQ=255)
                                                                                     cGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTaa                                                                <  1:10331894/150‑1 (MQ=255)
                                                                                      ggAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAg                                                               <  1:8496097/150‑1 (MQ=255)
                                                                                      ggAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAg                                                               >  1:5709997/1‑150 (MQ=255)
                                                                                       gAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAgg                                                              >  1:4581090/1‑149 (MQ=255)
                                                                                       gAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAgg                                                              >  1:631978/1‑149 (MQ=255)
                                                                                        aataatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaa                                                             <  2:4909275/150‑1 (MQ=255)
                                                                                        aataatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaa                                                             <  2:8969288/150‑1 (MQ=255)
                                                                                        aataatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaa                                                             <  2:9050395/150‑1 (MQ=255)
                                                                                        aataatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACAGATTTAAGaa                                                             <  1:7234539/150‑1 (MQ=255)
                                                                                        aataatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaa                                                             <  2:1326485/150‑1 (MQ=255)
                                                                                           aatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaacaa                                                          >  1:8584028/1‑150 (MQ=255)
                                                                                           aatCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGaacaa                                                          >  2:1488000/1‑150 (MQ=255)
                                                                                             tCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACg                                                        >  2:3548520/1‑150 (MQ=255)
                                                                                             tCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACg                                                        >  1:8076227/1‑150 (MQ=255)
                                                                                             tCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACg                                                        <  1:10559853/150‑1 (MQ=255)
                                                                                             tCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATAGTTAAGAACAACg                                                        >  1:9519853/1‑150 (MQ=255)
                                                                                                ggAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGttt                                                     <  2:2197104/150‑1 (MQ=255)
                                                                                                   aCGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCa                                                  <  1:7553034/150‑1 (MQ=255)
                                                                                                      cTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGa                                               <  2:3801564/150‑1 (MQ=255)
                                                                                                       ttttGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGAAGGat                                              >  1:1643565/1‑150 (MQ=255)
                                                                                                        tttGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGata                                             <  2:10555567/150‑1 (MQ=255)
                                                                                                           gCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCg                                          <  2:7385052/150‑1 (MQ=255)
                                                                                                            cATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGt                                         <  1:11142359/150‑1 (MQ=255)
                                                                                                              tCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCa                                       >  1:12257631/1‑150 (MQ=255)
                                                                                                                 gCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGa                                    >  2:12106096/1‑150 (MQ=255)
                                                                                                                 gCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGa                                    >  2:10960629/1‑150 (MQ=255)
                                                                                                                 gCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGa                                    >  2:3152325/1‑150 (MQ=255)
                                                                                                                    aaTCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATgg                                 >  1:2372133/1‑150 (MQ=255)
                                                                                                                        gCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGAc                             <  1:3762243/150‑1 (MQ=255)
                                                                                                                         cTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGAcc                            <  1:5260115/150‑1 (MQ=255)
                                                                                                                           gAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCa                          >  2:9302062/1‑150 (MQ=255)
                                                                                                                            aaTATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCaa                         >  2:3678849/1‑150 (MQ=255)
                                                                                                                            aaTATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCaa                         >  2:4591591/1‑150 (MQ=255)
                                                                                                                             atatCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAAc                        <  1:6647525/150‑1 (MQ=255)
                                                                                                                              tatCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAAcg                       >  2:2316322/1‑150 (MQ=255)
                                                                                                                               atCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAAcgc                      <  2:5179823/150‑1 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  <  1:11551714/150‑1 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  >  1:8773220/1‑150 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  <  1:9451539/150‑1 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  >  2:2571296/1‑150 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  >  1:4033141/1‑150 (MQ=255)
                                                                                                                                   ttCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATc                  <  1:2599108/150‑1 (MQ=255)
                                                                                                                                      tcGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtc               <  2:11702175/150‑1 (MQ=255)
                                                                                                                                      tcGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtc               <  2:6783949/150‑1 (MQ=255)
                                                                                                                                      tcGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtc               >  2:5962960/1‑150 (MQ=255)
                                                                                                                                      tcGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtc               <  1:5419662/150‑1 (MQ=255)
                                                                                                                                       cgcgGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcg              <  2:6351284/150‑1 (MQ=255)
                                                                                                                                         cgGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTtcgtc            >  2:7863185/1‑150 (MQ=255)
                                                                                                                                             aGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGccc        >  1:11467663/1‑150 (MQ=255)
                                                                                                                                                aGGTCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCa     <  2:1829487/150‑1 (MQ=255)
                                                                                                                                                 ggCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGACCAGATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGccccccg    >  1:2178852/1‑148 (MQ=255)
                                                                                                                                                  gCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGt   >  2:5368759/1‑150 (MQ=255)
                                                                                                                                                  gCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGt   >  2:2393009/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGtt  >  1:10215209/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGGACGAATAACGATAAGTTTGCCAACAACGCTCATCAGTTTCGGAGGCAAACCTTTCATCTAACCTCCGTAGAACATATTAAAGAACAACGTTTGCAGGATATCGTCACTATGGTGACCCAACGCGAACTTCGTCGCCACCAgat  >  2:567529/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCCTGGATAACGTTTATCCCAGTCACGCAACATGTCAGCAATCACCTGACGTTGCAGGTTAGGCTGTGAACCGCACAGGTTGCACGGAATAATCGGGAACGCTTTTGCATCGGCAAATCGCTGAATATCTTTCTCGCGGCAGTAGGCCAGCGGACGAATAACGATATGTTTGCCATCATCGCTCATCAGTTTCGGAGGCATACCTTTCATCTTACCGCCGTAGAACATATTTAAGAACAACGTTTGCAGGATATCGTCACGATGGTGACCCAACGCGATCTTCGTCGCCCCCAGTT  >  NC_000913/1411237‑1411533

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