Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 470,515 G→A 5.0% A549T (GCC→ACC)  mdlA → putative multidrug ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913470,5150GA 5.0% 268.7 / 2.7 79A549T (GCC→ACC) mdlAputative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base A (2/2);  total (42/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 = 9.79e-01

GCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGG  >  NC_000913/470374‑470659
                                                                                                                                             |                                                                                                                                                
gcggTGGACGGACGCACTGAGCACCAGCTCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGc                                                                                                                                            <  4:4099/148‑1 (MQ=255)
         ggacgCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATc                                                                                                                                   <  7:334356/148‑1 (MQ=255)
                cTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGAtgt                                                                                                                            >  3:135816/1‑148 (MQ=255)
                    gCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCatg                                                                                                                                <  3:113726/140‑1 (MQ=255)
                    gCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCatg                                                                                                                                >  4:113726/1‑140 (MQ=255)
                      aCCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTgg                                                                                                                       >  8:210955/1‑147 (MQ=255)
                                gCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGc                                                                                                            >  3:334351/1‑148 (MQ=255)
                                   tAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGc                                                                                                         >  5:279379/1‑148 (MQ=255)
                                    aaCCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCt                                                                                                        >  7:262523/1‑148 (MQ=255)
                                    aaCCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCt                                                                                                        >  8:225159/1‑148 (MQ=255)
                                              aGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGa                                                                                              <  3:353942/148‑1 (MQ=255)
                                               gTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCGAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCTCGat                                                                                             <  7:336314/148‑1 (MQ=255)
                                                          gaagaaCGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCACCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTa                                                                                  >  7:356939/1‑148 (MQ=255)
                                                             gaaCGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCacaaca                                                                                                                >  6:201824/1‑115 (MQ=255)
                                                             gaaCGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCacaaca                                                                                                                <  5:201824/115‑1 (MQ=255)
                                                                cGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaaca                                                                            <  6:279379/148‑1 (MQ=255)
                                                                cGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaaca                                                                            <  4:334351/148‑1 (MQ=255)
                                                                      tcatcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCatg                                                                                                                                <  3:260251/90‑1 (MQ=255)
                                                                      tcatcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCatg                                                                                                                                >  4:260251/1‑90 (MQ=255)
                                                                       catcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAg                                                                     >  6:266370/1‑148 (MQ=255)
                                                                        atcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAg                                                                     <  8:262523/147‑1 (MQ=255)
                                                                         tcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTa                                                                                        >  1:221971/1‑127 (MQ=255)
                                                                         tcaGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTa                                                                                        <  2:221971/127‑1 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACACAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTcgac                                                         >  5:167831/1‑147 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaaca                                                                            >  7:257443/1‑128 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaaca                                                                            <  8:257443/128‑1 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAg                                                                     >  3:192244/1‑135 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAg                                                                     <  4:192244/135‑1 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTcgacg                                                        >  3:23252/1‑148 (MQ=255)
                                                                                    cGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTcgac                                                         >  2:301298/1‑147 (MQ=255)
                                                                                         ttCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGt                                                                                                     >  6:163557/1‑98 (MQ=255)
                                                                                         ttCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGt                                                                                                     <  5:163557/98‑1 (MQ=255)
                                                                                          tCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTc                                                  >  2:108163/1‑148 (MQ=255)
                                                                                            tGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTg                                                                                                       >  3:63762/1‑93 (MQ=255)
                                                                                            tGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTg                                                                                                       <  4:63762/93‑1 (MQ=255)
                                                                                               aCTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGa                                              >  3:259256/1‑147 (MQ=255)
                                                                                                 tGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGaaaa                                           >  7:43245/1‑148 (MQ=255)
                                                                                                 tGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGaaaa                                           >  4:235591/1‑148 (MQ=255)
                                                                                                 tGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGaaaa                                           >  2:369672/1‑148 (MQ=255)
                                                                                                 tGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGaaaa                                           >  2:364179/1‑148 (MQ=255)
                                                                                                  gACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaacaa                                                                           <  3:18534/115‑1 (MQ=255)
                                                                                                  gACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATcaacaa                                                                           >  4:18534/1‑115 (MQ=255)
                                                                                                  gACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAAt                                          >  2:384769/1‑148 (MQ=255)
                                                                                                   aCGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAAt                                          <  4:259256/147‑1 (MQ=255)
                                                                                                     ggAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgc                                       <  5:266370/148‑1 (MQ=255)
                                                                                                      gAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      <  7:225159/148‑1 (MQ=255)
                                                                                                       aaGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCga                                     <  1:108163/148‑1 (MQ=255)
                                                                                                         gCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCgagg                                   <  1:369672/148‑1 (MQ=255)
                                                                                                         gCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCgagg                                   <  1:364179/148‑1 (MQ=255)
                                                                                                              tGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGcc                              >  5:321825/1‑148 (MQ=255)
                                                                                                               gAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      <  6:309001/139‑1 (MQ=255)
                                                                                                               gAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      >  5:309001/1‑139 (MQ=255)
                                                                                                               gAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGc                               >  6:277559/1‑146 (MQ=255)
                                                                                                               gAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGc                               <  5:277559/146‑1 (MQ=255)
                                                                                                                aaaTTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGt                            >  5:105466/1‑148 (MQ=255)
                                                                                                                  attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCg                          <  8:43245/148‑1 (MQ=255)
                                                                                                                  attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCg                          <  1:384769/148‑1 (MQ=255)
                                                                                                                  attattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCg                          <  7:210955/148‑1 (MQ=255)
                                                                                                                  attattGTGATGCAGCACGGACATATCACCCAGCGGGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCg                          <  8:356939/148‑1 (MQ=255)
                                                                                                                  attatgTTGATGCAGCAAGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGAGGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCTCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCg                          <  6:167831/148‑1 (MQ=255)
                                                                                                                   ttattGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCGCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGa                         <  1:301298/148‑1 (MQ=255)
                                                                                                                    tattGTGATGCAGCACGGACATATCACCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGt                            <  8:44557/144‑1 (MQ=255)
                                                                                                                    tattGTGATGCAGCACGGACATATCACCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGt                            >  7:44557/1‑144 (MQ=255)
                                                                                                                            tGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGt                 <  3:235591/147‑1 (MQ=255)
                                                                                                                                 caCGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGc           >  5:133749/1‑148 (MQ=255)
                                                                                                                                  agggACATATCGCCCAGCGTTGCAATCATTATGTGCTGGCACAACAAAGCGGCTGGTATCGCTATATGAATCGCTAGCAACAACTGCAGGCGGCGCTCTACGCCGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGcc          <  2:299179/146‑1 (MQ=255)
                                                                                                                                    ggACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCaa        >  3:148880/1‑148 (MQ=255)
                                                                                                                                       cATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACtg     >  5:269492/1‑148 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGCCGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCTTAGTTTTAGCCAACTGTgg  >  3:65426/1‑148 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      >  4:126705/1‑112 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      >  7:123003/1‑112 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      <  3:126705/112‑1 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgcg                                      <  8:123003/112‑1 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgc                                       <  2:83520/111‑1 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATcgc                                       >  1:83520/1‑111 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGtt                <  4:145142/134‑1 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGtt                >  3:145142/1‑134 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTgg  >  2:145193/1‑148 (MQ=255)
                                                                                                                                          atCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTgg  >  2:24815/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGG  >  NC_000913/470374‑470659

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: