Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,179,562 C→T 5.6% G228G (GGC→GGT)  rplA → 50S ribosomal subunit protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,179,5620CT 5.6% 238.4 / 4.5 91G228G (GGC→GGT) rplA50S ribosomal subunit protein L1
Reads supporting (aligned to +/- strand):  ref base C (49/34);  new base T (3/2);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGC  >  NC_000913/4179414‑4179709
                                                                                                                                                    |                                                                                                                                                   
aCTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCGCTGCTGGTTGCGCTGAAAACCGCAAATCCGACGCAGGCGAAAGGCGGGTACATCCAGAAAGTTAGGATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGc                                                                                                                                                     <  2:445106/149‑1 (MQ=255)
   ttGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTg                                                                                                                                                  >  2:596293/1‑149 (MQ=255)
    tGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGAGTGGACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGGGTTGCAGTTGACCAGGCTGGCCTGa                                                                                                                                                 >  2:535522/1‑149 (MQ=255)
      aCGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGCCTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCAGCTCCACCCCCATGGGTGAAGGTGTTGGAGTTGCCCAGGCTGGCCTGAgg                                                                                                                                               >  2:672540/1‑148 (MQ=255)
           gACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTc                                                                                                                                          >  1:632259/1‑149 (MQ=255)
           gACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGCAGAAATGCGTGTACATCAAGAAAGTTAGAATCTCCACCCCCATGGGTGCAGGGGTTGAAGTTGCCCACGCTGGCCTGTGAGCTTc                                                                                                                                          >  1:291958/1‑149 (MQ=255)
            aCAAACTGAAAGAACACCTGGAAGCTCTGCTGGTTGCGCTGAAAAACGCAAAACCGACTCAGGCGAAAGGCGGGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGGTGCAGTTGACCAGGCTGGCCTGAGCGCTTCt                                                                                                                                         <  2:1094823/149‑1 (MQ=255)
            aCAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTACt                                                                                                                                         >  1:131456/1‑149 (MQ=255)
              aaaCTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGt                                                                                                                                       >  1:204347/1‑149 (MQ=255)
               aaCTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTa                                                                                                                                      <  1:734413/149‑1 (MQ=255)
                 cTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGTGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTaaa                                                                                                                                    <  1:491897/149‑1 (MQ=255)
                  tGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAAc                                                                                                                                   <  1:501783/149‑1 (MQ=255)
                  tGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAAc                                                                                                                                   >  1:125828/1‑149 (MQ=255)
                  tGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAAc                                                                                                                                   <  1:13278/149‑1 (MQ=255)
                  tGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAAc                                                                                                                                   >  1:246080/1‑149 (MQ=255)
                      cgaaaACCTGGAAGCGCTGCTGGTGGCGCTGCAAACAGCAAACCCGACTCAGGCGAACGGCGTGTACATCCAGAAAGTTCGCATCTCCACCACCATGGGTGCAGGTGTGGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAAt                                                                                                                               <  1:182993/148‑1 (MQ=255)
                      agaaaACCTGGAAGCTCTGCTGGTTGCGCTGAAAACAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCGGAGCGCTTCTGTAAACTAAt                                                                                                                               <  2:545854/149‑1 (MQ=255)
                      agaaaACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAAt                                                                                                                               >  2:46006/1‑149 (MQ=255)
                      agaaaACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAAt                                                                                                                               >  2:848564/1‑149 (MQ=255)
                      agaaaACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGACTGAGCGCTTCTGTAAACTAAt                                                                                                                               >  1:407764/1‑149 (MQ=255)
                       gaaaACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAAGTAATg                                                                                                                              >  2:976293/1‑149 (MQ=255)
                            ccTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCttt                                                                                                                         <  1:446988/149‑1 (MQ=255)
                            ccTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCttt                                                                                                                         <  1:596293/149‑1 (MQ=255)
                             cTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTa                                                                                                                        >  2:270961/1‑149 (MQ=255)
                                gAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAATGCGTGTACATAAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGt                                                                                                                     <  2:122398/149‑1 (MQ=255)
                                  aGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCCTGAACATCAAGAAAGTTAGCATCTCCACCACCATCGCTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTgg                                                                                                                   >  2:1071667/1‑149 (MQ=255)
                                     tctGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCg                                                                                                                <  1:857912/149‑1 (MQ=255)
                                         ctgGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGa                                                                                                            >  1:423063/1‑149 (MQ=255)
                                         ctgGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATATCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGAGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGa                                                                                                            >  1:879634/1‑149 (MQ=255)
                                                cgcTGAAAAATGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGGAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTc                                                                                                     <  2:331612/149‑1 (MQ=255)
                                                  cTGAAAAAACCAAAACCGACTCAGGCGAAAGGCGTGGACATCAAGACAGGAAGCAGCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAATCTAATGCCTTTACGTGGGCGGGGATTTTGTCTa                                                                                                   <  1:707471/149‑1 (MQ=255)
                                                     aaaaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGCGCGCTTCTGTAATCTAATGCCTTTACGTGGGCGGTGATTTTGTCTACaa                                                                                                >  2:121556/1‑149 (MQ=255)
                                                       aaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              >  1:1391/1‑149 (MQ=255)
                                                       aaaaGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              >  1:115847/1‑149 (MQ=255)
                                                           gCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTAc                                                                                          >  2:224153/1‑149 (MQ=255)
                                                            cAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTAcc                                                                                         >  1:159469/1‑149 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGTCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              >  2:203302/1‑143 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGTCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              <  1:793897/143‑1 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGTCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              <  1:203302/143‑1 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGTCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATc                                                                                              >  2:793897/1‑143 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTAccc                                                                                        >  1:421900/1‑149 (MQ=255)
                                                             aaaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTAccc                                                                                        >  1:1161189/1‑149 (MQ=255)
                                                              aaaCCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTAcccc                                                                                       >  2:353220/1‑149 (MQ=255)
                                                                          gCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGCCCAGGCTGGCCGGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGc                                                                           >  2:214053/1‑149 (MQ=255)
                                                                              aaGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTaa                                                                       <  2:125828/149‑1 (MQ=255)
                                                                               aGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGGGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAAt                                                                      <  2:501492/149‑1 (MQ=255)
                                                                                ggCGTGTACATCAACAACGTTAGCATCCCCCCCACCATGGGTGCAGGGGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATg                                                                     <  2:700430/149‑1 (MQ=255)
                                                                                 gCGTGAACAGCACGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGc                                                                    <  1:121614/149‑1 (MQ=255)
                                                                                  cGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCa                                                                   >  2:841791/1‑149 (MQ=255)
                                                                                    tgtACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTAAAATCTTACCCCCACGTATAATGCTTAATGCAGa                                                                 <  2:660384/149‑1 (MQ=255)
                                                                                    tgtACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGAGTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGa                                                                 >  2:155020/1‑149 (MQ=255)
                                                                                     gtACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGAc                                                                >  1:840568/1‑149 (MQ=255)
                                                                                      tACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACg                                                               >  1:1001832/1‑149 (MQ=255)
                                                                                       aCATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGt                                                              <  2:421900/149‑1 (MQ=255)
                                                                                       aCATCAAGAAAGTTAGCATCGCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGt                                                              <  2:347968/149‑1 (MQ=255)
                                                                                         aTCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGtat                                                            >  1:508759/1‑149 (MQ=255)
                                                                                            aaGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATc                                                         >  2:29959/1‑149 (MQ=255)
                                                                                             aGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATcc                                                        <  1:630845/149‑1 (MQ=255)
                                                                                             aGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATcc                                                        >  2:739813/1‑149 (MQ=255)
                                                                                             aGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATcc                                                        >  2:394924/1‑149 (MQ=255)
                                                                                               aaaGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCga                                                      >  1:974241/1‑149 (MQ=255)
                                                                                                aaGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCgag                                                     <  1:705499/149‑1 (MQ=255)
                                                                                                      gCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTc                                               >  1:1131032/1‑149 (MQ=255)
                                                                                                        aTCTCCACCAACATGGGTGCAGGTGGTGCATCTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTTCGTGGGCGGGGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATTCAGACGTATATCCGAGATATTCgg                                             <  2:1069484/149‑1 (MQ=255)
                                                                                                            ccccccccaTGGGTGCAGGGGTTGCAGTTGACCAGGCGGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGGCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTtg                                         <  1:464164/143‑1 (MQ=255)
                                                                                                            ccaccacgcTGGGTGCAGGTGTTGCACTTGACCAGGCTGGCCTGCGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATAATCGGGTtg                                         <  2:122942/149‑1 (MQ=255)
                                                                                                            ccaccaccaTGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGAGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTtg                                         <  1:224153/149‑1 (MQ=255)
                                                                                                              cccaccaTGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTATTGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGATTAATGCAGACGTATATCCGAGATATTCGGGTtgtg                                       <  2:369858/148‑1 (MQ=255)
                                                                                                              accaccaTGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTtgtg                                       >  2:535935/1‑149 (MQ=255)
                                                                                                               ccaccaTGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTgg                                      <  1:1018344/149‑1 (MQ=255)
                                                                                                                       ggTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAggcgg                              >  1:225211/1‑149 (MQ=255)
                                                                                                                       ggTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAggcgg                              <  2:632259/149‑1 (MQ=255)
                                                                                                                         tGCAGGTGTTGCAGTTGACCAGGCTGGCCTCCGCGCGTCTGTAATCTAATGCCTTTACGAGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTAGAAGGCTTAATGCAGACGTATATCCGAGCTATTCGGGTGGTGTCAAGGCGCCa                            >  1:769172/1‑149 (MQ=255)
                                                                                                                            aGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACt                         <  2:115847/149‑1 (MQ=255)
                                                                                                                            aGGTGTTGCACTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACTATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGAAATTCGGGTTGTGGCAAGGCGGCAACt                         <  1:1119297/149‑1 (MQ=255)
                                                                                                                              gtgtTGCAGTTGACCAGGCTGTACTGAGCGCTTCTGTAAACTAATGCCTTTACTTGTGCGGTGATTTTGTCTACAATCTTACACCCACGTAAAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACtga                       >  2:1052241/1‑149 (MQ=255)
                                                                                                                               tgtTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACtgag                      >  1:1108310/1‑149 (MQ=255)
                                                                                                                                gtTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACtgagt                     <  1:739813/149‑1 (MQ=255)
                                                                                                                                gtTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCATTTACGTGGGCGGTGATTATGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACtgagt                     <  2:346505/149‑1 (MQ=255)
                                                                                                                                 ttGCAGGTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTATGTCTACTATCTTACCCCCACGTATAATGCTGAGTGCAGTCGTAGAGCCGAGATATTCGGGATGTGGCAAGGCGGCAACtgagtg                    >  2:657011/1‑149 (MQ=255)
                                                                                                                                   gcggTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACtgagtgag                  >  1:168006/4‑149 (MQ=255)
                                                                                                                                   gCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGACATATTCGGGTTGTGCCAAGGCGGCAACtgagtgat                  >  1:814432/1‑148 (MQ=255)
                                                                                                                                       ttGACAAGGCTGGTCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGAGTTTGTCTACAATCTTACCCCCACGAATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGc              >  2:395505/1‑149 (MQ=255)
                                                                                                                                         gACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCa            >  1:150501/1‑149 (MQ=255)
                                                                                                                                         gACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCAACGTATAATGCTTAATCCAGACGTATATCCGAGAAATTCGGGTTGTGGCAATGCGGCAACTGAGTGAGTAGcaa            >  2:7848/1‑147 (MQ=255)
                                                                                                                                            cAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGTCAAGGCGGCAACTGAGTGAGTCGCCAGGa         >  2:743607/1‑149 (MQ=255)
                                                                                                                                            cAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAggt         >  2:1072035/1‑148 (MQ=255)
                                                                                                                                            cAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAggg         >  2:847448/1‑148 (MQ=255)
                                                                                                                                              ggCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGc       >  1:557349/1‑149 (MQ=255)
                                                                                                                                              ggCTGGACTTAGCGCTTATGTAAATTAATTCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCCCGTATAATGCTTAATGCAGACGTATATCCGAGAGATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGc       <  2:879634/149‑1 (MQ=255)
                                                                                                                                                gtGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACTTATATCCTAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCAt     <  1:1071667/148‑1 (MQ=255)
                                                                                                                                                   gCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGc  >  1:1056933/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTTTGACGCTGACAAACTGAAAGAAAACCTGGAAGCTCTGCTGGTTGCGCTGAAAAAAGCAAAACCGACTCAGGCGAAAGGCGTGTACATCAAGAAAGTTAGCATCTCCACCACCATGGGTGCAGGTGTTGCAGTTGACCAGGCTGGCCTGAGCGCTTCTGTAAACTAATGCCTTTACGTGGGCGGTGATTTTGTCTACAATCTTACCCCCACGTATAATGCTTAATGCAGACGTATATCCGAGATATTCGGGTTGTGGCAAGGCGGCAACTGAGTGAGTCGCCAGGAGCATAGC  >  NC_000913/4179414‑4179709

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