Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 828,135 C→T 6.2% G283R (GGG→AGG)  ybhF ← putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913828,1350CT 6.2% 203.6 / 2.3 64G283R (GGG→AGG) ybhFputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (29/31);  new base T (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCACCGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTGGTCAT  >  NC_000913/827993‑828279
                                                                                                                                              |                                                                                                                                                
aTCACCGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTc                                                                                                                                               <  5:102071/146‑1 (MQ=255)
aTCACCGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCgg                                                                                                                                             >  3:163782/1‑148 (MQ=255)
aTCACCGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCAGTCATCCCTTc                                                                                                                                               >  6:102071/1‑146 (MQ=255)
    ccGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAACCGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCTGTCGGCATgg                                                                                                                                        >  3:45601/1‑149 (MQ=255)
    ccGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCTGTCGGCATgg                                                                                                                                        >  3:40234/1‑149 (MQ=255)
           gCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGaa                                                                                                                                    <  3:384124/146‑1 (MQ=255)
           gCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAata                                                                                                                                  <  2:211553/148‑1 (MQ=255)
           gCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAata                                                                                                                                  <  4:331654/148‑1 (MQ=255)
               ggTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTc                                                                                                                             <  7:79697/149‑1 (MQ=255)
                          tACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCct                                                                                                                  >  7:346326/1‑149 (MQ=255)
                                 tGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCtttttt                                                                                                           >  6:294080/1‑149 (MQ=255)
                                  gTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCtttttt                                                                                                           >  3:330559/1‑148 (MQ=255)
                                   tAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTg                                                                                                          <  4:163782/148‑1 (MQ=255)
                                     atatTGCGCCCCGCGGCGATTCCGAGGGTCCGGCACCGACCAGCAAAACAATACACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTTTTTTCCGGCATCCCGTCGGCATGACGAATATCGTCTGGTGTGGCCTCTTTTTTTAg                                                                                                        <  4:274502/148‑1 (MQ=255)
                                     atatTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAgg                                                                                                       <  5:354362/149‑1 (MQ=255)
                                      tatTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGa                                                                                                      <  5:294080/149‑1 (MQ=255)
                                           cgcCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAg                                                                                                  <  8:346326/148‑1 (MQ=255)
                                                     cgATTCCGAGGTTCCGGCACCGACCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGAtt                                                                                       <  5:281231/149‑1 (MQ=255)
                                                        ttCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGgt                                                                                                                         >  2:171758/1‑112 (MQ=255)
                                                        ttCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGgt                                                                                                                         <  1:171758/112‑1 (MQ=255)
                                                         tCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATc                                                                                                                                <  2:406791/104‑1 (MQ=255)
                                                         tCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATc                                                                                                                                >  1:406791/1‑104 (MQ=255)
                                                          ccGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTccc                                                                                   <  3:8347/148‑1 (MQ=255)
                                                          ccGAGGGTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGGTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGGCTGGTGGGGCCTCTTTTTTGAGGATCAGACGTACCGATTTccc                                                                                   <  6:414622/148‑1 (MQ=255)
                                                           cGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCTGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCt                                                                                 <  4:40234/149‑1 (MQ=255)
                                                           cGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCt                                                                                 <  5:213710/149‑1 (MQ=255)
                                                           cGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCt                                                                                 >  1:162683/1‑149 (MQ=255)
                                                           cGAGGCTCCGGCACCGCCAAGCAAATCAATAACCGCAGCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCTGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCt                                                                                 <  4:45601/149‑1 (MQ=255)
                                                            gAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGCTTTCCGGCATCCCGTCGGCATGCCGAATATCGTCTGGTGTGGCCTCGTTTTTGAGGATCAGACGTACCGATTTCCCCt                                                                                 >  6:113566/1‑148 (MQ=255)
                                                            gAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTg                                                                                <  4:330559/149‑1 (MQ=255)
                                                              ggTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGaa                                                                              <  7:263271/149‑1 (MQ=255)
                                                                 tCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCa                                                                           >  7:244321/1‑149 (MQ=255)
                                                                     gCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCg                                                                                                                                              <  1:15001/78‑1 (MQ=255)
                                                                     gCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCg                                                                                                                                              >  2:15001/1‑78 (MQ=255)
                                                                           cccAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCGTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATTCCGTCGCt                                                                 >  7:279803/1‑149 (MQ=255)
                                                                                   aTCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTCCCGATTTCCCCTGAATCATGCCGTCGCTGCCCTGCg                                                         >  3:218450/1‑149 (MQ=255)
                                                                                   aTCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCg                                                         <  7:25521/149‑1 (MQ=255)
                                                                                   aTCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCg                                                         <  8:279803/149‑1 (MQ=255)
                                                                                     cAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCgg                                                        >  8:123677/1‑148 (MQ=255)
                                                                                       aTAAACGCATCTTCAAAACGCGGCGTACTTTCCTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTAGAGATTTCCCCTGAATCATGCAGTCGCTGACCTGCGGCAg                                                     >  5:86853/1‑149 (MQ=255)
                                                                                         aaaCGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGtt                                                   >  8:218425/1‑149 (MQ=255)
                                                                                          aaCGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGttt                                                  <  2:171153/149‑1 (MQ=255)
                                                                                           aCGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGttt                                                  >  4:126677/1‑148 (MQ=255)
                                                                                               aTCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGAGAATATGGTCTGGTGTGGCCTCTTTTTTGAGGATCAGCCGTACCGAATTCCCCTGAATCATGCCGTCGCTGACCTTCGGCAGTTTCAAgg                                             >  2:377648/1‑149 (MQ=255)
                                                                                                   tCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGaa                                                                              <  7:305833/112‑1 (MQ=255)
                                                                                                   tCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGaa                                                                              >  8:305833/1‑112 (MQ=255)
                                                                                                    cAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGgcgc                                          <  7:123677/147‑1 (MQ=255)
                                                                                                      aaaCGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAg                                              <  5:289418/141‑1 (MQ=255)
                                                                                                      aaaCGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAg                                              >  6:289418/1‑141 (MQ=255)
                                                                                                                tAGTTTCGTTGATGTTGATTTCCGGCATCACGTCGGCATGGCGAATATAGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCAGAATCATTCCGTCGCTGAACTGCGGCAGTATCAAGGCGCGTTGAAACAGTTTg                            >  4:378359/1‑149 (MQ=255)
                                                                                                                   tttCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCgg                         >  3:326280/1‑149 (MQ=255)
                                                                                                                       gttgatgttgatTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTg                      <  3:378359/148‑1 (MQ=255)
                                                                                                                       gttgatgttgatTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGc                     >  7:285132/1‑149 (MQ=255)
                                                                                                                              ttgatTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCgt               <  3:126677/148‑1 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCg                          >  5:114515/1‑133 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCg                          <  6:114515/133‑1 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGa          <  4:376385/149‑1 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGa          <  8:244321/149‑1 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGa          >  7:8122/1‑149 (MQ=255)
                                                                                                                                  tttCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCACTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGa          >  1:64586/1‑149 (MQ=255)
                                                                                                                                        gCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTGGt     >  4:175134/1‑148 (MQ=255)
                                                                                                                                           tCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGCTTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACATTTTTCGGTTGCCCTCTTGTGGACTGGtcat  >  6:155409/1‑148 (MQ=255)
                                                                                                                                           tCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTg       >  1:273049/1‑143 (MQ=255)
                                                                                                                                           tCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTg       <  2:273049/143‑1 (MQ=255)
                                                                                                                                             ccGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGAc                                                              <  1:67967/86‑1 (MQ=255)
                                                                                                                                             ccGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGAc                                                              >  2:67967/1‑86 (MQ=255)
                                                                                                                                             ccGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTGGtcat  >  1:372483/1‑146 (MQ=255)
                                                                                                                                              |                                                                                                                                                
ATCACCGTCTCGCCGGGTGTGCCTTCTACCGTATGTAATATTGCGCCCAGCGGCGATTCCGAGGTTCCGGCACCGCCCAGCAAATCAATAAACGCATCTTCAAAACGCGGCGTAGTTTCGTTGATGTTGATTTCCGGCATCCCGTCGGCATGGCGAATATCGTCTGGTGTGGCCTCTTTTTTGAGGATCAGACGTACCGATTTCCCCTGAATCATGCCGTCGCTGACCTGCGGCAGTTTCAAGGCGCGTTGCAACAGTTTGCGGTTGCCCTCGTGTGGACTGGTCAT  >  NC_000913/827993‑828279

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: ‑