Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,158,885 C→G 5.3% S167W (TCG→TGG) ‡ mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,158,8850CG 5.3% 251.4 / 4.9 75S167W (TCG→TGG) ‡mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base C (43/28);  new base G (2/2);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GCTGCCCAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAAC  >  NC_000913/2158739‑2159031
                                                                                                                                                  |                                                                                                                                                  
gctgcCCAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGc                                                                                                                                                  >  2:1460592/1‑149 (MQ=255)
gctgcCCAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTAAGATCCTCACGCTGACGTCCGATACTTATGCGCAGGGTGAACTCTACTATATCGCCTCGACGCAGCTTGCTCAGACGATTTCGc                                                                                                                                                  >  2:1460600/1‑149 (MQ=255)
    cccAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAAt                                                                                                                                              >  1:1435112/1‑149 (MQ=255)
    cccAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAAt                                                                                                                                              >  1:1515675/1‑149 (MQ=255)
     ccAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATc                                                                                                                                             >  2:1123368/1‑149 (MQ=255)
      cAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCg                                                                                                                                            >  1:1273306/1‑149 (MQ=255)
       aGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGTCGATTTGCGAAATCgt                                                                                                                                           <  2:1516645/149‑2 (MQ=255)
         tGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACg                                                                                                                                         >  2:1064982/1‑149 (MQ=255)
           ggATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGTCGATTTGCGAAATCGTCgga                                                                                                                                       <  1:873400/149‑2 (MQ=255)
                cccAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCCCGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTgg                                                                                                                                  >  2:1437353/1‑149 (MQ=255)
                 ccAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCCCGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGt                                                                                                                                 <  1:1086865/149‑1 (MQ=255)
                  cAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTg                                                                                                                                >  2:508842/1‑149 (MQ=255)
                    gccgccCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATCCTAATTCGCAGGGGGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATATCGCAAATCGACGGTGTTGGTgat                                                                                                                              <  2:815944/149‑1 (MQ=255)
                        cccGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcg                                                                                                                          >  2:1369477/1‑149 (MQ=255)
                        cccGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcg                                                                                                                          >  2:176335/1‑149 (MQ=255)
                         ccGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcga                                                                                                                         <  1:1343935/149‑1 (MQ=255)
                         ccGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTgatgtcga                                                                                                                         >  2:1464490/1‑149 (MQ=255)
                              cTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCGGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGatgtcgatgtcg                                                                                                                    <  2:1124257/149‑1 (MQ=255)
                                aTCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGa                                                                                                                  >  2:1360713/1‑149 (MQ=255)
                                aTCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGa                                                                                                                  >  2:543514/1‑149 (MQ=255)
                                 tCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGCCGGTGTTGGTGATGTCGATGTCGGAg                                                                                                                 >  1:376332/1‑149 (MQ=255)
                                  cGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGAGGGTGTTGGTGATGTCGATGTCGGAgg                                                                                                                >  2:50215/1‑149 (MQ=255)
                                   gCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGc                                                                                                               >  2:866254/1‑149 (MQ=255)
                                   gCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGc                                                                                                               >  1:1482699/1‑149 (MQ=255)
                                    cAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCa                                                                                                              >  2:711694/1‑149 (MQ=255)
                                     aaaGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAg                                                                                                             <  1:343031/149‑1 (MQ=255)
                                          gAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCAc                                                                                                        <  1:404735/149‑1 (MQ=255)
                                          gAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTAATCGCAGGGGGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCAc                                                                                                        <  1:373040/149‑1 (MQ=255)
                                            aCCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTg                                                                                                      <  2:230775/149‑1 (MQ=255)
                                                   ggATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCccgccg                                                                                               <  2:1395556/149‑1 (MQ=255)
                                                    gATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGt                                                                                              >  2:307720/1‑149 (MQ=255)
                                                     aTGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTa                                                                                             <  1:426996/149‑1 (MQ=255)
                                                      tGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTAc                                                                                            <  1:1312058/149‑1 (MQ=255)
                                                        cgcCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTAGATGTCGGAGGCAGCTCACTGCCCGCCGTAcgc                                                                                          >  2:1504197/1‑149 (MQ=255)
                                                         gcCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTAcgcg                                                                                         >  1:1547857/1‑149 (MQ=255)
                                                          ccAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGt                                                                                        >  2:717651/1‑149 (MQ=255)
                                                           cAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTc                                                                                       >  1:1389075/1‑149 (MQ=255)
                                                           cAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTc                                                                                       >  1:295155/1‑149 (MQ=255)
                                                            aaTTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCg                                                                                      >  2:738503/1‑149 (MQ=255)
                                                             aTTATGATCCTCACGCTGACGTCCGATACTTAATCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCgg                                                                                     <  1:1216631/149‑1 (MQ=255)
                                                               tATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGc                                                                                   <  1:109849/149‑1 (MQ=255)
                                                                 tGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGcag                                                                                 <  2:1482699/149‑3 (MQ=255)
                                                                    tCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAAt                                                                              >  1:1334591/1‑149 (MQ=255)
                                                                       tCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCg                                                                           >  1:379729/1‑149 (MQ=255)
                                                                         aCGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCa                                                                         >  1:757153/1‑149 (MQ=255)
                                                                         aCGCTGACGTACGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTAGACGCAGCAGGCTCCGACGATTTCGCAAATCGAAGGTGATGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCTTCGGGCTGAATCCGcc                                                                         >  2:62517/1‑148 (MQ=255)
                                                                               aCGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCt                                                                   <  1:1540647/149‑1 (MQ=255)
                                                                                 gTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGt                                                                 <  1:147118/149‑1 (MQ=255)
                                                                                  tCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCGCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGtt                                                                >  2:775700/1‑149 (MQ=255)
                                                                                    cGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTa                                                              >  2:486673/1‑149 (MQ=255)
                                                                                      aTACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAAt                                                            <  2:1515675/149‑1 (MQ=255)
                                                                                      aTACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCGGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAAt                                                            <  1:107656/149‑1 (MQ=255)
                                                                                          ttATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAgg                                                        >  1:1143026/1‑149 (MQ=255)
                                                                                           tATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAggg                                                       <  2:424328/149‑1 (MQ=255)
                                                                                                gCAGGGTTAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCgtgt                                                  <  2:1172498/149‑1 (MQ=255)
                                                                                                    ggTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCt                                              >  2:37831/1‑149 (MQ=255)
                                                                                                          cTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGgacga                                        <  1:176335/149‑1 (MQ=255)
                                                                                                           tGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGgacgac                                       >  2:801298/1‑149 (MQ=255)
                                                                                                           tGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGgacgac                                       >  1:1031442/1‑149 (MQ=255)
                                                                                                            gTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGacgacg                                      >  1:683062/1‑149 (MQ=255)
                                                                                                             tACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGt                                     >  1:918650/1‑149 (MQ=255)
                                                                                                                gATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACg                                  >  1:1430426/1‑149 (MQ=255)
                                                                                                                       ccTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGCTGTCCGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCAccgcc                           >  1:87970/1‑149 (MQ=255)
                                                                                                                       ccTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCAccgcc                           <  1:775700/149‑1 (MQ=255)
                                                                                                                        cTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCAccgccg                          >  2:1365422/1‑149 (MQ=255)
                                                                                                                           gACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCa                       <  2:1389075/149‑1 (MQ=255)
                                                                                                                                aGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCGCCGCCGTCAGCAAt                  >  1:1314304/1‑149 (MQ=255)
                                                                                                                                      cTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAAc            >  2:1104936/1‑149 (MQ=255)
                                                                                                                                        ccGTCGATTTGCGAAATCGTCGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCCGCAATGCCAACGt          >  2:1074138/1‑149 (MQ=255)
                                                                                                                                        ccGTCGATTTGCGAAATCGTCGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGt          >  2:1079686/1‑149 (MQ=255)
                                                                                                                                        ccGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCGGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGt          <  1:1360713/149‑1 (MQ=255)
                                                                                                                                        ccGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGt          <  1:1464490/149‑1 (MQ=255)
                                                                                                                                          gacgaTTTCGCAAATCGACGGTGGTGGTGCTGTCGATGGCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGc        <  1:1437353/149‑1 (MQ=255)
                                                                                                                                                ttCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAAc  <  1:1064982/149‑1 (MQ=255)
                                                                                                                                                ttCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAAc  <  1:1027906/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GCTGCCCAGTGGGATGCCCAGCCGCCCGACCTATCGCAAAGCGAACCCGTCGGATGCGCCAATTATGATCCTCACGCTGACGTCCGATACTTATTCGCAGGGTGAACTGTACGATTTCGCCTCGACGCAGCTGGCTCCGACGATTTCGCAAATCGACGGTGTTGGTGATGTCGATGTCGGAGGCAGCTCACTGCCCGCCGTACGCGTCGGGCTGAATCCGCAGGCGCTGTTTAATCAGGGCGTGTCGCTGGACGACGTACGCACCGCCGTCAGCAATGCCAACGTGCGTAAAC  >  NC_000913/2158739‑2159031

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: