Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,248,023 G→A W84* (TGG→TAG)  yqjA → general envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,248,0230GA98.8% 295.7 / NA 85W84* (TGG→TAG) yqjAgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/47);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTGGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAACGCGCACACCA  >  NC_000913/3247926‑3248119
                                                                                                 |                                                                                                
aGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCACGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTagg                                                                                                <  1:2109442/100‑1 (MQ=255)
   ttACTGGTATTGGTCGGCGTGTTGATTGCGAATGGCGCGCTGCGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCGGCCTCGGCTGCTAGGTCa                                                                                             <  1:2104077/100‑1 (MQ=255)
   ttACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCa                                                                                             <  2:2403509/100‑1 (MQ=255)
     aCTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGTTGCTAGGTCAGc                                                                                           <  2:1971582/100‑1 (MQ=255)
     aCTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGc                                                                                           <  2:1058099/100‑1 (MQ=255)
     aCTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGc                                                                                           <  2:2781424/100‑1 (MQ=255)
      cTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGc                                                                                           >  1:2002681/1‑99 (MQ=255)
      cTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCt                                                                                          >  2:2752422/1‑100 (MQ=255)
      cTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCt                                                                                          <  1:1374051/100‑1 (MQ=255)
       tGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCta                                                                                         >  1:2577158/1‑100 (MQ=255)
          tATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTatat                                                                                      <  1:2564987/100‑1 (MQ=255)
          tATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTatat                                                                                      <  1:2615449/100‑1 (MQ=255)
              ggTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAg                                                                                  <  2:2632640/100‑1 (MQ=255)
               gTCGGCGAGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAgg                                                                                 <  1:2311626/100‑1 (MQ=255)
                 cGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAgggg                                                                               >  1:646043/1‑100 (MQ=255)
                 cGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAggg                                                                                >  1:1446937/1‑99 (MQ=255)
                   gCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCg                                                                             >  1:2491265/1‑100 (MQ=255)
                   gCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCg                                                                             <  1:1339533/100‑1 (MQ=255)
                      tgtTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATg                                                                          <  2:688474/100‑1 (MQ=255)
                       gtTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGAATCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATgg                                                                         >  1:1845187/1‑100 (MQ=255)
                         tGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCt                                                                       >  1:1913194/1‑100 (MQ=255)
                         tGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCt                                                                       >  2:2983737/1‑100 (MQ=255)
                          gATTGCGAAAGGCGCGATGGGCTATCTGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTg                                                                      >  2:2896082/1‑100 (MQ=255)
                            ttGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTggg                                                                    >  1:2048928/1‑100 (MQ=255)
                             tGCGAAAGGCGCGATGGGCTGTCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCACGCGCCATGGcccg                                                                     >  1:1833344/1‑95 (MQ=255)
                             tGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGc                                                                   <  2:2348852/100‑1 (MQ=255)
                             tGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGc                                                                   <  2:746738/100‑1 (MQ=255)
                             tGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGc                                                                   <  2:2513486/100‑1 (MQ=255)
                             tGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGc                                                                   <  2:1064334/100‑1 (MQ=255)
                              gCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCa                                                                  <  1:859989/100‑1 (MQ=255)
                              gCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCa                                                                  >  2:899362/1‑100 (MQ=255)
                              gCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCa                                                                  >  2:650472/1‑100 (MQ=255)
                                gAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAAt                                                                >  1:1978535/1‑100 (MQ=255)
                                gAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAAt                                                                <  1:1585932/100‑1 (MQ=255)
                                 aaaGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATa                                                               >  1:2190183/1‑100 (MQ=255)
                                    ggCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCCATAcac                                                            >  2:2883739/1‑98 (MQ=255)
                                    ggCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATAccc                                                            <  1:242834/100‑1 (MQ=255)
                                      cgcgATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGc                                                          <  2:1519832/100‑1 (MQ=255)
                                        cgATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCAc                                                        >  2:855874/1‑100 (MQ=255)
                                         gATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCAcc                                                       <  1:1181403/100‑1 (MQ=255)
                                          aTGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGTGCGATGGCTGGGCAATACCCGCACCg                                                      <  1:8705/100‑1 (MQ=255)
                                          aTGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCg                                                      <  2:2485541/100‑1 (MQ=255)
                                          aTGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCg                                                      <  2:347858/100‑1 (MQ=255)
                                          aTGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCg                                                      <  2:943335/100‑1 (MQ=255)
                                           tGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGt                                                     <  2:2769594/100‑1 (MQ=255)
                                            gggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTa                                                    >  2:2312099/1‑100 (MQ=255)
                                             ggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTAc                                                   >  1:1259272/1‑100 (MQ=255)
                                             ggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTAc                                                   >  2:1206975/1‑100 (MQ=255)
                                             ggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTAc                                                   >  2:1241176/1‑100 (MQ=255)
                                             ggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTAc                                                   >  2:2490324/1‑100 (MQ=255)
                                             ggCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTAc                                                   >  1:1761986/1‑100 (MQ=255)
                                              gCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACa                                                  <  2:2975127/100‑1 (MQ=255)
                                               cTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCATGTGCATTGGCTGGGCAATACCCGCACCg                                                      <  1:1399333/95‑1 (MQ=255)
                                                tATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACaaa                                                <  2:2852693/100‑1 (MQ=255)
                                                   ccGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACt                                             >  1:524317/1‑100 (MQ=255)
                                                    cGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTg                                            >  1:1660254/1‑100 (MQ=255)
                                                    cGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTg                                            >  1:2783545/1‑100 (MQ=255)
                                                    cGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTg                                            <  2:1340161/100‑1 (MQ=255)
                                                    cGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTg                                            >  1:2769622/1‑100 (MQ=255)
                                                     gCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTgg                                           >  2:919430/1‑100 (MQ=255)
                                                     gCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTgg                                           <  2:361443/100‑1 (MQ=255)
                                                            aTTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTAtctc                                    >  1:1782364/1‑100 (MQ=255)
                                                              tCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCa                                   <  1:1098250/99‑1 (MQ=255)
                                                               ctgctgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCAtt                                 <  2:849281/100‑1 (MQ=255)
                                                               ctgctgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCAtt                                 >  2:2713251/1‑100 (MQ=255)
                                                               ctgctgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCAtt                                 <  1:1895193/100‑1 (MQ=255)
                                                                tgctgctgAACGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCAttt                                <  2:2393007/100‑1 (MQ=255)
                                                                 gctgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTa                               >  1:2151390/1‑100 (MQ=255)
                                                                 gctgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTa                               <  1:338447/100‑1 (MQ=255)
                                                                   tgctgACCGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTAcc                             <  1:150730/100‑1 (MQ=255)
                                                                       gACCGTTGCCGCCAGCCTCGGCTGCTGGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCcgcg                         >  2:1827761/1‑100 (MQ=255)
                                                                          cGTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCAt                      <  2:646045/100‑1 (MQ=255)
                                                                           gTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCAtt                     >  2:1993339/1‑100 (MQ=255)
                                                                           gTTGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCAtt                     >  1:1623545/1‑100 (MQ=255)
                                                                           gTTGCCGCCAGCCTCGGCTGATAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCAtt                     >  1:1829580/1‑100 (MQ=255)
                                                                            ttGCCGCCAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTa                    <  2:2290275/100‑1 (MQ=255)
                                                                               ccgccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatca                 <  1:2143514/100‑1 (MQ=255)
                                                                               ccgccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatca                 <  2:150746/100‑1 (MQ=255)
                                                                               ccgccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGAAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatca                 <  1:1016504/100‑1 (MQ=255)
                                                                                cgccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatcat                <  1:1428055/100‑1 (MQ=255)
                                                                                cgccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatcat                >  2:453516/1‑100 (MQ=255)
                                                                                 gccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatcatc               <  2:1845195/100‑1 (MQ=255)
                                                                                 gccAGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTatcatc               <  1:855871/100‑1 (MQ=255)
                                                                                    aGCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAAc            <  2:2577170/100‑1 (MQ=255)
                                                                                     gCCTCGGCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAAcg           <  1:650470/100‑1 (MQ=255)
                                                                                          ggCTGCTAGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAACGCGcac      >  2:565273/1‑100 (MQ=255)
                                                                                             tgctaGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAACGCGCACAcc   >  1:2517180/1‑100 (MQ=255)
                                                                                              gctaGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAACGCGCACACCa  >  2:682767/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGTTTACTGGTATTGGTCGGCGTGTTGATTGCGAAAGGCGCGATGGGCTATCCGCAAACGATTCTGCTGCTGACCGTTGCCGCCAGCCTCGGCTGCTGGGTCAGCTATATTCAGGGGCGATGGCTGGGCAATACCCGCACCGTACAAAACTGGCTATCTCATTTACCCGCGCATTATCATCAACGCGCACACCA  >  NC_000913/3247926‑3248119

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