Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 872,398 C→T 5.5% A170V (GCA→GTA)  gsiD → glutathione ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913872,3980CT 5.5% 239.3 / 2.3 73A170V (GCA→GTA) gsiDglutathione ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (38/31);  new base T (2/2);  total (40/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 = 4.96e-01

TGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTC  >  NC_000913/872251‑872546
                                                                                                                                                   |                                                                                                                                                    
tGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCa                                                                                                                                                     <  4:229559/149‑1 (MQ=255)
      gCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGc                                                                                                                                                >  1:80735/1‑148 (MQ=255)
        ttGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCAt                                                                                                                                             >  3:296554/1‑149 (MQ=255)
            tCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTTCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTATCGTGATTATTGCAGTCGCCAttttt                                                                                                                                         >  5:287095/1‑149 (MQ=255)
            tCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCAttttt                                                                                                                                         >  5:123031/1‑149 (MQ=255)
              gCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCAtttttt                                                                                                                                        <  2:80735/148‑1 (MQ=255)
                tGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTc                                                                                                                                       <  5:241106/147‑1 (MQ=255)
                 ggATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCAt                                                                                                                                    >  1:55169/1‑149 (MQ=255)
                    tattatGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGGCATTTTTTCCATccc                                                                                                                                 >  5:196432/1‑149 (MQ=255)
                      ttatGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAAAGTGATTATTGCAGTCGCCATTTTTTCCATCCCcg                                                                                                                               >  5:294337/1‑149 (MQ=255)
                        atGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGTGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCcgc                                                                                                                              <  2:413863/148‑1 (MQ=255)
                                               aTCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCaa                                                                                                      >  8:300079/1‑149 (MQ=255)
                                               aTCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCaa                                                                                                      >  7:211003/1‑149 (MQ=255)
                                                  aTGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAacac                                                                                                   >  1:30499/1‑149 (MQ=255)
                                                     cgcATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGc                                                                                                 <  7:97037/148‑1 (MQ=255)
                                                      gcATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTg                                                                                               <  4:327929/149‑1 (MQ=255)
                                                      gcATTTGCGATGTGCTGTTTGCCTTCACGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGAATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTg                                                                                               >  7:413342/1‑149 (MQ=255)
                                                         tttGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGgtg                                                                                            <  4:296554/149‑1 (MQ=255)
                                                         tttGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGgtg                                                                                            >  1:305315/1‑149 (MQ=255)
                                                             cGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGa                                                                                        <  3:253483/149‑1 (MQ=255)
                                                               aTGTGCTGGTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGATGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGGCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGaaa                                                                                      <  5:75635/149‑1 (MQ=255)
                                                                    cTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGGTGTTGCGGTGGTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTTCGCGGCAACACGCTGGTGTTGAAACAGc                                                                                  <  6:287095/148‑1 (MQ=255)
                                                                     tGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTATCAATCCCCGCGATTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACATCa                                                                                 >  1:338623/1‑148 (MQ=255)
                                                                       tttGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGTAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                              >  8:283515/1‑149 (MQ=255)
                                                                       tttGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGTTGCGCGGCAACACGCTGGTGTTGAAACAGCaaa                                                                               <  3:157418/148‑1 (MQ=255)
                                                                       tttGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                              >  4:11974/1‑149 (MQ=255)
                                                                       tttGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                              >  2:251476/1‑149 (MQ=255)
                                                                       tttGCCTTCCCGGGTATTTTAATGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGTAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                              >  8:290457/1‑149 (MQ=255)
                                                                        ttGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAAc                                                                              <  7:300079/148‑1 (MQ=255)
                                                                            cTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTa                                                                         <  2:30499/149‑1 (MQ=255)
                                                                            cTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCAGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTa                                                                         >  2:248476/1‑149 (MQ=255)
                                                                             ttCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAt                                                                        >  2:97230/1‑149 (MQ=255)
                                                                               cccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAtt                                                                       <  8:98454/148‑1 (MQ=255)
                                                                               cccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAtt                                                                       <  3:301982/148‑1 (MQ=255)
                                                                                ccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGttt                                                                                                                          >  1:277778/1‑96 (MQ=255)
                                                                                ccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGttt                                                                                                                          <  2:277778/96‑1 (MQ=255)
                                                                                ccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAtt                                                                       >  5:26004/1‑147 (MQ=255)
                                                                                ccGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGa                                                                     <  1:251476/149‑1 (MQ=255)
                                                                                 cGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAg                                                                    >  1:291904/1‑149 (MQ=255)
                                                                                 cGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAg                                                                    >  2:225840/1‑149 (MQ=255)
                                                                                  gggTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGTAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGt                                                                   <  7:290457/149‑1 (MQ=255)
                                                                                  gggTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGTAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGt                                                                   <  7:283515/149‑1 (MQ=255)
                                                                                  gggTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGa                                                                     <  8:15398/147‑1 (MQ=255)
                                                                                  gggTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGa                                                                     >  7:15398/1‑147 (MQ=255)
                                                                                                cGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTAtt                                                     <  6:196432/149‑1 (MQ=255)
                                                                                                  aTCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTgg                                                   <  1:20747/149‑1 (MQ=255)
                                                                                                     gCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGc                                                >  3:68969/1‑149 (MQ=255)
                                                                                                      cTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGc                                                >  7:154918/1‑148 (MQ=255)
                                                                                                       tgttgtTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCt                                                                            >  6:348040/1‑119 (MQ=255)
                                                                                                       tgttgtTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCt                                                                            <  5:348040/119‑1 (MQ=255)
                                                                                                       tgttgtTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCAGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTGTATTGAGTCAGCACGCAGTATTGGTGc                                                >  2:77720/1‑147 (MQ=255)
                                                                                                         ttgttgCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGc                                            <  8:211003/149‑1 (MQ=255)
                                                                                                              gCGGAGTTGGGAAGCGGCATTTCTAACGTGATTATTGCAGTCGCCATTATTTCCATCCCCGCGTTTGCCCGCCTGCTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGatat                                       >  2:81365/1‑149 (MQ=255)
                                                                                                               cGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATg                                      >  8:263275/1‑149 (MQ=255)
                                                                                                                    ttGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATCGGTGCCAGCGATATGACCGt                                 >  8:188946/1‑149 (MQ=255)
                                                                                                                      gggAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAt                                                                        >  5:243660/1‑108 (MQ=255)
                                                                                                                      gggAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTAt                                                                        <  6:243660/108‑1 (MQ=255)
                                                                                                                       ggAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGtttt                              <  4:68969/149‑1 (MQ=255)
                                                                                                                        gAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttg                             <  1:248476/149‑1 (MQ=255)
                                                                                                                        gAAGCGGCATTGCTAACGTGATTAATGCAGTCGCCATTTTTTACATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGACACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttg                             >  5:338134/1‑149 (MQ=255)
                                                                                                                         aaGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttgt                            <  7:263275/149‑1 (MQ=255)
                                                                                                                         aaGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttgt                            >  3:360531/1‑149 (MQ=255)
                                                                                                                          aGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTGCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttgtt                           <  6:294337/149‑1 (MQ=255)
                                                                                                                          aGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTttgtt                           <  8:331775/149‑1 (MQ=255)
                                                                                                                                 ttGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGa                                          <  4:169355/127‑1 (MQ=255)
                                                                                                                                 ttGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGa                                          >  3:169355/1‑127 (MQ=255)
                                                                                                                                     tAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGatat                                       <  8:431610/126‑1 (MQ=255)
                                                                                                                                     tAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGatat                                       >  7:431610/1‑126 (MQ=255)
                                                                                                                                     tAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATcc                >  3:29839/1‑149 (MQ=255)
                                                                                                                                           gATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATcctgcctg          >  2:28087/1‑149 (MQ=255)
                                                                                                                                              tattGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGa       >  4:174734/1‑149 (MQ=255)
                                                                                                                                               attGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGAc      >  4:135562/1‑149 (MQ=255)
                                                                                                                                                ttGCAGTCGCCATGTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGGATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGAcc     <  5:276670/149‑1 (MQ=255)
                                                                                                                                                   cAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGtc  <  3:11974/149‑1 (MQ=255)
                                                                                                                                                   cAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGGGTTGAAACAGCAAACCTTTATTGAGGCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGt   <  1:28087/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTCAGCACGCAGTATTGGTGCCAGCGATATGACCGTTTTGTTGCGTCATATCCTGCCTGGGACCGTC  >  NC_000913/872251‑872546

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: