Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,447 A→G 5.1% I250T (ATA→ACA)  cydC ← glutathione/cysteine ABC transporter export permease/ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,4470AG 5.1% 407.2 / 7.7 117I250T (ATA→ACA) cydCglutathione/cysteine ABC transporter export permease/ATPase
Reads supporting (aligned to +/- strand):  ref base A (57/54);  new base G (3/3);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 120 total aligned reads displayed.

CCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTG  >  NC_000913/928301‑928595
                                                                                                                                                  |                                                                                                                                                    
ccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATc                                                                                                                                                    >  3:1147077/1‑149 (MQ=255)
 cGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCg                                                                                                                                                   >  1:931136/1‑149 (MQ=255)
        ggTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGa                                                                                                                                            >  4:753203/1‑149 (MQ=255)
              aGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGc                                                                                                                                      >  4:935600/1‑149 (MQ=255)
               gTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCg                                                                                                                                     <  3:1349976/149‑1 (MQ=255)
               gTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCg                                                                                                                                     >  4:892950/1‑149 (MQ=255)
                tGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATGctg                                                                                                                                    <  2:836549/149‑3 (MQ=255)
                   ttCAAACGCGGCTAACGCTCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTa                                                                                                                                                      >  1:423364/1‑128 (MQ=255)
                   ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTa                                                                                                                                                      <  2:423364/128‑1 (MQ=255)
                   ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATg                                                                                                                                       >  2:87175/1‑143 (MQ=255)
                   ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTGTCGCAAGCGATAATg                                                                                                                                       <  1:87175/143‑1 (MQ=255)
                   ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGt                                                                                                                                   >  1:1007142/1‑147 (MQ=255)
                   ttCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCa                                                                                                                                 >  4:117547/1‑149 (MQ=255)
                    tCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAg                                                                                                                                >  3:83052/1‑149 (MQ=255)
                     cAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGt                                                                                                                               >  1:1341894/1‑149 (MQ=255)
                       aaCGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTc                                                                                                                             >  4:455389/1‑149 (MQ=255)
                          gcgGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGa                                                                                                                          <  1:412139/149‑1 (MQ=255)
                          gcgGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGa                                                                                                                          <  1:322923/149‑1 (MQ=255)
                             gCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTg                                                                                                                       >  1:565348/1‑149 (MQ=255)
                              cTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGAt                                                                                                                         <  3:1061141/146‑1 (MQ=255)
                              cTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGAAAATGCGGTCAGTTCAGAt                                                                                                                         >  4:1061141/1‑146 (MQ=255)
                               tAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGAc                                                                                                                     <  3:1317709/149‑1 (MQ=255)
                                       agaagaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGa                                                                                                                          >  1:356392/1‑136 (MQ=255)
                                       agaagaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGa                                                                                                                          <  2:356392/136‑1 (MQ=255)
                                       agaagaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTgc                                                                                                             <  3:117547/149‑1 (MQ=255)
                                         aagaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTgcgc                                                                                                           <  2:209098/149‑1 (MQ=255)
                                         aagaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTgcgc                                                                                                           <  3:892950/149‑1 (MQ=255)
                                          agaCAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCt                                                                                                          >  1:835740/1‑149 (MQ=255)
                                              aaaCAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCa                                                                                                      >  4:334711/1‑149 (MQ=255)
                                                  aGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCa                                                                                                  >  3:1279831/1‑149 (MQ=255)
                                                    ggCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAtt                                                                                                >  4:1039744/1‑149 (MQ=255)
                                                       aaTTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGTTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGaa                                                                                             <  3:334711/149‑1 (MQ=255)
                                                            aCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAt                                                                                                 <  3:59391/140‑1 (MQ=255)
                                                            aCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAt                                                                                                 >  4:59391/1‑140 (MQ=255)
                                                            aCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCt                                                                                        <  2:1341894/149‑1 (MQ=255)
                                                              gcgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGt                                                                                      >  1:285333/1‑149 (MQ=255)
                                                               cgcCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTa                                                                                     <  3:455389/149‑1 (MQ=255)
                                                                    ggTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtctc                                                                                <  3:776371/149‑1 (MQ=255)
                                                                      ttGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTa                                                                              >  4:1375109/1‑149 (MQ=255)
                                                                       tGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAg                                                                             >  4:1463513/1‑149 (MQ=255)
                                                                         aGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGt                                                                                                                >  4:96597/1‑112 (MQ=255)
                                                                         aGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGt                                                                                                                <  3:96597/112‑1 (MQ=255)
                                                                          gCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAt                                                                                            >  2:982495/1‑131 (MQ=255)
                                                                          gCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAt                                                                                            <  1:982495/131‑1 (MQ=255)
                                                                           cATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAg                                                                                                                                <  4:489933/94‑1 (MQ=255)
                                                                            aTTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcg                                                                        >  1:478275/1‑149 (MQ=255)
                                                                             ttGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACaa                                                                                                                                         >  4:329273/1‑83 (MQ=255)
                                                                             ttGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACaa                                                                                                                                         <  3:329273/83‑1 (MQ=255)
                                                                                    cAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATaa                                                                <  4:1279831/149‑1 (MQ=255)
                                                                                     aaCGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAAc                                                               <  3:935600/149‑1 (MQ=255)
                                                                                      aCGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACg                                                              <  4:83052/149‑1 (MQ=255)
                                                                                        gccgccAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGAt                                                            >  4:133869/1‑149 (MQ=255)
                                                                                             cAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTg                                                       <  3:1432600/149‑1 (MQ=255)
                                                                                               gACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGc                                                     >  1:1122528/1‑149 (MQ=255)
                                                                                                aCGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCAtt                                                                                                >  1:800049/1‑105 (MQ=255)
                                                                                                aCGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCAAGCCAtt                                                                                                <  2:800049/105‑1 (MQ=255)
                                                                                                  gCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCAACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCt                                                                                        <  1:807242/111‑1 (MQ=255)
                                                                                                     aTCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAAtt                                                                                           >  4:142104/1‑105 (MQ=255)
                                                                                                     aTCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCaaa                                               >  1:818097/1‑149 (MQ=255)
                                                                                                         acaGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtctc                                                                                >  3:1091728/1‑112 (MQ=255)
                                                                                                         acaGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATtctc                                                                                <  4:1091728/112‑1 (MQ=255)
                                                                                                         acaGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATg                                           >  3:396554/1‑149 (MQ=255)
                                                                                                             cATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtca                                       <  3:1375109/149‑1 (MQ=255)
                                                                                                              aTCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcag                                      >  2:391728/1‑149 (MQ=255)
                                                                                                                cAGGATCACCGCTAACGCGCCAATGAGTAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagct                                    <  3:1463513/149‑1 (MQ=255)
                                                                                                                  ggATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagctca                                  >  3:727976/1‑149 (MQ=255)
                                                                                                                     tCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTa                                                                              <  3:609805/102‑1 (MQ=255)
                                                                                                                     tCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTa                                                                              >  4:609805/1‑102 (MQ=255)
                                                                                                                        ccGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAg                                                                             <  1:709673/100‑1 (MQ=255)
                                                                                                                        ccGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAg                                                                             >  2:709673/1‑100 (MQ=255)
                                                                                                                        ccGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGc                            <  4:727976/149‑1 (MQ=255)
                                                                                                                         cGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGAtt                                                                                                                        >  2:577514/1‑56 (MQ=255)
                                                                                                                         cGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGt                                                                                                                >  3:144722/1‑64 (MQ=255)
                                                                                                                         cGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGt                                                                                                                <  4:144722/64‑1 (MQ=255)
                                                                                                                            tAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCtt                        >  1:7220/1‑149 (MQ=255)
                                                                                                                             aaCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTg                       >  2:117162/1‑149 (MQ=255)
                                                                                                                              aCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                      <  2:338485/149‑1 (MQ=255)
                                                                                                                              aCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                      <  2:285333/149‑1 (MQ=255)
                                                                                                                              aCGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGc                      <  2:835740/149‑1 (MQ=255)
                                                                                                                               cgcgCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagct                                    <  2:608126/134‑1 (MQ=255)
                                                                                                                               cgcgCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcagct                                    >  1:608126/1‑134 (MQ=255)
                                                                                                                               cgcgCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCa                     <  2:7220/149‑1 (MQ=255)
                                                                                                                                gcgcCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCt                         >  2:808511/1‑144 (MQ=255)
                                                                                                                                 cgcCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGc                   <  2:818097/149‑1 (MQ=255)
                                                                                                                                   ccAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCa                 >  3:420646/1‑149 (MQ=255)
                                                                                                                                    cAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATc                                                           <  3:566428/106‑1 (MQ=255)
                                                                                                                                     aaTGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAg                                                                             <  4:1207041/87‑1 (MQ=255)
                                                                                                                                     aaTGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAg                                                                             >  3:1207041/1‑87 (MQ=255)
                                                                                                                                      aTGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATgg                                          >  4:1193075/1‑121 (MQ=255)
                                                                                                                                          gcagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc          >  3:277786/1‑149 (MQ=255)
                                                                                                                                          gcagcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTc          >  3:29147/1‑149 (MQ=255)
                                                                                                                                             gcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcag                                      >  4:659747/1‑118 (MQ=255)
                                                                                                                                             gcaTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGtcag                                      <  3:659747/118‑1 (MQ=255)
                                                                                                                                              caTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt      <  2:1122528/149‑1 (MQ=255)
                                                                                                                                              caTTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt      <  3:803337/149‑1 (MQ=255)
                                                                                                                                              caTTATCGCTTGCGACAATGAGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgtt      <  4:420646/149‑1 (MQ=255)
                                                                                                                                               aTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCa                                                                                                      <  2:324358/52‑1 (MQ=25)
                                                                                                                                               aTTGTCGCAAGCGATAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAgttg     >  4:1205718/1‑149 (MQ=255)
                                                                                                                                                ttATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTtgcgtgc                                                                     <  3:982333/84‑1 (MQ=255)
                                                                                                                                                  aTCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGttgttg  <  4:953664/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATCGCTTGCGACAATGCGGTCAGTTCAGATTGACGGCGTTGCGCTTCCAGCCATTGAATTTCTGTATTCTCTAGTTGCGTGCGATAACGATCGCTGGCACCAAAAATGGTCAGCTCAGCTTGCCCTTGCAGCCAGGCCGTCAGTTGTTG  >  NC_000913/928301‑928595

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