Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 868,920 T→A 5.4% G467G (GGT→GGA gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,9200TA 5.4% 220.0 / 2.8 74G467G (GGT→GGAgsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base A (2/2);  total (36/38)
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.01e-01

CCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCG  >  NC_000913/868775‑869067
                                                                                                                                                 |                                                                                                                                                   
ccATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGTTAAAGACGCGGCTGCACGCGTTGCGTGGGTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCTTTACCCGCATGAGTTTTCCGGCTGTca                                                                                                                                                   <  3:321601/148‑1 (MQ=255)
       tGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgccag                                                                                                                                            <  7:324754/148‑1 (MQ=255)
          ttCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgccagcg                                                                                                                                          >  2:61038/1‑147 (MQ=255)
                 aTCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGc                                                                                                                                   <  2:356037/147‑1 (MQ=255)
                  tCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCa                                                                                                                                  <  3:272692/147‑1 (MQ=255)
                     aaCCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGc                                                                                                                              >  2:289990/1‑148 (MQ=255)
                             cGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGTTTGCTGGAGCGCGTGGGCCTGTTACCTGACCATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGt                                                                                                                      >  6:79480/1‑148 (MQ=255)
                               tgtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTg                                                                                                                    <  7:415842/148‑1 (MQ=255)
                                gtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgcc                                                                                                                                              >  6:335721/1‑121 (MQ=255)
                                gtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTcagcgcc                                                                                                                                              <  5:335721/121‑1 (MQ=255)
                                gtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGtt                                                                                                                     <  6:87501/146‑1 (MQ=255)
                                gtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGtt                                                                                                                     >  5:87501/1‑146 (MQ=255)
                                gtACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTgg                                                                                                                   >  6:83873/1‑148 (MQ=255)
                                 tACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGc                                                                                                                  <  4:160110/148‑1 (MQ=255)
                                     cGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATc                                                                                                                                      <  5:66583/124‑1 (MQ=255)
                                     cGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATc                                                                                                                                      >  6:66583/1‑124 (MQ=255)
                                                cAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGt                                                                                                   >  5:241471/1‑148 (MQ=255)
                                                cAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGt                                                                                                   >  7:18940/1‑148 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTTCCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGAt                                                                                                >  3:3390/1‑147 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGAt                                                                                                >  1:43285/1‑147 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGAt                                                                                                >  1:38845/1‑147 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATc                                                                                               <  5:83873/148‑1 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATc                                                                                               <  3:132922/148‑1 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATc                                                                                               <  8:139855/148‑1 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGACAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCaaa                                                                                                      >  7:348698/1‑141 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGACAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCaaa                                                                                                      <  8:348702/141‑1 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGACAGCGCCAGCGCATCTGCATTGCTCGCGCGGTGGCATTGAATCCaaa                                                                                                      <  8:348698/141‑1 (MQ=255)
                                                    tAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCAGCATGAGTTTTCCGGCGGACAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCaaa                                                                                                      >  7:348702/1‑141 (MQ=255)
                                                         aCGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGc                                                                                          >  7:242908/1‑148 (MQ=255)
                                                                tGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCcgacg                                                                                     >  1:299477/1‑146 (MQ=255)
                                                                         tGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGc                                                                                                                                   <  3:283177/91‑1 (MQ=255)
                                                                         tGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGc                                                                                                                                   >  4:283177/1‑91 (MQ=255)
                                                                           cGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCgg                                                                         <  7:176959/147‑1 (MQ=255)
                                                                           cGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGgc                                                                        >  4:122685/1‑148 (MQ=255)
                                                                                   cTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGAt                                                                 <  2:299477/147‑1 (MQ=255)
                                                                                   cTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATg                                                                <  8:242908/148‑1 (MQ=255)
                                                                                   cTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATg                                                                <  3:70412/148‑1 (MQ=255)
                                                                                      gAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGtt                                                              <  2:38845/147‑1 (MQ=255)
                                                                                            gTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTc                                                       <  1:61038/148‑1 (MQ=255)
                                                                                                    gtaaCCTGAACATGCCAGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCCTCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGACGTTTCGGCGCTGGATGTTTCTATTCGCggg                                                  >  6:135492/4‑145 (MQ=255)
                                                                                                    gTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGtt                                                              <  2:39423/133‑1 (MQ=255)
                                                                                                    gTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGtt                                                              >  1:39423/1‑133 (MQ=255)
                                                                                                        ccTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTc                                                            >  1:287488/1‑131 (MQ=255)
                                                                                                        ccTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTc                                                            <  2:287488/131‑1 (MQ=255)
                                                                                                         cTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAg                                               <  2:72279/143‑1 (MQ=255)
                                                                                                         cTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAg                                               >  1:72279/1‑143 (MQ=255)
                                                                                                             aCATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAAc                                      >  1:7126/1‑148 (MQ=255)
                                                                                                               aTGCCTGGCTTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGGGCATCTGCATTGCTCGCGCGTTGGCATTGATTCCAAAAGTTATCATTGCCTACGAAGCCGTTTCTGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACtt                                    <  1:181736/148‑1 (MQ=255)
                                                                                                               aTGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACtt                                    >  4:373049/1‑148 (MQ=255)
                                                                                                                 gCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACttgt                                  >  4:23702/1‑148 (MQ=255)
                                                                                                                    tGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTg                                                                                                  <  4:220172/81‑1 (MQ=255)
                                                                                                                    tGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTg                                                                                                  >  3:220172/1‑81 (MQ=255)
                                                                                                                     ggCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATc                                                                                               >  7:377681/1‑83 (MQ=255)
                                                                                                                     ggCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATc                                                                                               <  8:377681/83‑1 (MQ=255)
                                                                                                                          tACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTgg                                                                   >  3:271622/1‑106 (MQ=255)
                                                                                                                          tACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTgg                                                                   <  4:271622/106‑1 (MQ=255)
                                                                                                                            cccGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGAtctc                       <  6:241471/148‑1 (MQ=255)
                                                                                                                               gCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAg                    <  2:43285/148‑1 (MQ=255)
                                                                                                                                   gAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCg                  >  5:322653/1‑146 (MQ=255)
                                                                                                                                   gAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGt                 >  3:271406/1‑147 (MQ=255)
                                                                                                                                   gAGCTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCg                  <  6:322653/146‑1 (MQ=255)
                                                                                                                                    aGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTg                <  1:289990/147‑1 (MQ=255)
                                                                                                                                     gTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCcg                                                                                        >  7:108020/1‑74 (MQ=255)
                                                                                                                                     gTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCcg                                                                                        <  8:108020/74‑1 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCCGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCTGCGCTGGATTTATCTATTCGAGGGTAGATTATCAACTTGTTGCTCGATCTCCAGCGTGAttt            >  2:296600/1‑148 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTg                                                                                                  <  2:318830/62‑1 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTg                                                                                                  >  1:318830/1‑62 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGAt                                                                 <  2:24558/95‑1 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGAt                                                                 >  1:24558/1‑95 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGtt                                                              <  5:260942/98‑1 (MQ=255)
                                                                                                                                       tttCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGtt                                                              >  6:260942/1‑98 (MQ=255)
                                                                                                                                            ggcggTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGTCTTTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCa       <  5:79480/148‑1 (MQ=255)
                                                                                                                                              cggTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCAtt     >  4:325567/1‑148 (MQ=255)
                                                                                                                                               ggTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTCTTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTg    <  8:18940/148‑1 (MQ=255)
                                                                                                                                                 tCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTg                                                                    <  2:392743/82‑1 (MQ=255)
                                                                                                                                                 tCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTg                                                                    >  1:392743/1‑82 (MQ=255)
                                                                                                                                                 tCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCg  >  6:289454/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCATCGGTGATTCGATTATCGAACCGCTGCGTGTACACGGTTTATTGCCAGGTAAAGACGCGGCTGCACGCGTTGCGTGGTTGCTGGAGCGCGTGGGCCTGTTACCTGAACATGCCTGGCGTTACCCGCATGAGTTTTCCGGCGGTCAGCGCCAGCGCATCTGCATTGCTCGCGCGTTGGCATTGAATCCAAAAGTGATCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCG  >  NC_000913/868775‑869067

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