Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,704,843 A→C 8.2% M316R (ATG→AGG)  dppB ← dipeptide/heme ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,704,8430AC 8.2% 161.4 / 3.0 61M316R (ATG→AGG) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (40/15);  new base C (2/4);  total (42/19)
Fisher's exact test for biased strand distribution p-value = 6.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

ACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCTGCAGGCCG  >  NC_000913/3704697‑3704989
                                                                                                                                                  |                                                                                                                                                  
aCTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAatcatc                                                                                                                                                  <  1:451046/149‑1 (MQ=255)
  tCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCAtcgt                                                                                                                                                >  1:619872/1‑149 (MQ=255)
   ccggtgaCGGGGTCATCGGCACTGGTGCGCTAATCGCTTTATTTTTAGTTACCTGTGAGCTGATGTGCTTTTACTTCTTATGACGACTACGCGGGTTCACCACGCCGTACAGCAGAGCGACCAGCAAGTTGATGAGGATAATCAtcgtc                                                                                                                                               <  2:305843/143‑1 (MQ=255)
   ccggTAACGGGGTCACCGGCACCGGTGCGCTAATCATTTGATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCAtcgtc                                                                                                                                               <  2:415332/146‑1 (MQ=255)
     tGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGATGATTATCCtcgtcaa                                                                                                                                             <  1:350334/149‑3 (MQ=255)
     tGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGATGATTATCCtcgtcaa                                                                                                                                             <  1:202758/149‑3 (MQ=255)
     ggTAACGGGGTCATCCGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACACGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGATGATTATCCtcgtcaa                                                                                                                                             <  1:21370/148‑3 (MQ=255)
        gaCGGGNTCATCGGCACCCGTGCGCTACTCACTTTATTTTCAGGAACCTGTGGCATGATGTCCCCTTACTTCTTATGACGAATACGCGTGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGccac                                                                                                                                          <  2:241066/148‑1 (MQ=255)
         aCGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCACCCAGTTGACGAGGATAATCCTCGCCGccccc                                                                                                                                         >  2:743774/1‑149 (MQ=255)
          cGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCGAGTTGACGAGGATAATCATCGTCGccacca                                                                                                                                        <  1:407826/149‑1 (MQ=255)
             ggTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTGCTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCGGCAAGGTGAGGAGGATAATCATCGTCGACACCCGCa                                                                                                                                     >  2:206689/1‑149 (MQ=255)
             ggTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCCGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCGACCAgaa                                                                                                                                     >  2:667743/1‑147 (MQ=255)
               tCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTTTGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAAGATAATCATCGTCGCCACCAGCAAt                                                                                                                                   >  1:480848/1‑149 (MQ=255)
               tCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGGTAATCATCGTCGCCACCAGCAAt                                                                                                                                   >  1:46526/1‑149 (MQ=255)
                   cGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATAcgc                                                                                                                               >  1:580318/1‑149 (MQ=255)
                   cGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATAcgc                                                                                                                               >  1:474719/1‑149 (MQ=255)
                    ggCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAGTACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCCTCGTCGCCACCAGCAATAcgca                                                                                                                              >  1:659865/1‑148 (MQ=255)
                            tgaacTAATCACTTTATTTCCAGGAACCCGGGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCa                                                                                                                      <  2:390575/145‑1 (MQ=255)
                                cTAATCACTTTATTTTCCGTATCCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTAc                                                                                                                  <  1:501752/149‑1 (MQ=255)
                                        tttATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGt                                                                                                          >  1:366462/1‑149 (MQ=255)
                                        tttATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGt                                                                                                          >  1:604632/1‑149 (MQ=255)
                                           actttCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGCCCAGCAAGTTGACGCGGATAATCCTCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTcgc                                                                                                       <  2:561260/147‑1 (MQ=255)
                                             tttCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCgcgg                                                                                                     >  1:176412/1‑149 (MQ=255)
                                                   tAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGc                                                                                               >  1:414818/1‑149 (MQ=255)
                                                   tAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGc                                                                                               >  2:352630/1‑149 (MQ=255)
                                                        tgtgACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCCCGGCGTTGCAGTGc                                                                                          >  1:406798/1‑149 (MQ=255)
                                                        tgtgACATGATGTCCCCTTACTTCTTAAGACGAATAAGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCAACGTCGCCACCAGCAAAACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCaaggc                                                                                          >  1:66496/1‑145 (MQ=255)
                                                                 atgTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGGGTTGCAGGGCGtcaatcaa                                                                                 >  1:435634/1‑149 (MQ=255)
                                                                     ccccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCGGTGCGTCAATCAACCAg                                                                             >  2:416229/1‑149 (MQ=255)
                                                                     ccccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGGCAATCAACCAg                                                                             >  2:610333/1‑149 (MQ=255)
                                                                      cccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAGTCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTCGAGTGCGTCAATCAACCAGc                                                                            >  1:227987/1‑149 (MQ=255)
                                                                      cccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGc                                                                            >  2:605885/1‑149 (MQ=255)
                                                                      cccTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGc                                                                            >  1:295833/1‑149 (MQ=255)
                                                                                ttATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGCAAGTCGCGGCGTTGCAGGGCGTCAATCAACCAGCGCCCCAGAcc                                                                  >  1:83036/1‑149 (MQ=255)
                                                                                ttATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGAcc                                                                  >  1:228530/1‑149 (MQ=255)
                                                                                ttATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGAcc                                                                  >  2:685869/1‑149 (MQ=255)
                                                                                ttATGACGAATACGCGGGCTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCACAACGCCGCCCTGCCCGCCCGCATAGTCGGGGCGTTGCAATCCATCAATCAACCATCATACCAAAcc                                                                  >  1:228528/1‑149 (MQ=255)
                                                                                  aTGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTCCCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCCGACCgg                                                                >  2:49843/1‑149 (MQ=255)
                                                                                   tGACCAATACGCGGGTTCCCCACGCCGTACAGCAGATCGACCAGCAAGTGGACGGGTGTAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCggg                                                               <  1:105620/149‑1 (MQ=255)
                                                                                     aCGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGcc                                                             >  2:534094/1‑149 (MQ=255)
                                                                                     aCGAATACGCGGGATCACCACGCCGTACAGCAGAGCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACCACCGGATAGTCGCGGCGTTGCAGTGCGACAATCTCCCAGCGCCCCAGCCCGGGcc                                                             >  2:688201/1‑149 (MQ=255)
                                                                                            cgcgGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGagaa                                                      <  1:661378/149‑1 (MQ=255)
                                                                                             gcgGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGagaag                                                     >  2:309774/1‑149 (MQ=255)
                                                                                                ggTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATg                                                  <  1:694454/149‑1 (MQ=255)
                                                                                                 gTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGGTgg                                                 >  1:46899/1‑149 (MQ=255)
                                                                                                   tCACCACGCCGTACAGCAGATCGACCGGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGTGTTGCAGTGCGTCAATCAACCAGCGGCCCAGACCGGGCCACGCGAAGATGGtg                                               >  1:87519/1‑148 (MQ=255)
                                                                                                       cacGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATCGTCGCGGCGCTGCAGTCCGTCAATCAACCAGCGTCCCGGACCGGGCCACGAGAAGAGGGGTTCgg                                           >  2:347882/1‑149 (MQ=255)
                                                                                                         cGCCGTACAGCAGATCGCCCAGCAAGCTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTGGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTc                                         <  2:187830/149‑1 (MQ=255)
                                                                                                            cGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGa                                      >  2:97072/1‑149 (MQ=255)
                                                                                                                         gtccAGTAAGTTGGCGAGGATAATCACCGTCGCCACCAGCAATACGCCGCCCTGCATTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCa                         <  1:57078/147‑1 (MQ=255)
                                                                                                                         gACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGCGAACATCGTTTCCGTCAGAATCGCCCCCGcaa                         >  1:364824/1‑147 (MQ=255)
                                                                                                                         gACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGccg                         >  1:364745/1‑148 (MQ=255)
                                                                                                                         gACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCa                         >  2:200983/1‑149 (MQ=255)
                                                                                                                         gACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCa                         <  1:759694/149‑1 (MQ=255)
                                                                                                                         gACCAGAAAGTTGACGAGGATAATCACCGTCGCCACCAGCAATACGCCACCCTGCAGTACCGGATATTAGCGACGTTGCAGTGCGTCAACCAAGAAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCa                         >  2:337147/1‑149 (MQ=255)
                                                                                                                            cAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCa                      <  1:459038/149‑1 (MQ=255)
                                                                                                                             aGCAAGTCGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAAAAGATGGTTTCGGTCAGAACCGCCCCCGCCCGCaa                     >  2:482843/1‑149 (MQ=255)
                                                                                                                              gCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAAt                    >  2:632705/1‑149 (MQ=255)
                                                                                                                                      aCGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCCGCAATGTTCCCAc            >  1:139598/1‑149 (MQ=255)
                                                                                                                                        gAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCt          <  2:292898/149‑1 (MQ=255)
                                                                                                                                        gAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCt          <  2:722652/149‑1 (MQ=255)
                                                                                                                                              aaTCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCTGCAGGc    >  1:351881/1‑149 (MQ=255)
                                                                                                                                               atcatcGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCTGCAGGcc   <  1:605885/149‑1 (MQ=255)
                                                                                                                                                tcatcGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCTGCAGGCCg  >  1:434606/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACTCCTGTAACGGGGTCATCGGCACCGGTGCGCTAATCACTTTATTTTCAGTAACCTGTGACATGATGTCCCCTTACTTCTTATGACGAATACGCGGGTTCACCACGCCGTACAGCAGATCGACCAGCAAGTTGACGAGGATAATCATCGTCGCCACCAGCAATACGCCGCCCTGCACTACCGGATAGTCGCGGCGTTGCAGTGCGTCAATCAACCAGCGTCCCAGACCGGGCCACGAGAAGATGGTTTCGGTCAGAATCGCCCCCGCCAGCAATGTTCCCACCTGCAGGCCG  >  NC_000913/3704697‑3704989

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