Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,206,179 T→G 8.9% A26A (GCA→GCC ttdR ← DNA‑binding transcriptional activator Dan

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,1790TG 8.9% 257.7 / 16.9 90A26A (GCA→GCCttdRDNA‑binding transcriptional activator Dan
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base G (2/6);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 2.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

ATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCA  >  NC_000913/3206030‑3206326
                                                                                                                                                     |                                                                                                                                                   
aTTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCt                                                                                                                                                     >  2:971262/1‑150 (MQ=255)
   tCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcg                                                                                                                                                  >  1:919073/1‑150 (MQ=255)
   tCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcg                                                                                                                                                  >  2:644137/1‑150 (MQ=255)
   tCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcg                                                                                                                                                  >  2:649603/1‑150 (MQ=255)
     cAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTgcggc                                                                                                                                                >  2:413204/1‑150 (MQ=255)
         gCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGAAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTa                                                                                                                                            <  2:62725/150‑1 (MQ=255)
          cATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                           <  1:886449/150‑3 (MQ=255)
          cATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                           <  1:324024/150‑3 (MQ=255)
          cATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttg                                                                                                                                           <  2:1116794/150‑3 (MQ=255)
            tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttggg                                                                                                                                         <  1:1145424/150‑5 (MQ=255)
            tGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttggg                                                                                                                                         <  1:371214/150‑5 (MQ=255)
             ggtCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGGTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGCCGCAGCGGCGACGCttgggc                                                                                                                                        <  1:1622780/148‑6 (MQ=255)
             gTTTGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAc                                                                                                                                        >  2:1491378/1‑150 (MQ=255)
             gTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAc                                                                                                                                        >  1:1274366/1‑150 (MQ=255)
               tCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAAcac                                                                                                                                      >  1:903086/1‑148 (MQ=255)
                 gTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTccc                                                                                                                                    >  2:924485/1‑150 (MQ=255)
                  tAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCg                                                                                                                                   >  1:141177/1‑150 (MQ=255)
                  tAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCg                                                                                                                                   >  1:994604/1‑150 (MQ=255)
                   aacaacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGc                                                                                                                                  <  1:772495/150‑1 (MQ=255)
                     caacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCtg                                                                                                                                <  2:466784/150‑1 (MQ=255)
                     caacGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAACAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCtg                                                                                                                                >  1:1406492/1‑150 (MQ=255)
                         gTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                            <  2:423561/150‑1 (MQ=255)
                         gTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGaa                                                                                                                            <  2:919073/150‑1 (MQ=255)
                          ttGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAAc                                                                                                                           <  1:1312120/150‑1 (MQ=255)
                            ggCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGa                                                                                                                         >  1:191082/1‑150 (MQ=255)
                             gCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAt                                                                                                                        <  1:644137/150‑1 (MQ=255)
                             gCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAt                                                                                                                        <  1:649603/150‑1 (MQ=255)
                              ccACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAtt                                                                                                                       <  2:842058/150‑1 (MQ=255)
                               cACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAttt                                                                                                                      >  2:9032/1‑150 (MQ=255)
                                aCTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACCAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTc                                                                                                                     <  1:504424/150‑1 (MQ=255)
                                 cTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCt                                                                                                                    <  2:711556/150‑1 (MQ=255)
                                 cTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCt                                                                                                                    >  1:773143/1‑150 (MQ=255)
                                  tttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTa                                                                                                                   >  1:234246/1‑150 (MQ=255)
                                   ttCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAc                                                                                                                  >  1:255955/1‑150 (MQ=255)
                                    tcagtcagCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTAcc                                                                                                                 <  1:762182/150‑1 (MQ=255)
                                          agCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTAcc                                                                                                           <  1:183540/150‑1 (MQ=255)
                                             gcCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGc                                                                                                        <  2:1630404/150‑1 (MQ=255)
                                              ccACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCa                                                                                                       >  1:1025399/1‑150 (MQ=255)
                                                aCGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAg                                                                                                     >  2:427683/1‑150 (MQ=255)
                                                   ccgcGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCt                                                                                                  <  2:1268180/150‑1 (MQ=255)
                                                      cgAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTg                                                                                               >  2:1400494/1‑150 (MQ=255)
                                                       gagcggagcgGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTgg                                                                                              <  2:234246/150‑1 (MQ=255)
                                                        agcggagcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGc                                                                                             <  1:971262/150‑1 (MQ=255)
                                                             agcggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAg                                                                                        <  1:584552/150‑1 (MQ=255)
                                                                ggTTCCGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGcc                                                                                     >  2:821269/1‑150 (MQ=255)
                                                                ggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGGAAGTCTTTGGCTAAAGGcc                                                                                     >  1:989125/1‑150 (MQ=255)
                                                                ggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGcc                                                                                     >  2:298816/1‑150 (MQ=255)
                                                                ggTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGcc                                                                                     <  2:467428/150‑1 (MQ=255)
                                                                  ttCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAg                                                                                   >  2:12291/1‑150 (MQ=255)
                                                                   tCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGc                                                                                  <  1:924485/150‑1 (MQ=255)
                                                                    cagcagCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTaaa                                                                                         >  1:1301953/1‑142 (MQ=255)
                                                                    cagcagCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTaaa                                                                                         <  2:1301953/142‑1 (MQ=255)
                                                                    cagcagCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCt                                                                                 <  1:829906/150‑1 (MQ=255)
                                                                     agcagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTa                                                                                >  2:977059/1‑150 (MQ=255)
                                                                      gcagcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAt                                                                               >  1:1608081/1‑150 (MQ=255)
                                                                         gcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCa                                                                            <  1:1555963/150‑1 (MQ=255)
                                                                         gcGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCa                                                                            <  1:530011/150‑1 (MQ=255)
                                                                            ttgttgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCa                                                                         <  2:1425747/150‑1 (MQ=255)
                                                                                tgCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGt                                                                     >  1:22246/1‑150 (MQ=255)
                                                                                  ccAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTa                                                                   <  1:415216/150‑1 (MQ=255)
                                                                                   cAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTaa                                                                  >  2:235824/1‑150 (MQ=255)
                                                                                    aGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                                 >  2:1173208/1‑150 (MQ=255)
                                                                                    aGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAt                                                                 <  1:1235485/150‑1 (MQ=255)
                                                                                     gCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAAtt                                                                >  1:120554/1‑150 (MQ=255)
                                                                                       gaaTTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCACGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTAGTCAGCATCGTTAATTAt                                                              >  2:765076/3‑150 (MQ=255)
                                                                                         aTTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATcc                                                            >  2:1533998/1‑150 (MQ=255)
                                                                                          ttttCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCg                                                           <  1:9032/150‑1 (MQ=255)
                                                                                           tttCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGc                                                          <  2:1608081/150‑1 (MQ=255)
                                                                                             tCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAg                                                        <  1:1629088/150‑1 (MQ=255)
                                                                                              cAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGt                                                       <  1:156130/150‑1 (MQ=255)
                                                                                               aaGGATTTGGATCCGCTTGGTGACAAATGTAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGAGTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGtt                                                      >  1:678569/1‑150 (MQ=255)
                                                                                                aGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTtg                                                     <  2:1574889/150‑1 (MQ=255)
                                                                                                    tttGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGAt                                                 <  1:1308741/150‑1 (MQ=255)
                                                                                                       ggATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAg                                              <  1:1556067/150‑1 (MQ=255)
                                                                                                                 ggTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAt                                    >  2:1135168/1‑150 (MQ=255)
                                                                                                                 ggTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAt                                    >  2:1591572/1‑150 (MQ=255)
                                                                                                                 ggTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTAt                                    >  2:1586087/1‑150 (MQ=255)
                                                                                                                   tGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCACAGTTGTGATAAGCGCAGTGTAttt                                  >  1:1293734/1‑150 (MQ=255)
                                                                                                                     aCAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCg                                <  2:501741/150‑1 (MQ=255)
                                                                                                                      cAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTACAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGc                               >  2:850101/1‑150 (MQ=255)
                                                                                                                          tGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCaaaa                           <  2:1406492/150‑1 (MQ=255)
                                                                                                                            cAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACa                         >  1:493968/1‑150 (MQ=255)
                                                                                                                            cAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACa                         >  2:388818/1‑150 (MQ=255)
                                                                                                                               gCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTg                      <  2:255955/150‑1 (MQ=255)
                                                                                                                                cGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACACCCGCTGTGCACGAATTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGc                     <  2:1530930/150‑1 (MQ=255)
                                                                                                                                     gCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTc                <  1:427683/150‑1 (MQ=255)
                                                                                                                                     gCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTc                >  1:960446/1‑150 (MQ=255)
                                                                                                                                     gCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTc                <  1:1267355/150‑1 (MQ=255)
                                                                                                                                       tttaGCGCCGCAGCGGCGACGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTcca              >  1:927428/4‑149 (MQ=255)
                                                                                                                                              ccgcagcGGCGACGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTg       >  1:248349/2‑150 (MQ=255)
                                                                                                                                                   gcTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCa  <  1:977059/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ATTTCCAACGCATGTTCGTAACAACGTTGGCCACTTTCAGTCAGCGCCACGCCGCGAGCGGAGCGGTTCAGCAGCGTTGTTGCCAGCGTATTTTCAAGGATTTGGATCCGCTTGGTGACAAATGCAGGCGTCTGCCCAAGCGTCGCCGCTGCGGCGCTAAAACTCCCGCTGTGAACGATTTCTACCAGTACCTGCAAGTCTTTGGCTAAAGGCCAGCTATTCAGCATCGTTAATTATCCGCAGTTGTGATAAGCGCAGTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCA  >  NC_000913/3206030‑3206326

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