Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,953,532 C→T 5.3% Q31* (CAG→TAG)  cmoB → tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,953,5320CT 5.3% 241.7 / 2.3 75Q31* (CAG→TAG) cmoBtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base C (42/29);  new base T (2/2);  total (44/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTG  >  NC_000913/1953387‑1953679
                                                                                                                                                 |                                                                                                                                                   
ggTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTTCCCGCGCAGATTGCTAACTGGCAgc                                                                                                                                                  >  3:135662/1‑149 (MQ=255)
     cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGag                                                                                                                                             >  2:396966/1‑149 (MQ=255)
     cAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGa                                                                                                                                              <  6:411875/148‑1 (MQ=255)
      aGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGag                                                                                                                                             >  8:30423/1‑148 (MQ=255)
             tAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGagcagcag                                                                                                                                       <  4:280347/147‑1 (MQ=255)
              aaCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCAc                                                                                                                                    >  3:37025/1‑149 (MQ=255)
               aCTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACg                                                                                                                                   >  7:120112/1‑149 (MQ=255)
                    ggTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGgtt                                                                                                                               <  7:207510/148‑1 (MQ=255)
                               ggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTg                                                                                                                   <  4:37025/149‑1 (MQ=255)
                               ggCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGGTGTTTAAGCAGTg                                                                                                                   <  2:96275/149‑1 (MQ=255)
                                gCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTg                                                                                                                   <  1:396966/148‑1 (MQ=255)
                                 cATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGt                                                                                                                 >  8:104072/1‑149 (MQ=255)
                                  aTTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGt                                                                                                                 >  5:412964/1‑148 (MQ=255)
                                  aTTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTc                                                                                                                <  3:750/149‑1 (MQ=255)
                                   ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTc                                                                                                                >  7:314128/1‑148 (MQ=255)
                                   ttAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTc                                                                                                                >  6:199711/1‑148 (MQ=255)
                                      aaaGCAGATGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGAGGTCaaa                                                                                                             <  8:342601/148‑1 (MQ=255)
                                         gCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCg                                                                                                         >  5:281404/1‑149 (MQ=255)
                                                  gCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATTTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTg                                                                                                      >  3:210878/1‑143 (MQ=255)
                                                  gCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATTTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTg                                                                                                      <  4:210878/143‑1 (MQ=255)
                                                  gCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGttt                                                                                                <  8:314128/149‑1 (MQ=255)
                                                   cTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTc                                                                                               >  3:239177/1‑149 (MQ=255)
                                                    tGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTc                                                                                               >  7:169936/1‑148 (MQ=255)
                                                     gCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGGTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTa                                                                                             >  4:227364/1‑149 (MQ=255)
                                                      cATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCt                                                                                              <  8:169936/147‑1 (MQ=255)
                                                       aTGATCGACTTTTGTAACTTTTATTCTCTGATTGCCAAAAATAATCTTTCACACTGGCTAGAAACGCTGCCCGCGAAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTAGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTAcc                                                                                           >  5:163276/1‑149 (MQ=255)
                                                       aTGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTAcc                                                                                           >  7:140414/1‑149 (MQ=255)
                                                        tGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCt                                                                                          >  6:155573/1‑149 (MQ=255)
                                                         gATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCt                                                                                          >  7:428628/1‑148 (MQ=255)
                                                         gATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCt                                                                                          >  7:108661/1‑148 (MQ=255)
                                                          aTCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTg                                                                                         >  8:3302/1‑148 (MQ=255)
                                                                tttGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAgcagcagca                                                                                                                                     <  8:304098/98‑1 (MQ=255)
                                                                tttGGTAAATTTTATTCTCTGATTGCCAAAAATCATCTTTCAAACTGGCTCGAAACGCTTCCCGCGCAGATTCCTAACTGGCACAGCGAGCAGCAGCAAGGGTTGTTTAACCAGAGCTCAAAAGCTCTGGATTTTCTACCTGAAATTa                                                                                   >  6:217940/1‑148 (MQ=255)
                                                                  tGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTaa                                                                                  <  8:140414/147‑1 (MQ=255)
                                                                    gTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAAcc                                                                               <  3:237430/148‑1 (MQ=255)
                                                                       aCTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTAt                                                                           <  5:155573/149‑1 (MQ=255)
                                                                         ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCg                                                                         >  6:93659/1‑149 (MQ=255)
                                                                         ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCg                                                                         >  1:412319/1‑149 (MQ=255)
                                                                         ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCg                                                                         >  2:23898/1‑149 (MQ=255)
                                                                         ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCg                                                                         >  2:105823/1‑149 (MQ=255)
                                                                         ttttATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACAGTATCg                                                                         >  8:28195/1‑149 (MQ=255)
                                                                             aTTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCt                                                                      <  7:3302/148‑1 (MQ=255)
                                                                               tctcTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGa                                                                   >  7:14140/1‑149 (MQ=255)
                                                                                     aTTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTAtt                                                             <  1:105823/149‑1 (MQ=255)
                                                                                       tGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGAGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTg                                                            >  5:207841/1‑148 (MQ=255)
                                                                                         ccAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGTAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAAGTAAACCGTATCGTCTTGATTTATTGCa                                                          >  7:27948/1‑148 (MQ=255)
                                                                                         ccAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGTAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAACTAAACCGTATCGTCTTGATTTATTGCa                                                          >  7:21574/1‑148 (MQ=255)
                                                                                           aaaaaTCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGAGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAg                                                       <  6:207841/149‑1 (MQ=255)
                                                                                            aaaaTCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGc                                                      <  7:30423/149‑1 (MQ=255)
                                                                                             aaaTCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCg                                                     <  5:217940/149‑1 (MQ=255)
                                                                                                  atcTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCTAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACg                                                <  5:199711/149‑1 (MQ=255)
                                                                                                  atcTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACg                                                >  2:397669/1‑149 (MQ=255)
                                                                                                   tcTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACgg                                               >  1:177164/1‑149 (MQ=255)
                                                                                                     tttCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACgg                                               >  1:366129/1‑147 (MQ=255)
                                                                                                     tttCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACgg                                               >  1:371971/1‑147 (MQ=255)
                                                                                                     tttCACACAGGCTCGAAACGCTACCCGCTCAGATTGCAAACTCGCAGCGCGAGCAGCAGCACGGGATGTTTAAGCAGTTGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGGATCGTCTTGAGTTATTGCATAGCGTTACGGCa                                             >  4:204584/1‑149 (MQ=255)
                                                                                                            cTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAg                                                       >  7:126442/1‑132 (MQ=255)
                                                                                                            cTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAg                                                       <  8:126442/132‑1 (MQ=255)
                                                                                                             tGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCg                                      >  6:52111/1‑148 (MQ=255)
                                                                                                               gCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGaa                                    <  6:281404/148‑1 (MQ=255)
                                                                                                                   gAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGc                                >  7:190820/1‑148 (MQ=255)
                                                                                                                      aCGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACt                            >  7:354647/1‑149 (MQ=255)
                                                                                                                          tGCCCGCGCAGATTGCTGACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAgc                        >  2:16739/1‑149 (MQ=255)
                                                                                                                            cccGCGCAGATTGCTAACTGGTAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAACTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAgcgc                      <  8:21574/149‑1 (MQ=255)
                                                                                                                            cccGCGCAGATTGCTAACTGGTAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAAATAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAgcgc                      <  8:27948/149‑1 (MQ=255)
                                                                                                                             ccGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCa                     <  7:85826/149‑1 (MQ=255)
                                                                                                                              cgcgCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATCGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAg                    >  1:309370/1‑149 (MQ=255)
                                                                                                                              cgcgCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAg                    <  5:93659/149‑1 (MQ=255)
                                                                                                                               gcgcAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTAt                                                              <  4:367895/106‑1 (MQ=255)
                                                                                                                               gcgcAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTAt                                                              >  3:367895/1‑106 (MQ=255)
                                                                                                                                         cTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTaaa         <  4:51671/149‑1 (MQ=255)
                                                                                                                                         cTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTaaa         <  8:108661/149‑1 (MQ=255)
                                                                                                                                          tAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAAc        >  6:43584/1‑149 (MQ=255)
                                                                                                                                                gCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCAtt   >  8:175491/1‑148 (MQ=255)
                                                                                                                                                gCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCAt    <  1:16739/147‑1 (MQ=255)
                                                                                                                                                 cAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTg  >  6:213782/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGTTCCAGTGCTTTAACTTTGGTTCACTGGTGGCATTAAAAGCAGAGGACGCTGCATGATCGACTTTGGTAACTTTTATTCTCTGATTGCCAAAAATCATCTTTCACACTGGCTCGAAACGCTGCCCGCGCAGATTGCTAACTGGCAGCGCGAGCAGCAGCACGGGTTGTTTAAGCAGTGGTCAAATGCGGTGGAGTTTCTACCTGAAATTAAACCGTATCGTCTTGATTTATTGCATAGCGTTACGGCAGAAAGCGAAGAGCCACTGAGCGCAGGGCAAATTAAACGCATTG  >  NC_000913/1953387‑1953679

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