Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 868,748 Δ1 bp 100% coding (1229/1872 nt) gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,7460C.56.0% 25.4 / 308.5 141coding (1227/1872 nt)gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base . (40/39);  total (73/68)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

Only 100 of 148 total aligned reads displayed.

GCAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCA  >  NC_000913/868647‑868842
                                                                                                   |                                                                                                
gCAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAc                                                                                                  >  1:2074885/1‑100 (MQ=255)
gCAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAc                                                                                                  >  2:756685/1‑100 (MQ=255)
 cAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAcc                                                                                                 <  1:2733269/100‑1 (MQ=255)
 cAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAATTTATTTTTCAGGAc                                                                                                  >  2:2419743/1‑99 (MQ=255)
  aGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAcc                                                                                                 >  2:3732170/1‑99 (MQ=255)
     gcggcgAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGAccttac                                                                                             >  1:3789645/1‑96 (MQ=255)
         cgAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACCCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCtt                                                                                        >  2:2583871/1‑100 (MQ=255)
            aaTTATCTTTAACGGTCAGCGAATCGATACCGTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGc                                                                                     <  2:4096772/100‑1 (MQ=255)
             aTTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGc                                                                                     >  1:2529698/1‑100 (MQ=255)
              ttATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTg                                                                                   <  2:1898478/100‑1 (MQ=255)
               tATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTg                                                                                   >  1:1672320/1‑100 (MQ=255)
               tATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTg                                                                                   <  1:2917467/100‑1 (MQ=255)
               tATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTgg                                                                                  >  1:2653420/1‑100 (MQ=255)
                aTCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGa                                                                                 <  1:2877883/100‑1 (MQ=255)
                 tCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCCGGATATTCAGTTTATTTTTCAGGA‑CCTTGCGCTTCGCTgggc                                                                                <  1:2797165/100‑3 (MQ=255)
                  cTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAACCTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGa                                                                                 <  2:674396/99‑1 (MQ=255)
                  cTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGAcc                                                                               <  2:4118333/100‑1 (MQ=255)
                   tttAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGAccc                                                                              >  2:3574187/1‑100 (MQ=255)
                   tttAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGAccc                                                                              >  1:2289902/1‑100 (MQ=255)
                    ttAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGAccc                                                                              >  2:3840510/1‑100 (MQ=255)
                    ttAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCa                                                                             <  1:2732905/100‑1 (MQ=255)
                     tAACTGTCAGCAAATCGATAACTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGAcc                                                                               <  2:2799868/98‑1 (MQ=255)
                     tAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCAc                                                                            >  1:1090395/1‑100 (MQ=255)
                      aaCGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCAc                                                                            >  2:670492/1‑100 (MQ=255)
                       aCGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACg                                                                           >  1:582719/1‑100 (MQ=255)
                        cGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGt                                                                          >  1:2135023/1‑100 (MQ=255)
                        cGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTc                                                                         <  1:25478/100‑1 (MQ=255)
                        cGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTc                                                                         >  1:1491275/1‑100 (MQ=255)
                         ggTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTc                                                                         <  2:2958494/100‑1 (MQ=255)
                         ggTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTc                                                                         <  1:3558970/100‑1 (MQ=255)
                          gTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCan                                                                       <  2:3260355/100‑2 (MQ=255)
                          gTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAg                                                                       <  1:4092834/100‑1 (MQ=255)
                            cAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGAc                                                                     >  1:2191296/1‑100 (MQ=255)
                              gCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGAcc                                                                    >  2:1698504/1‑100 (MQ=255)
                              gCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCa                                                                   >  1:2216573/1‑100 (MQ=255)
                               cGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCAt                                                                  >  2:4010071/1‑100 (MQ=255)
                                 aaTCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATc                                                                 <  1:2852181/100‑1 (MQ=255)
                                 aaTCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTTGACCCACGTCAGACCATCg                                                                <  2:1275163/100‑1 (MQ=255)
                                   ttgATACCTTGTCACCGGGCAGCCTCCACGCATTACGCCGTGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGTCCAtgtg                                                               <  2:3485310/97‑4 (MQ=255)
                                    cGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGTACCCACGTCAGACCATCGGTg                                                             >  2:318516/1‑100 (MQ=255)
                                    cGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTg                                                             <  2:590478/100‑1 (MQ=255)
                                     gATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTg                                                             >  2:240105/1‑100 (MQ=255)
                                     gATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTg                                                             >  1:337377/1‑99 (MQ=255)
                                       tACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGAt                                                           >  1:1752520/1‑100 (MQ=255)
                                        aacTTGTAACCCGGCAAACTTCAGGCATAACGCCGCGATATCCAGCCTATTCTTCACG‑CCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGAttg                                                         <  2:341444/98‑2 (MQ=255)
                                         ccTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCg                                                        >  1:3943489/1‑100 (MQ=255)
                                           ttGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGCTTCGa                                                       >  2:2201098/1‑100 (MQ=255)
                                            tatCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTACAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGAtt                                                     <  2:1149940/98‑1 (MQ=255)
                                             gTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGAtt                                                     >  2:2051265/1‑100 (MQ=255)
                                              tCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTa                                                    <  1:2343712/100‑1 (MQ=255)
                                               cACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                  <  1:3144949/100‑1 (MQ=255)
                                                aCCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                  <  2:1997376/100‑1 (MQ=255)
                                                aCCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCg                                                 >  1:3563679/1‑100 (MQ=255)
                                                aCCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCg                                                 <  2:183816/100‑1 (MQ=255)
                                                 cccGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                  >  1:1845664/1‑99 (MQ=255)
                                                 cccGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGa                                                >  1:730834/1‑100 (MQ=255)
                                                 accGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATc                                                  >  2:1012150/2‑98 (MQ=255)
                                                  ccGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGaa                                               >  2:1817615/1‑100 (MQ=255)
                                                    ggCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAAc                                              <  1:2603846/100‑1 (MQ=255)
                                                      cAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCg                                            >  1:3972911/1‑100 (MQ=255)
                                                         aCTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTg                                         <  2:2937615/100‑1 (MQ=255)
                                                         aCTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGc                                        >  1:463959/1‑100 (MQ=255)
                                                          cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTg                                         >  2:2450196/1‑99 (MQ=255)
                                                          cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGc                                        >  2:569178/1‑100 (MQ=255)
                                                          cTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCg                                       <  1:1817628/100‑1 (MQ=255)
                                                          cTTCAGGCATTACGCCGGGATATTCAGTCTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCACACCATCGGTGATTCCATTATCCAACCGCtcc                                        >  2:3309893/1‑98 (MQ=255)
                                                             ccggCATTACGCCGTGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGCCCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCgtg                                     <  2:3018248/98‑1 (MQ=255)
                                                             cAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCgtg                                     >  2:2357798/1‑100 (MQ=255)
                                                               ggCATTACGCCGGGATATTCAGTTTATTTTTCAGTA‑CCTTACGCTTCGCTGGGCCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                  <  2:1920329/100‑1 (MQ=255)
                                                               ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTa                                   >  2:3654481/1‑100 (MQ=255)
                                                               ggCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                  >  1:1634416/1‑100 (MQ=255)
                                                                gCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                  <  2:2453144/100‑1 (MQ=255)
                                                                gCATTACGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTac                                  <  1:73954/99‑1 (MQ=255)
                                                                 cATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTaca                                 <  2:295473/100‑1 (MQ=255)
                                                                     aCGCCGGGATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                            <  2:1485590/100‑1 (MQ=255)
                                                                      cGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGtt                            >  2:1338472/1‑100 (MQ=255)
                                                                          gggATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTAtt                        <  2:1760529/100‑1 (MQ=255)
                                                                          gggATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTAtt                        <  2:3641016/99‑1 (MQ=255)
                                                                          gggATATTCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTg                       <  2:145443/100‑1 (MQ=255)
                                                                              tatTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCa                    >  2:1595360/1‑100 (MQ=255)
                                                                               atTCAGTTTATTTTTCAGGA‑ACTTACGCTTCGCTGGACCCACGTCAGATCATCGGTGATTCGATTATCGAACCGCTTCGTGTACACAGTTTATTGCCAgg                  >  2:2509140/1‑100 (MQ=255)
                                                                                ttCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAgg                  >  1:1345024/1‑100 (MQ=255)
                                                                                ttCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGt                 >  2:925813/1‑100 (MQ=255)
                                                                                ttCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGt                 <  1:1989344/100‑1 (MQ=255)
                                                                                 tCAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTa                <  2:3544620/100‑1 (MQ=255)
                                                                                  cAGTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCGACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTa                >  2:2654239/1‑99 (MQ=255)
                                                                                   aGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTaa               >  1:2328212/1‑100 (MQ=255)
                                                                                    gTTTATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAg             >  1:2961217/1‑100 (MQ=255)
                                                                                     tttATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGa            <  2:1530408/100‑1 (MQ=255)
                                                                                      ttATTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAc           >  1:3651227/1‑100 (MQ=255)
                                                                                        aTTTTTCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgc         >  2:3987923/1‑100 (MQ=255)
                                                                                         tttttCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgcg        <  2:2517879/100‑1 (MQ=255)
                                                                                         tttttCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgcg        >  1:3409838/1‑100 (MQ=255)
                                                                                          ttttCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGAcgcg        >  2:1457178/1‑100 (MQ=255)
                                                                                           tttCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCgg       >  1:2884815/1‑100 (MQ=255)
                                                                                            ttCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGc      >  2:3186808/1‑100 (MQ=255)
                                                                                            ttCAGGA‑CCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGCCGCGGCt     <  1:2175846/100‑1 (MQ=255)
                                                                                              cAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTg    >  2:1479772/1‑100 (MQ=255)
                                                                                                ggACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCa  >  2:1493029/1‑100 (MQ=255)
                                                                                                ggACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCa  <  1:625416/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GCAGGGCGGCGAAATTATCTTTAACGGTCAGCGAATCGATACCTTGTCACCCGGCAAACTTCAGGCATTACGCCGGGATATTCAGTTTATTTTTCAGGACCCTTACGCTTCGCTGGACCCACGTCAGACCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCA  >  NC_000913/868647‑868842

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