Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,229,322 C→G 5.5% intergenic (+257/+116) yohP → / ← dusC uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,229,3220CG 5.5% 623.7 / 22.5 183intergenic (+257/+116)yohP/dusCuncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand):  ref base C (86/87);  new base G (5/5);  total (91/92)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 183 total aligned reads displayed.

AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAG  >  NC_000913/2229177‑2229469
                                                                                                                                                 |                                                                                                                                                   
  ggggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgt                                                                                                                                                >  4:420437/1‑149 (MQ=32)
   gggCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtgtg                                                                                                                                               <  4:567026/149‑1 (MQ=32)
    ggcggcCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGTGGTGCCCTCGGCTTATCCCGTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGGCCTCAGCCCACATCTATGTgg                                                                                                                                              <  2:31085/149‑1 (MQ=255)
         ccAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTga                                                                                                                                         <  1:247319/149‑2 (MQ=17)
           aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                                                                                                                                       <  4:273707/149‑1 (MQ=32)
           aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCACGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGccc                                                                                                                                       >  1:985741/1‑149 (MQ=17)
           aaTCGCCGCCGCCCCCTGCACCCCCGGGCTCGGGCGAACAAAATCGCCGCTGCTCGGTGCTCTCGGCTTATCCCTTACTGCTACCGGGTCGGGCGCGAGTTAACATTCCTGTAAAACGCGCCCTCATCCCACATCCATGTGCGCTGccc                                                                                                                                       >  2:1324361/1‑149 (MQ=11)
            aTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGcccc                                                                                                                                      >  4:1295839/1‑149 (MQ=32)
             tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                                                                                                                                     >  2:973257/1‑149 (MQ=32)
             tCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCg                                                                                                                                     <  3:325404/149‑1 (MQ=32)
               gccgccgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGc                                                                                                                                   <  1:44143/149‑1 (MQ=32)
                    cgccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa                                                                                                                              <  4:1061/149‑1 (MQ=32)
                     gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAg                                                                                                                             <  3:1092386/149‑1 (MQ=32)
                     gccCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAg                                                                                                                             >  2:116112/1‑149 (MQ=32)
                      cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg                                                                                                                            >  1:209391/1‑149 (MQ=32)
                      cccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAgg                                                                                                                            <  3:1119101/149‑1 (MQ=32)
                       ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg                                                                                                                           <  2:432064/149‑1 (MQ=32)
                       ccccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAggg                                                                                                                           >  4:1259594/1‑149 (MQ=32)
                        cccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGa                                                                                                                          >  2:727001/1‑149 (MQ=32)
                         ccTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCCGCCCCGGCCTTCAGGGaa                                                                                                                         <  4:670878/149‑1 (MQ=17)
                              aCCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCt                                                                                                                    >  3:294532/1‑149 (MQ=32)
                                ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCa                                                                                                                              <  4:1130861/137‑1 (MQ=32)
                                ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCg                                                                                                                  <  1:783843/149‑1 (MQ=32)
                                ccccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCg                                                                                                                  >  1:1253502/1‑149 (MQ=32)
                                 cccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCgg                                                                                                                 >  3:183914/1‑149 (MQ=32)
                                  ccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTgag                                                                                                                                        <  1:1451904/125‑3 (MQ=12)
                                  ccGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATGGATGTGGGCTgag                                                                                                                                        >  2:1451904/1‑123 (MQ=12)
                                        tctGGCGAACAAAATCGCCGCTGCTCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAGCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAAtttt                                                                                                          <  3:1447234/149‑1 (MQ=21)
                                          tGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGa                                                                                                        <  2:600315/149‑1 (MQ=34)
                                           ggCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGAc                                                                                                       <  1:169674/149‑1 (MQ=35)
                                               aaCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCCTGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGcca                                                                                                   >  4:718/1‑149 (MQ=21)
                                               aaCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGcca                                                                                                   <  3:525725/149‑1 (MQ=35)
                                                aCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGccac                                                                                                  >  2:1247212/1‑149 (MQ=35)
                                                  aaaaTCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGccacca                                                                                                >  1:499704/1‑149 (MQ=34)
                                                      tCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAaaca                                                                                            >  4:1439846/1‑149 (MQ=37)
                                                             tGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGc                                                                                     <  1:1039697/149‑1 (MQ=255)
                                                                gcgGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGc                                                                                  <  3:1330451/149‑1 (MQ=255)
                                                                   gTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTa                                                                               <  2:764714/149‑1 (MQ=255)
                                                                    tGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTAt                                                                              >  2:1199609/1‑149 (MQ=255)
                                                                    tGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGACCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTCCAGGGAACGCCTCGGCAATTTTGACGCCAACAAACAACCGTGCGGCCTAt                                                                              <  4:801593/149‑1 (MQ=255)
                                                                        cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                          >  1:114951/1‑149 (MQ=255)
                                                                        cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                          <  2:497315/149‑1 (MQ=255)
                                                                        cTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGAt                                                                          >  1:1013294/1‑149 (MQ=255)
                                                                             cTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAaga                                                                     <  4:147380/149‑1 (MQ=255)
                                                                              ttATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAagag                                                                    <  1:973257/149‑1 (MQ=255)
                                                                                  cccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTaa                                                                >  1:1211964/1‑149 (MQ=255)
                                                                                  cccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTaa                                                                <  3:162918/149‑1 (MQ=255)
                                                                                   ccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAac                                                               >  4:194245/1‑149 (MQ=255)
                                                                                   ccTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAac                                                               >  2:1257881/1‑149 (MQ=255)
                                                                                     ttACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAacac                                                             >  2:571099/1‑149 (MQ=255)
                                                                                      tACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAAcaca                                                            <  2:209391/149‑1 (MQ=255)
                                                                                           cTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTc                                                       <  3:853108/149‑1 (MQ=255)
                                                                                           cTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTc                                                       >  2:532039/1‑149 (MQ=255)
                                                                                            tACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCa                                                      >  3:775929/1‑149 (MQ=255)
                                                                                             aCCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCaa                                                     >  1:1213351/1‑149 (MQ=255)
                                                                                                      ggCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCa                                            <  2:353994/149‑1 (MQ=255)
                                                                                                       gcgcGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCAc                                           <  2:1083430/149‑1 (MQ=255)
                                                                                                       gcgcGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCGTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGCCCTATTGATAAAGGGCTAACCCATTGTCAAAAAACATCAc                                           >  3:377644/1‑149 (MQ=255)
                                                                                                           gAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                       <  1:34959/149‑1 (MQ=255)
                                                                                                            aGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATgg                                      >  1:874613/1‑149 (MQ=255)
                                                                                                             ggTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCTCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGt                                     >  2:1432755/1‑149 (MQ=255)
                                                                                                              gTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGtt                                    >  1:560624/1‑149 (MQ=255)
                                                                                                                aaCATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCACAAAACATCACTATGGtttt                                  >  4:744409/1‑149 (MQ=255)
                                                                                                                   aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCTCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTa                               <  1:1432755/149‑1 (MQ=255)
                                                                                                                   aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACa                                               <  4:833204/133‑1 (MQ=255)
                                                                                                                   aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTa                               >  4:921857/1‑149 (MQ=255)
                                                                                                                   aTCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTa                               <  2:316427/149‑1 (MQ=255)
                                                                                                                     cccTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGGCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGATCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                             <  4:350631/149‑1 (MQ=255)
                                                                                                                     cccTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTaga                             <  2:847823/149‑1 (MQ=255)
                                                                                                                        tGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGtt                          >  2:169883/1‑149 (MQ=255)
                                                                                                                        tGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGtt                          >  4:712055/1‑149 (MQ=255)
                                                                                                                         gTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGttt                         >  4:1359811/1‑149 (MQ=255)
                                                                                                                            aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtctc                      <  2:1405243/149‑1 (MQ=255)
                                                                                                                            aaaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTtctc                      >  4:456628/1‑149 (MQ=255)
                                                                                                                             aaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                     >  1:1086026/1‑149 (MQ=255)
                                                                                                                             aaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                     <  2:629681/149‑1 (MQ=255)
                                                                                                                             aaCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCg                     <  2:1253502/149‑1 (MQ=255)
                                                                                                                              aCGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGa                    >  4:534263/1‑149 (MQ=255)
                                                                                                                                   ggggcAGCCCACATGGATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                       <  4:1118265/121‑1 (MQ=255)
                                                                                                                                     cTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                       <  1:137295/91‑1 (MQ=255)
                                                                                                                                     cTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAAt             >  4:619781/1‑149 (MQ=255)
                                                                                                                                       cAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGGTTTTTAGAGTTTCTCGATATCAATTg           >  3:21511/1‑149 (MQ=255)
                                                                                                                                        aGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCa                                            >  1:59840/1‑115 (MQ=255)
                                                                                                                                        aGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc          <  1:532039/149‑1 (MQ=255)
                                                                                                                                        aGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc          <  3:636269/149‑1 (MQ=255)
                                                                                                                                        aGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGc          <  3:815687/149‑1 (MQ=255)
                                                                                                                                          cccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCt        <  3:712055/149‑1 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                       >  3:452907/1‑85 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                       <  4:452907/85‑1 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATaaa                                                                       <  3:536245/85‑1 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACAtt                                                          >  1:114704/1‑98 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTc                                                       <  4:595959/101‑1 (MQ=255)
                                                                                                                                           ccACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATg                                       <  3:1121896/117‑1 (MQ=255)
                                                                                                                                            cacaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGa      >  3:79750/1‑149 (MQ=255)
                                                                                                                                             acaTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAATCATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGaa     >  3:1266949/1‑149 (MQ=255)
                                                                                                                                               aTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAGCACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATAGCAATTGCCTGAATa   >  4:1236639/1‑149 (MQ=255)
                                                                                                                                               aTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATa   >  3:518930/1‑149 (MQ=255)
                                                                                                                                               aTCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATa   <  2:408031/149‑1 (MQ=255)
                                                                                                                                                tCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAg  >  4:1234727/1‑149 (MQ=255)
                                                                                                                                                tCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAg  <  1:416355/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATGTGGGCTGCCCCGGCCTTCAGGGAACGCCTCGGCAATTTTGACGCCACCAAACAACCGTGCGGCCTATTGATAAAGAGCTAACACATTGTCAAAAAACATCACTATGGTTTTTTAGAGTTTCTCGATATCAATTGCCTGAATAG  >  NC_000913/2229177‑2229469

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