Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 470,045 A→T 5.2% D392V (GAT→GTT)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,0450AT 5.2% 257.6 / 2.1 76D392V (GAT→GTT) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (37/35);  new base T (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAC  >  NC_000913/469898‑470183
                                                                                                                                                   |                                                                                                                                          
cGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGa                                                                                                                                            <  3:241894/148‑1 (MQ=255)
 gTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGat                                                                                                                                           <  7:186746/148‑1 (MQ=255)
  tATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCAGAGGGCCGACTGGTTCCGGCAAAAGTACCATGTTGTCTCTCATTAAGCGTAATTTCGACGTCATCGAGGGGGata                                                                                                                                          >  2:221439/1‑148 (MQ=255)
     ccGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTACGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTAGACGTCAGCGAGGGGGATATTc                                                                                                                                       >  1:94093/1‑148 (MQ=255)
        cAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCt                                                                                                                                    >  1:365039/1‑148 (MQ=255)
            cTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCa                                                                                                                                >  3:297275/1‑148 (MQ=255)
              gACCATCCTGCGCTGGAAAACGTCAATTTCTCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCAt                                                                                                                               >  1:130572/1‑147 (MQ=255)
                cgatgcTGCGCTGGAAAACGGCAATTTCGGGCTGAAACCCGGTCAGATGCTGGTTATCTGCGGGCCGACTGGTTCCGGCAAAAGTAGCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGat                                                                                                                            <  1:221439/143‑1 (MQ=255)
                   tCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATg                                                                                                                              >  1:214631/1‑143 (MQ=255)
                   tCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATg                                                                                                                              <  2:214631/143‑1 (MQ=255)
                    ccTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTc                                                                                                                        >  1:370676/1‑148 (MQ=255)
                    ccTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTc                                                                                                                        >  1:45814/1‑148 (MQ=255)
                    ccTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCATCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTc                                                                                                                        >  6:378575/1‑148 (MQ=255)
                       gcgcTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGTTATTCGCTTTCATGATATTCctc                                                                                                                     >  7:311666/1‑148 (MQ=255)
                         gcTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCTTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTg                                                                                                                   <  2:130572/148‑1 (MQ=255)
                          cTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGa                                                                                                                  >  4:119824/1‑148 (MQ=255)
                             gAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGa                                                                                                               <  2:365039/148‑1 (MQ=255)
                             gAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGa                                                                                                               <  2:3326/148‑1 (MQ=255)
                              aaaaCGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACg                                                                                                                >  8:41880/1‑146 (MQ=255)
                              aaaaCGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGaa                                                                                                              <  7:41880/148‑1 (MQ=255)
                                 aCGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGt                                                                                                            >  7:391220/1‑147 (MQ=255)
                                  cGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTa                                                                                                          <  1:165461/148‑1 (MQ=255)
                                    tCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCttt                                                                                                                                  <  1:36829/122‑1 (MQ=255)
                                    tCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCttt                                                                                                                                  >  2:36829/1‑122 (MQ=255)
                                      aaTTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCctc                                                                                                                     <  8:142575/133‑1 (MQ=255)
                                      aaTTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCctc                                                                                                                     >  7:142575/1‑133 (MQ=255)
                                         ttCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTc                                                                                                    >  7:332741/1‑147 (MQ=255)
                                         ttCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTc                                                                                                    <  8:332741/147‑1 (MQ=255)
                                          tCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCATCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCg                                                                                                   <  5:378575/147‑1 (MQ=255)
                                              ccTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGc                                                                                              <  3:119824/148‑1 (MQ=255)
                                                 gAAACCCGGTCAGATGCTGGGGATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTgg                                                                                           >  8:326977/1‑148 (MQ=255)
                                                  aaaCCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGa                                                                                                               <  6:174676/127‑1 (MQ=255)
                                                  aaaCCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGa                                                                                                               >  5:174676/1‑127 (MQ=255)
                                                  aaaCCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGc                                                                                          <  8:391220/148‑1 (MQ=255)
                                                     cccGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCg                                                                                         >  5:342929/1‑146 (MQ=255)
                                                       cGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTcc                                                                                                                       >  4:6937/1‑114 (MQ=255)
                                                       cGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTcc                                                                                                                       <  3:6937/114‑1 (MQ=255)
                                                       cGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTAc                                                                                                         <  2:106344/128‑1 (MQ=255)
                                                       cGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTAc                                                                                                         >  1:106344/1‑128 (MQ=255)
                                                       cGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAgc                                                                                     <  6:342929/148‑1 (MQ=255)
                                                               tGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGTTATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGc                                                                             >  4:257072/1‑148 (MQ=255)
                                                               tGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGc                                                                             >  2:77261/1‑148 (MQ=255)
                                                                         cTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGc                                                                   >  1:125206/1‑148 (MQ=255)
                                                                            cGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGACGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAg                                                                <  2:125206/148‑1 (MQ=255)
                                                                                ccGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGc                                                            <  2:94093/148‑1 (MQ=255)
                                                                                ccGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGc                                                            <  1:77261/148‑1 (MQ=255)
                                                                                 cGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGcc                                                           <  2:370676/148‑1 (MQ=255)
                                                                                    cTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCAtt                                                        >  3:21484/1‑148 (MQ=255)
                                                                                     tGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTc                                                       >  5:103951/1‑148 (MQ=255)
                                                                                      ggTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGTTAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTcc                                                      >  1:211492/1‑148 (MQ=255)
                                                                                      ggTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAgccgcc                                                                                 >  2:168552/1‑121 (MQ=255)
                                                                                      ggTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAgccgcc                                                                                 <  1:168552/121‑1 (MQ=255)
                                                                                            ggCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGTTATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCtttttt                                                <  8:311666/148‑1 (MQ=255)
                                                                                               aaaaGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGTTATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTg                                             <  3:257072/148‑1 (MQ=255)
                                                                                                aaaGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGa                                            <  6:103951/148‑1 (MQ=255)
                                                                                                aaaGTACCCTGTTGGCGCTGATTCAGGGTCATTTCGTCGTCAGGGAGGTGGATTTTCGCTTTCATGATATTCCTCTGACGTAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTg                                             <  2:386993/147‑1 (MQ=255)
                                                                                                  aGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGa                                            >  1:386993/1‑146 (MQ=255)
                                                                                                    tACACTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTATCCGCCTGGGGGAAGTTAGCCAGACGCCATTCCTTTTTTCTGACACt                                        <  6:304496/148‑1 (MQ=255)
                                                                                                     aCCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACtg                                       >  1:25431/1‑148 (MQ=255)
                                                                                                        cTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTg                                                                                            <  7:87069/92‑1 (MQ=255)
                                                                                                        cTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTg                                                                                            >  8:87069/1‑92 (MQ=255)
                                                                                                           ttgtCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCg                                  >  7:198329/1‑147 (MQ=255)
                                                                                                                gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATaa                             >  6:123238/1‑147 (MQ=255)
                                                                                                                gCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAAc                            <  7:326977/148‑1 (MQ=255)
                                                                                                                      tCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGa                                 <  6:173709/137‑1 (MQ=255)
                                                                                                                      tCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGa                                 >  5:173709/1‑137 (MQ=255)
                                                                                                                      tCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATcgcg                      >  8:56566/1‑148 (MQ=255)
                                                                                                                        aGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttttt                                                 <  2:177509/119‑1 (MQ=255)
                                                                                                                        aGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCttttt                                                 >  1:177509/1‑119 (MQ=255)
                                                                                                                        aGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCt                    <  4:21484/148‑1 (MQ=255)
                                                                                                                             cATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCAt                                                         >  4:303980/1‑106 (MQ=255)
                                                                                                                             cATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCAt                                                         <  3:303980/106‑1 (MQ=255)
                                                                                                                                ttCGTCGTCAGCGAGGGGGATATTCGCTTTGATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGCCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGAGAATAACATCGCGCTGGGTTGcc            >  5:373443/1‑148 (MQ=255)
                                                                                                                                ttCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGcc                                                                  <  7:220246/94‑1 (MQ=255)
                                                                                                                                ttCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGcc                                                                  >  8:220246/1‑94 (MQ=255)
                                                                                                                                   gACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGa         <  4:297275/148‑1 (MQ=255)
                                                                                                                                   gACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGa         <  5:123238/148‑1 (MQ=255)
                                                                                                                                   gACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGa         <  8:198329/148‑1 (MQ=255)
                                                                                                                                          gCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAc  >  2:278999/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
CGTATCCGCAGACTGACCATCCTGCGCTGGAAAACGTCAATTTCGCCCTGAAACCCGGTCAGATGCTGGGTATCTGCGGGCCGACTGGTTCCGGCAAAAGTACCCTGTTGTCGCTCATTCAGCGTCATTTCGACGTCAGCGAGGGGGATATTCGCTTTCATGATATTCCTCTGACGAAGTTACAACTCGATAGCTGGCGTAGCCGCCTGGCGGTAGTTAGCCAGACGCCATTCCTTTTTTCTGACACTGTGGCGAATAACATCGCGCTGGGTTGCCCGAATGCCAC  >  NC_000913/469898‑470183

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: