Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,332,822 G→A 6.0% intergenic (‑87/+40) yhbE ← / ← rpmA EamA family inner membrane putative transporter/50S ribosomal subunit protein L27

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,332,8220GA 6.0% 409.2 / 7.6 131intergenic (‑87/+40)yhbE/rpmAEamA family inner membrane putative transporter/50S ribosomal subunit protein L27
Reads supporting (aligned to +/- strand):  ref base G (64/59);  new base A (5/3);  total (69/62)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CGCCCCCCAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACG  >  NC_000913/3332673‑3332970
                                                                                                                                                     |                                                                                                                                                    
cGCCCCCCAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACg                                                                                                                                                      <  2:526896/150‑1 (MQ=255)
  ccccccAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACgtg                                                                                                                                                    <  2:872225/150‑1 (MQ=255)
     cccAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg                                                                                                                                                 >  2:1875486/1‑150 (MQ=255)
     cccAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg                                                                                                                                                 <  1:243457/150‑1 (MQ=255)
     cccAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTg                                                                                                                                                 >  1:1475022/1‑150 (MQ=255)
       cAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGAGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCg                                                                                                                                               >  1:248307/1‑150 (MQ=255)
           aaaTTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACACGTTGCGGGGc                                                                                                                                           <  2:494088/150‑1 (MQ=255)
             aTTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCtt                                                                                                                                         <  1:1544306/150‑1 (MQ=255)
              ttGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGGGTTGCGGGGCttt                                                                                                                                        >  1:842790/1‑150 (MQ=255)
                gCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                                                                      >  1:1505465/1‑150 (MQ=255)
                gCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                                                                      <  1:1733079/150‑1 (MQ=255)
                gCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                                                                      <  1:1757364/150‑1 (MQ=255)
                gCTGTGGGGAGCGCCAAAAGAATGCCCATGCCTGCCGGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCa                                                                                                                                      <  1:833051/150‑1 (MQ=255)
                 cTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCGTGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCAt                                                                                                                                     >  1:1272771/1‑150 (MQ=255)
                  tgtgGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATc                                                                                                                                    >  1:1623930/1‑150 (MQ=255)
                   gtggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATcc                                                                                                                                   >  1:1825214/1‑150 (MQ=255)
                   gtggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATcc                                                                                                                                   <  2:1351503/150‑1 (MQ=255)
                    tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCg                                                                                                                                  <  2:1458847/150‑1 (MQ=255)
                    tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCg                                                                                                                                  <  1:749106/150‑1 (MQ=255)
                    tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCACATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGGTGCGGGGCTTTCATCCg                                                                                                                                  <  1:1005790/150‑1 (MQ=255)
                    tggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGAAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCg                                                                                                                                  >  1:1656734/1‑150 (MQ=255)
                     ggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGt                                                                                                                                 >  2:1861342/1‑150 (MQ=255)
                     ggtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGt                                                                                                                                 >  1:216587/1‑150 (MQ=255)
                      gtgAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGAGGGGCTTTCATCCGtt                                                                                                                                >  2:1001780/1‑150 (MQ=255)
                          gcgcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCAAACCGTATTCCCTGCAAAGCGCATTTGCGGAGTACCGTCGCCCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTGCCg                                                                                                                            <  1:357615/150‑1 (MQ=255)
                           cgacAAACGAATGCCAATGCCAGCCTGCTGCTTCATACCGTATACCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCAAGCGCATTTTTACCGTTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCgg                                                                                                                           <  1:965089/147‑1 (MQ=255)
                            gcCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCggg                                                                                                                          <  1:799228/150‑1 (MQ=255)
                              cAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAc                                                                                                                        >  2:1109548/1‑150 (MQ=255)
                                      tgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaaa                                                                                                                >  2:722238/1‑150 (MQ=255)
                                      tgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaaa                                                                                                                <  2:220932/150‑1 (MQ=255)
                                      tgccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGaaaaa                                                                                                                >  1:491276/1‑150 (MQ=255)
                                       gccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                                                                               <  2:43269/150‑1 (MQ=255)
                                       gccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                                                                               <  1:1269579/150‑1 (MQ=255)
                                       gccaaggccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAAc                                                                                                               <  2:322770/150‑1 (MQ=255)
                                        ccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACt                                                                                                              >  2:2015606/1‑150 (MQ=255)
                                        ccaatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGAAGCGAAAAACt                                                                                                              >  2:1711953/1‑150 (MQ=255)
                                          aatgccTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTa                                                                                                            <  1:642337/150‑1 (MQ=255)
                                            tgcctgccTTCTGCTTCATACCGTATTCCCTTCAAAGCGCATTTGCTGATTAACGTCTCACCATGCGCATTTTTACCGTTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCTGTGCTTTCATCCGTTACCGGGACGCGAAAAACTTAtt                                                                                                          <  1:1811138/150‑1 (MQ=255)
                                              cctgcctgCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCa                                                                                                        <  1:1266884/150‑1 (MQ=255)
                                                 gcctgcTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGGGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCt                                                                                                     >  1:1572002/1‑150 (MQ=255)
                                                    tgctgcTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGATACCGGGACGCGAAAAACTTATTCAGCTTCg                                                                                                  >  2:470584/1‑150 (MQ=255)
                                                     gctgctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGAAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGa                                                                                                 <  1:727231/150‑1 (MQ=255)
                                                        gctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGc                                                                                              >  2:111599/1‑150 (MQ=255)
                                                         ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCt                                                                                             <  2:568420/150‑1 (MQ=255)
                                                         ctTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCt                                                                                             <  2:926297/150‑1 (MQ=255)
                                                          ttCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTg                                                                                            >  1:17242/1‑150 (MQ=255)
                                                             aTACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATa                                                                                         >  1:178235/1‑150 (MQ=255)
                                                                 cGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTATCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAtt                                                                                     <  2:1475022/150‑1 (MQ=255)
                                                                 cGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAAtt                                                                                     >  2:1123635/1‑150 (MQ=255)
                                                                    aTTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTAc                                                                                  <  1:98776/150‑1 (MQ=255)
                                                                       cccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGATAAACTTATTCAGCTTCGATGCTGATAAATTTACGGt                                                                               >  2:962186/1‑150 (MQ=255)
                                                                        ccTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGtt                                                                              >  2:1781109/1‑150 (MQ=255)
                                                                         cTGCAAAGCGCATTTGCGGATTAACGTCGCCCCATGCGCATGTTTACCGGTTATCGTATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCTTCCGTTACCGGGACGCGAAAAACTTATTCAGCTTAGATGCTGAGAAATTTACGGttt                                                                             >  2:250889/1‑150 (MQ=255)
                                                                           gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                                                                           >  2:628091/1‑150 (MQ=255)
                                                                           gCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGttttt                                                                           >  1:142780/1‑150 (MQ=255)
                                                                              aaGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCgg                                                                        <  2:1024692/150‑1 (MQ=255)
                                                                                  gcATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCt                                                                    >  1:42182/1‑150 (MQ=255)
                                                                                   cATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCtt                                                                   >  2:1035405/1‑150 (MQ=255)
                                                                                      ttGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTaa                                                                >  2:642936/1‑150 (MQ=255)
                                                                                          ggATTAACGTCGCACCATGCGCATTTTTACCGGTTGTCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTc                                                            <  2:54803/150‑1 (MQ=255)
                                                                                           gATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCg                                                           >  2:646624/1‑150 (MQ=255)
                                                                                              tAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAt                                                        <  2:1620232/150‑1 (MQ=255)
                                                                                                 cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc                                                     <  2:1783506/150‑1 (MQ=255)
                                                                                                 cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTc                                                     <  2:1656734/150‑1 (MQ=255)
                                                                                                 cGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACATATTCAGCTTCGATGCTGATAAATTTACGGTTATTCGGGCCTTTAACTTCGAATTTc                                                     >  1:47005/1‑150 (MQ=255)
                                                                                                  gTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCa                                                    >  2:1570510/1‑150 (MQ=255)
                                                                                                   tCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGTGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATTCTGTTAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCAc                                                   <  1:1998115/150‑1 (MQ=255)
                                                                                                     gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACtt                                                 <  2:590970/150‑1 (MQ=255)
                                                                                                     gCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACtt                                                 >  2:420193/1‑150 (MQ=255)
                                                                                                       aCCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTa                                               >  1:687314/1‑150 (MQ=255)
                                                                                                         cctGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTAcc                                             <  1:1522232/148‑1 (MQ=255)
                                                                                                          aTGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCg                                            >  1:1657757/1‑150 (MQ=255)
                                                                                                                      aCCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCa                                >  1:1083246/1‑150 (MQ=255)
                                                                                                                        cGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCGTTAGCaaa                              >  2:380728/1‑150 (MQ=255)
                                                                                                                          gTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCTAACa                            >  1:1219473/1‑150 (MQ=255)
                                                                                                                           ttATCGAATGTAAAAAGCCCCGAAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACag                           <  2:79609/150‑1 (MQ=255)
                                                                                                                              acGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCTGGAAGAGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTAAGCAAACagaga                        <  2:1468189/149‑2 (MQ=255)
                                                                                                                               cGATTGTATACATCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAATCAGAgtg                       <  1:346580/150‑1 (MQ=255)
                                                                                                                               cGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAgtt                       >  2:1916468/1‑149 (MQ=255)
                                                                                                                                gAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCttt                                    >  1:1771249/1‑136 (MQ=255)
                                                                                                                                 aaTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGtgtg                     <  2:1572002/150‑1 (MQ=255)
                                                                                                                                  aTGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTgg                    <  1:854518/150‑1 (MQ=255)
                                                                                                                                   tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGt                   >  2:1291323/1‑150 (MQ=255)
                                                                                                                                   tGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGt                   <  2:283228/150‑1 (MQ=255)
                                                                                                                                    gTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAAGTTCACTTTACCGTCTGCTTTAGCAATCAGAGTGTGGTc                  >  1:537796/1‑150 (MQ=255)
                                                                                                                                      tgaaaGCCCCGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCAc                >  2:1448954/3‑150 (MQ=255)
                                                                                                                                      aaaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACg               <  1:366564/150‑1 (MQ=255)
                                                                                                                                       aaaaGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACg               >  2:1954863/1‑150 (MQ=255)
                                                                                                                                          aGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGTTAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCGGCTTTAGCAAACAGAGTGTGGTCACGAcc            >  1:1928053/1‑150 (MQ=255)
                                                                                                                                           gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCg           <  2:178235/150‑1 (MQ=255)
                                                                                                                                           gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATGTCACTTTACCGTCTGATTTAGCAAACAGAGTGTGGTCACGcccc           >  2:1284947/1‑149 (MQ=255)
                                                                                                                                           gCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTAAGCTTCGATGCTGATAAATTTACGGATTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACTAccc           >  2:221397/1‑149 (MQ=255)
                                                                                                                                            ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATTAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGc          >  1:836209/1‑150 (MQ=255)
                                                                                                                                            ccccGCAACACGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGc          >  1:450918/1‑150 (MQ=255)
                                                                                                                                             cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCa         >  2:1685312/1‑150 (MQ=255)
                                                                                                                                             cccGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCa         >  1:755757/1‑150 (MQ=255)
                                                                                                                                                gCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAAcc      <  2:1944169/150‑1 (MQ=255)
                                                                                                                                                gCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAAcc      <  1:2015606/150‑1 (MQ=255)
                                                                                                                                                    ggtgtTTCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACg  <  1:963510/149‑1 (MQ=255)
                                                                                                                                                    cGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACg  >  2:1529785/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCCCCCAGCAAATTGCTGTGGTGAGCGCCAAAAGAATGCCAATGCCTGCCTGCTGCTTCATACCGTATTCCCTGCAAAGCGCATTTGCGGATTAACGTCGCACCATGCGCATTTTTACCGGTTATCGAATGTAAAAAGCCCCGCAACGTGTTGCGGGGCTTTCATCCGTTACCGGGACGCGAAAAACTTATTCAGCTTCGATGCTGATAAATTTACGGTTTTTCGGGCCTTTAACTTCGAATTTCACTTTACCGTCTGCTTTAGCAAACAGAGTGTGGTCACGACCGCAACCTACG  >  NC_000913/3332673‑3332970

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