Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 784,012 T→A 6.1% Q271L (CAG→CTG)  zitB ← zinc efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913784,0120TA 6.1% 982.2 / 36.8 264Q271L (CAG→CTG) zitBzinc efflux system
Reads supporting (aligned to +/- strand):  ref base T (128/120);  new base A (9/7);  total (137/127)
Fisher's exact test for biased strand distribution p-value = 8.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.98e-01

Only 100 of 264 total aligned reads displayed.

CACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCG  >  NC_000913/783866‑784158
                                                                                                                                                  |                                                                                                                                                  
cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   <  1:805808/148‑1 (MQ=255)
cACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   >  4:1124360/1‑148 (MQ=255)
      cgcTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaa                                                                                                                                             <  2:649843/148‑1 (MQ=255)
                ttAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGa                                                                                                                                                  >  3:230802/1‑133 (MQ=255)
                  aaTGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtg                                                                                                                                 <  2:764727/148‑1 (MQ=255)
                                 aaTGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaac                                                                                                                                            >  3:393115/1‑122 (MQ=255)
                                      ccAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgat                                                                                                                               >  4:801019/1‑130 (MQ=255)
                                          acacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   >  3:1015644/1‑106 (MQ=255)
                                          acacCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   <  4:1015644/106‑1 (MQ=255)
                                             ccccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   >  4:519431/1‑103 (MQ=255)
                                             ccccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTg                                                                                                            >  3:51703/1‑142 (MQ=255)
                                               ccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtga                                                                                                                                >  4:460060/1‑120 (MQ=255)
                                               ccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGc                                                                                                     >  4:299084/1‑147 (MQ=255)
                                               ccTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                    <  4:495744/148‑1 (MQ=255)
                                                 ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   >  4:83435/1‑99 (MQ=255)
                                                 ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   <  3:83435/99‑1 (MQ=255)
                                                 ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtga                                                                                                                                >  3:34561/1‑118 (MQ=255)
                                                 ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCAcc                                                                                                              >  4:391486/1‑136 (MQ=255)
                                                 ttcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAAc                                                                                                  >  2:471019/1‑148 (MQ=255)
                                                  tcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                                                                >  2:505247/1‑133 (MQ=255)
                                                  tcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                                                                <  1:1121287/133‑1 (MQ=255)
                                                  tcattcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACg                                                                                                 >  4:1089027/1‑148 (MQ=255)
                                                     ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGAt                                                                                                                                                 >  3:800543/1‑97 (MQ=255)
                                                     ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                                                                              >  1:612824/1‑148 (MQ=255)
                                                     ttcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtca                                                                                              <  4:172361/148‑1 (MQ=255)
                                                      tcaGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGtcat                                                                                             >  4:334648/1‑148 (MQ=255)
                                                         ggTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtg                                                                                                                                 >  1:968079/1‑109 (MQ=255)
                                                            ggCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCgg                                                                                                                     >  1:393297/1‑118 (MQ=255)
                                                             gCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGa                                                                                                                   >  1:1012910/1‑119 (MQ=255)
                                                                gTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                      >  2:114491/1‑113 (MQ=255)
                                                                gTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                      <  1:114491/113‑1 (MQ=255)
                                                                 tCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTtct                                                                                  >  4:633693/1‑148 (MQ=255)
                                                                  cTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcgtg                                                                                                                           >  2:838755/1‑106 (MQ=255)
                                                                    ggCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGc                                                                               <  3:636627/148‑1 (MQ=255)
                                                                      cccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                      <  3:6831/107‑1 (MQ=255)
                                                                       ccATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCgg                                                                                                                     <  2:396048/107‑1 (MQ=255)
                                                                        cATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGa                                                                                                                   >  2:29634/1‑108 (MQ=255)
                                                                        cATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                                                                              <  3:131046/145‑1 (MQ=255)
                                                                         aTGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAcc                                                                          <  1:471019/148‑1 (MQ=255)
                                                                           gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   <  2:1095263/73‑1 (MQ=255)
                                                                           gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTg                                                                                                                                                   >  1:1095263/1‑73 (MQ=255)
                                                                           gACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCt                                                                                     <  3:368465/135‑1 (MQ=255)
                                                                            aCAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatc                                                                       >  4:187104/1‑148 (MQ=255)
                                                                             cAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatcc                                                                      >  4:426817/1‑148 (MQ=255)
                                                                               aGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccatccat                                                                    >  3:488019/1‑148 (MQ=255)
                                                                                  ttGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgat                                                                                                                               <  2:501576/86‑1 (MQ=255)
                                                                                  ttGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                           >  3:1120691/1‑138 (MQ=255)
                                                                                  ttGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                 >  2:870657/1‑148 (MQ=255)
                                                                                   tGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacat                                                                <  3:453587/148‑1 (MQ=255)
                                                                                    gATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatg                                                               >  4:929213/1‑148 (MQ=255)
                                                                                        ttCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtaca                                                           <  3:187104/148‑1 (MQ=255)
                                                                                         tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCa                                                                                                                >  1:915042/1‑94 (MQ=255)
                                                                                         tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacat                                                          >  1:398976/1‑148 (MQ=255)
                                                                                         tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtacat                                                          <  1:496094/148‑1 (MQ=255)
                                                                                         tCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtaca                                                                 >  1:888365/1‑141 (MQ=255)
                                                                                           cATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgta                                                             >  2:722898/1‑143 (MQ=255)
                                                                                            aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCt                                                                             >  4:480084/1‑126 (MQ=255)
                                                                                            aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCt                                                                             <  3:480084/126‑1 (MQ=255)
                                                                                            aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGt                                                       >  1:1042108/1‑148 (MQ=255)
                                                                                            aTCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGt                                                       <  3:122023/148‑1 (MQ=255)
                                                                                              cTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCg                                                                                                                      <  4:891078/83‑1 (MQ=255)
                                                                                                 aaTGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCatcgtgatcg                                                                                                                             >  4:685547/1‑73 (MQ=255)
                                                                                                     gtggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGc                                              <  2:574905/148‑1 (MQ=255)
                                                                                                       ggCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCgg                                            <  4:488019/148‑1 (MQ=255)
                                                                                                              tCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAtacatgtaca                                                           <  4:122519/126‑1 (MQ=255)
                                                                                                                 aTCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCcaa                                                                                                                                             >  1:862999/1‑41 (MQ=255)
                                                                                                                   cTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGAtt                                <  3:929213/148‑1 (MQ=255)
                                                                                                                     gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGc                                              >  3:284273/1‑132 (MQ=255)
                                                                                                                     gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTc                              >  4:672147/1‑148 (MQ=255)
                                                                                                                     gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTc                              >  1:510261/1‑148 (MQ=255)
                                                                                                                     gATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTc                              >  2:416040/1‑148 (MQ=255)
                                                                                                                       tAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTccc                            <  1:696202/148‑1 (MQ=255)
                                                                                                                        aGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCg                           <  4:612821/148‑1 (MQ=255)
                                                                                                                         gTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCgg                          >  3:632713/1‑148 (MQ=255)
                                                                                                                           gATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATgg                                                        <  4:395548/116‑1 (MQ=255)
                                                                                                                            aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                    <  2:213291/71‑1 (MQ=255)
                                                                                                                            aTCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCa                                                                                                    >  1:213291/1‑71 (MQ=255)
                                                                                                                             tCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATgg                                                        >  2:112760/1‑114 (MQ=255)
                                                                                                                                aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                                                                   <  4:261827/100‑1 (MQ=255)
                                                                                                                                aTCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTcc                             >  2:41031/1‑138 (MQ=255)
                                                                                                                                  cAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccat                                                                        <  3:745541/93‑1 (MQ=255)
                                                                                                                                  cAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAccat                                                                        >  4:745541/1‑93 (MQ=255)
                                                                                                                                   aGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATacgac                <  1:416040/148‑1 (MQ=255)
                                                                                                                                     gTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGAttt                               >  3:318225/1‑131 (MQ=255)
                                                                                                                                      ccttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTg                                               >  3:1041796/4‑114 (MQ=255)
                                                                                                                                      ccttgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTg                                               <  4:1041796/111‑1 (MQ=255)
                                                                                                                                        gtgtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc           <  2:888749/148‑1 (MQ=255)
                                                                                                                                        gtgtTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttc           <  3:897177/148‑1 (MQ=255)
                                                                                                                                         tgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCg                                             >  3:1083264/1‑113 (MQ=255)
                                                                                                                                         tgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTcc                                      <  1:501600/120‑1 (MQ=255)
                                                                                                                                         tgtTGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGAt                                 <  3:576454/125‑1 (MQ=255)
                                                                                                                                           ttGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagt        >  1:289609/1‑148 (MQ=255)
                                                                                                                                            tGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCAta                                                                   >  2:922747/1‑88 (MQ=255)
                                                                                                                                            tGGATCAGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCg                                     >  3:417509/1‑118 (MQ=255)
                                                                                                                                             ggATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGt                                                                                                >  2:1033425/1‑58 (MQ=255)
                                                                                                                                                tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGcc                                                                              >  1:751566/1‑73 (MQ=255)
                                                                                                                                                tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTAc                                                                           <  4:358586/76‑1 (MQ=255)
                                                                                                                                                tCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCttcagttcag    >  4:1068810/1‑147 (MQ=255)
                                                                                                                                                 cTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGca                                                    >  3:1062579/1‑98 (MQ=255)
                                                                                                                                                 cTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCg  >  1:385923/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CACGCGCGCTGTCCCATTAATGGTGATGATGTGAATGGCCAGACACCCCTTCATTCAGGTGGCAGTCTGGCCCATGACAAGGTTGATATTCCATCTGAATGGTGGCGTGCTCAATCTGATAGTGATCCATCAGGTAGTGTTGGATCTGATCCAACAAGGCATCGTGATCGTGTGGCGGGATCACCTGCACATGCAACGTCATCACCGGCTTCTCGCCTACCATCCATACATGTACATGGTGCACATTGCGGACTTCCGGGATTTCCCGGCACATACGACGCTTCAGTTCAGCG  >  NC_000913/783866‑784158

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