Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 870,063 T→C 5.1% I218T (ATA→ACA)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913870,0630TC 5.1% 243.8 / 2.8 78I218T (ATA→ACA) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

AGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCCGTCAATTATG  >  NC_000913/869917‑870211
                                                                                                                                                  |                                                                                                                                                    
aGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATaa                                                                                                                                                    <  6:122811/149‑1 (MQ=255)
 gcgcTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAAc                                                                                                                                                   <  4:324458/149‑1 (MQ=255)
     tGGAATAATATCGCATGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGTAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCt                                                                                                                                                 <  6:285624/147‑1 (MQ=255)
       gAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGc                                                                                                                                              <  5:48494/148‑1 (MQ=255)
       gAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCg                                                                                                                                             >  1:275214/1‑149 (MQ=255)
          aaaTATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTcc                                                                                                                                          >  6:360711/1‑149 (MQ=255)
          aaaTATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTcc                                                                                                                                          <  4:551/149‑1 (MQ=255)
             tatGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGt                                                                                                                                       <  8:419141/149‑1 (MQ=255)
                   aaGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGa                                                                                                                                 <  8:368345/149‑1 (MQ=255)
                   aaGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGa                                                                                                                                 >  8:303987/1‑149 (MQ=255)
                    aGGAGATTTGTTTTTATAGGGTGGGTACCGGAAAGTATGAACTGGATACCTGGATTCAGACCTATTTTGTGAAGGTGAAAAAATTCGCTGGTTACTGGCAGCCAGGATTGCCCTAACTGGACAGCATAACCTGGCGTCCGGTGGCGGAt                                                                                                                                <  5:266042/149‑1 (MQ=255)
                    aGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGAt                                                                                                                                <  1:257276/149‑1 (MQ=255)
                    aGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGGTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGAt                                                                                                                                <  5:369851/149‑1 (MQ=255)
                     ggAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATa                                                                                                                               <  7:303987/149‑1 (MQ=255)
                      gagaTTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaa                                                                                                                              <  8:61338/149‑1 (MQ=255)
                        gaTTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaaca                                                                                                                            <  2:275214/149‑1 (MQ=255)
                         aTTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGAGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaacaa                                                                                                                           >  8:120290/1‑149 (MQ=255)
                           tGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaacaacc                                                                                                                         >  5:173794/1‑148 (MQ=255)
                              tttttATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACAccc                                                                                                                      >  7:411138/1‑149 (MQ=255)
                               ttttATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTc                                                                                                                                           >  8:94340/1‑127 (MQ=255)
                               ttttATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGAt                                                                                                                                <  6:99384/138‑1 (MQ=255)
                               ttttATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGAt                                                                                                                                >  5:99384/1‑138 (MQ=255)
                                    tCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCgg                                                                                                                >  8:262607/1‑149 (MQ=255)
                                         tGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAAtg                                                                                                           <  1:68466/149‑1 (MQ=255)
                                                       tATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaacaac                                                                                                                          <  6:277921/120‑1 (MQ=255)
                                                       tATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATaacaac                                                                                                                          >  5:277921/1‑120 (MQ=255)
                                                       tATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTg                                                                                              >  6:124663/1‑148 (MQ=255)
                                                        aTGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCACAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTg                                                                                              >  5:209392/1‑147 (MQ=255)
                                                           aaCTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAgcg                                                                                         >  5:316274/1‑149 (MQ=255)
                                                             cTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCa                                                                                       <  8:411138/149‑1 (MQ=255)
                                                               ggATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCt                                                                                                                                                 <  6:246805/89‑1 (MQ=255)
                                                               ggATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCt                                                                                                                                                 >  5:246805/1‑89 (MQ=255)
                                                                gATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGtt                                                                                    >  5:242236/1‑149 (MQ=255)
                                                                    ccTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCt                                                                                <  1:425050/149‑1 (MQ=255)
                                                                     cTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCACTGCTGCAAACCGGTGAAGCGCAGTTTGCtt                                                                               >  5:101502/1‑149 (MQ=255)
                                                                      tGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCttt                                                                              >  7:208862/1‑149 (MQ=255)
                                                                       ggAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGAATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTc                                                                             >  4:390110/1‑149 (MQ=255)
                                                                              gACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTc                                                                      <  3:4040/149‑1 (MQ=255)
                                                                                 cGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCt                                                                    >  8:110724/1‑148 (MQ=255)
                                                                                   aTTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTAc                                                                 <  7:262607/149‑1 (MQ=255)
                                                                                   aTTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCACAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTAc                                                                 <  6:209392/149‑1 (MQ=255)
                                                                                        gtgAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAAc                                                                                                   >  5:70331/1‑110 (MQ=255)
                                                                                        gtgAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAAc                                                                                                   <  6:70331/110‑1 (MQ=255)
                                                                                         tgAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCa                                                            <  6:316274/148‑1 (MQ=255)
                                                                                         tgAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGGGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAg                                                           <  5:144104/149‑1 (MQ=255)
                                                                                          gAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTccc                                                                           <  8:166876/132‑1 (MQ=255)
                                                                                          gAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTccc                                                                           >  7:166876/1‑132 (MQ=255)
                                                                                               tGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCcaca                                                     <  6:371726/149‑1 (MQ=255)
                                                                                                  aaaaaTTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCa                                                            >  4:234861/1‑139 (MQ=255)
                                                                                                  aaaaaTTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCa                                                            <  3:234861/139‑1 (MQ=255)
                                                                                                   aaaaTTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACACCACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACActgc                                                 >  5:125207/1‑149 (MQ=255)
                                                                                                       ttCGCGGGTTACTGGCAGCCAGAATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGga                                             <  3:390110/149‑1 (MQ=255)
                                                                                                         cgcgGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGgag                                            <  7:120290/148‑1 (MQ=255)
                                                                                                             gggtACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCATTGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACAGTGCTGGAGaaaaa                                       <  6:173794/146‑1 (MQ=255)
                                                                                                             ggTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGaaaaa                                       <  5:124663/149‑1 (MQ=255)
                                                                                                                  cTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACaaa                                   <  8:208862/148‑1 (MQ=255)
                                                                                                                   tGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTgg                                              >  6:250467/1‑136 (MQ=255)
                                                                                                                   tGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTgg                                              <  5:250467/136‑1 (MQ=255)
                                                                                                                   tGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACaaaaa                                 >  1:323503/1‑149 (MQ=255)
                                                                                                                     gCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAata                               >  7:35117/1‑149 (MQ=255)
                                                                                                                         ccAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGa                           <  5:360711/149‑1 (MQ=255)
                                                                                                                          cAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGa                                                               >  5:334566/1‑112 (MQ=255)
                                                                                                                          cAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGa                                                               <  6:334566/112‑1 (MQ=255)
                                                                                                                             gATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGtt                        >  6:279563/1‑148 (MQ=255)
                                                                                                                              aTTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGa                      >  2:314493/1‑149 (MQ=255)
                                                                                                                              aTTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGa                      <  6:242236/149‑1 (MQ=255)
                                                                                                                               ttGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGAt                     <  1:314493/149‑1 (MQ=255)
                                                                                                                                  cccAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  2:4031/1‑149 (MQ=255)
                                                                                                                                  cccAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  6:278831/1‑149 (MQ=255)
                                                                                                                                  cccAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  3:87115/1‑149 (MQ=255)
                                                                                                                                  cccAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATCCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  4:110254/1‑149 (MQ=255)
                                                                                                                                   ccAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  6:94761/1‑148 (MQ=255)
                                                                                                                                   ccAAACTGGACAGCACAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATg                    <  7:103700/146‑1 (MQ=255)
                                                                                                                                   ccAAACTGGACAGCACAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATg                    >  8:103700/1‑146 (MQ=255)
                                                                                                                                    cAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGc                  >  2:119884/1‑147 (MQ=255)
                                                                                                                                    cAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCa                <  4:87115/149‑1 (MQ=255)
                                                                                                                                             cAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCCGTCaa       <  3:128444/149‑1 (MQ=255)
                                                                                                                                                cATAACCTGGCGTCCGGTGGCGGATAACAAAACCCGCGCGGCAATGCTGCAAACCGGTGAATCGCAGTTTGCTTTCACCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCAGTCAATTa    >  6:389336/1‑149 (MQ=255)
                                                                                                                                                  tAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCCGTCAATTATg  >  2:38936/1‑149 (MQ=255)
                                                                                                                                                  tAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCCGTCAATTATg  >  4:173068/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAACTGGACAGCATAACCTGGCGTCCGGTGGCGGATAACAACACCCGCGCGGCAATGCTGCAAACCGGTGAAGCGCAGTTTGCTTTCCCCATTCCTTACGAGCAGGCCACACTGCTGGAGAAAAACAAAAATATCGAGTTGATGGCCAGTCCGTCAATTATG  >  NC_000913/869917‑870211

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