Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,303,758 A→G 5.3% S118P (TCG→CCG)  mdtN ← membrane fusion protein of efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,303,7580AG 5.3% 234.8 / 2.9 76S118P (TCG→CCG) mdtNmembrane fusion protein of efflux pump
Reads supporting (aligned to +/- strand):  ref base A (33/39);  new base G (2/2);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

TCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTT  >  NC_000913/4303616‑4303899
                                                                                                                                              |                                                                                                                                             
tCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCgg                                                                                                                                          >  8:374868/1‑148 (MQ=255)
 cTGCGGCGCTCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCGGAGACAAACCCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCgg                                                                                                                                          <  5:75559/147‑1 (MQ=255)
      gcgcgcTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCaaa                                                                                                                                   <  6:178781/149‑1 (MQ=255)
          gcTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGt                                                                                                                               <  7:336358/149‑1 (MQ=255)
           cTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGtt                                                                                                                              <  8:231963/149‑1 (MQ=255)
           cTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGtt                                                                                                                              <  4:299844/149‑1 (MQ=255)
           cTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGtt                                                                                                                              <  4:314080/149‑1 (MQ=255)
            tGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTg                                                                                                                             >  2:306232/1‑149 (MQ=255)
            tGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTg                                                                                                                             >  3:321830/1‑149 (MQ=255)
              cgcCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTgc                                                                                                                            <  8:381743/148‑1 (MQ=255)
                  gTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGt                                                                                                                        >  4:101618/1‑148 (MQ=255)
                     cTTGCACGGTCAACATCTTCCGCTTAAACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCaa                                                                                                                     <  5:167432/148‑1 (MQ=255)
                        gCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCAGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAacgc                                                                                                                  >  4:267896/1‑148 (MQ=255)
                        gCACGGTCAACATCTTCCGCGGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAac                                                                                                                    <  2:311005/146‑1 (MQ=255)
                                 aCATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGGTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGGCAACGCTACGCTGGGt                                                                                                        <  1:262930/149‑1 (MQ=255)
                                  cATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGt                                                                                                        >  1:13611/1‑148 (MQ=255)
                                      ttCCGCTGAGACAAAACCATCTTTCATTAATGGCTCGGTGCGGCGTAATGTATCTGTGCCCTGTTTCGCGGCGGCACGGGCATTTACTACCGCGGCATTAACCGACTCGGCACCAAACTGTTCCGCGTCAACGCTACGCTGGGTGAGCa                                                                                                   >  5:296247/1‑149 (MQ=255)
                                          gCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAAt                                                                                               >  2:353372/1‑149 (MQ=255)
                                           cTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAAtt                                                                                              <  7:373593/149‑1 (MQ=255)
                                            tGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAAttt                                                                                             >  5:80378/1‑149 (MQ=255)
                                              agaCAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGc                                                                                           <  3:129700/149‑1 (MQ=255)
                                                  aaaaCCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCCCGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCCCCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTAt                                                                                       >  2:329467/1‑149 (MQ=255)
                                                      ccTTCTTTACGTAACGGCTCGGTGCGGCGTAATGTATCTGTGGACTGTTGCGCGGCGGCACGGGCTTTTTCTACCTTGGCATTAACCGAGTCGGCACCAAACAGTTGCGCGGCAACGCTACGCTGGGTTAGCATAATTTGCTTATCCAg                                                                                   <  4:233629/149‑1 (MQ=255)
                                                       cttcttTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAgc                                                                                  >  3:338957/1‑149 (MQ=255)
                                                       cttcttTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAgc                                                                                  >  3:18962/1‑149 (MQ=255)
                                                       cttcttTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGAGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAgc                                                                                  >  1:216166/1‑149 (MQ=255)
                                                         tcttTCAGTAATGGCTCGGTGAGGCGTAATCTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAgcgc                                                                                >  6:121398/1‑149 (MQ=255)
                                                               aGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgag                                                                          <  4:242664/149‑1 (MQ=255)
                                                               aGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCga                                                                           <  3:101618/148‑1 (MQ=255)
                                                                 tAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgaggg                                                                        <  4:321830/149‑1 (MQ=255)
                                                                   aTGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCgagggag                                                                      <  2:356214/149‑1 (MQ=255)
                                                                     ggCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGc                                                                    >  7:362657/1‑149 (MQ=255)
                                                                     ggCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGc                                                                    >  7:357677/1‑149 (MQ=255)
                                                                           gTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGc                                                              >  8:303862/1‑149 (MQ=255)
                                                                            tGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCt                                                             >  5:145555/1‑149 (MQ=255)
                                                                              cggcgTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCttt                                                           <  5:339079/149‑1 (MQ=255)
                                                                               gtcttcaTGAATCTGTGGCCTTTTTCGCTGGGGCACGGTCTTTGTCTACCGTGGCATTACCCGAGTCGGCACCAAACAGTTGCGCGTCAACGCTACGCTGGGTTAGCATAATTTGCTTATCCAGCGGCGCGAGGGAGGCTTCAGCTTTc                                                          <  8:364449/143‑1 (MQ=255)
                                                                               ggcgTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTc                                                          <  2:13611/149‑1 (MQ=255)
                                                                                      tGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGAt                                                   <  6:193034/149‑1 (MQ=255)
                                                                                       gTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTc                                                          >  1:346568/1‑141 (MQ=255)
                                                                                       gTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTc                                                          <  2:346568/141‑1 (MQ=255)
                                                                                        tATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTg                                                 >  1:248601/1‑149 (MQ=255)
                                                                                        tATCTGTGGCCTGGTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGAtt                                                  <  7:340140/148‑1 (MQ=255)
                                                                                          tCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCtt                                                                  >  7:172080/1‑130 (MQ=255)
                                                                                          tCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCtt                                                                  <  8:172080/130‑1 (MQ=255)
                                                                                              tgGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGa                                                    >  2:112807/1‑140 (MQ=255)
                                                                                              tgGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGa                                                    <  1:112807/140‑1 (MQ=255)
                                                                                              tgGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTc                                            <  5:113039/148‑1 (MQ=255)
                                                                                                gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGGGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGc                                                        <  4:306892/134‑1 (MQ=255)
                                                                                                gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGGGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGc                                                        >  3:301898/1‑134 (MQ=255)
                                                                                                gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGGGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGc                                                        >  3:306892/1‑134 (MQ=255)
                                                                                                gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGGGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGc                                                        <  4:301898/134‑1 (MQ=255)
                                                                                                gCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTa                                         >  3:130719/1‑149 (MQ=255)
                                                                                                  cTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACg                                       <  5:129420/149‑1 (MQ=255)
                                                                                                  cTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACg                                       >  1:338112/1‑149 (MQ=255)
                                                                                                     tttCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACggg                                     <  4:130719/148‑1 (MQ=255)
                                                                                                      ttCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCTCCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGgcg                                   >  2:330215/1‑149 (MQ=255)
                                                                                                           ggcggcACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGt                              <  5:128344/149‑1 (MQ=255)
                                                                                                             cggcACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCg                            <  6:80378/149‑1 (MQ=255)
                                                                                                                 aCGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGc                        <  4:228140/149‑1 (MQ=255)
                                                                                                                 aCGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGc                        >  6:43707/1‑149 (MQ=255)
                                                                                                                   gggCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCg                       >  5:245278/1‑148 (MQ=255)
                                                                                                                    ggCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCAGCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGGCGATGCGGa                     <  7:374868/149‑1 (MQ=255)
                                                                                                                        ttttCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAg                 >  1:122906/1‑149 (MQ=255)
                                                                                                                         tttCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACa                  <  5:43707/147‑1 (MQ=255)
                                                                                                                         tttCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACa                  <  5:47366/147‑1 (MQ=255)
                                                                                                                          ttCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGc                <  3:267896/148‑1 (MQ=255)
                                                                                                                               ccGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCaaa             <  3:192185/146‑1 (MQ=255)
                                                                                                                                      agtAACCGAGTCGGAACCAAACTGTTGCCCCTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGAGCCGCGAGGGAGCCATCAGCTTTCGCCAGATTGGGTTCGAACGGGCGCGGCACGATTCCGAACAGCAAATCGCCCTCtt   >  2:80661/3‑149 (MQ=255)
                                                                                                                                      aTTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAAt            >  4:64824/1‑140 (MQ=255)
                                                                                                                                      aTTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAAt            <  3:64824/140‑1 (MQ=255)
                                                                                                                                       ttAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGttt  <  4:338957/149‑1 (MQ=255)
                                                                                                                                       ttAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGttt  >  7:78755/1‑149 (MQ=255)
                                                                                                                                       ttAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGttt  >  6:184118/1‑149 (MQ=255)
                                                                                                                                       ttAACCGAGTCGGCACCAAACTGTTGCGCGGCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTTGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGttt  <  6:251062/149‑1 (MQ=255)
                                                                                                                                        tAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGttt  >  1:232775/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                             
TCTGCGGCGCGCTGCGCCGTTCTTGCACGGTCAACATCTTCCGCTGAGACAAAACCTTCTTTCAGTAATGGCTCGGTGCGGCGTAATGTATCTGTGGCCTGTTTCGCGGCGGCACGGGCTTTTTCTACCGTGGCATTAACCGAGTCGGCACCAAACTGTTGCGCGTCAACGCTACGCTGGGTGAGCATAATTTGCTTATCCAGCGCCGCGAGGGAGGCTTCAGCTTTCGCCAGATTGGCTTCGTACGGGCGCGGGTCGATGCGGAACAGCAAATCGCCCTGTTT  >  NC_000913/4303616‑4303899

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: