Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,182,806 C→A 5.6% P86Q (CCA→CAA)  zupT → heavy metal divalent cation transporter ZupT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,182,8060CA 5.6% 236.9 / 2.8 71P86Q (CCA→CAA) zupTheavy metal divalent cation transporter ZupT
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base A (2/2);  total (32/39)
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.79e-01

CTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCT  >  NC_000913/3182658‑3182953
                                                                                                                                                    |                                                                                                                                                   
cTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTTTGTTGTGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGTCCTTGACCGCATGTTGCca                                                                                                                                                    <  2:449807/150‑1 (MQ=255)
 tAGGTTTTGCGGCGGGGATCATGTTGCTGATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcac                                                                                                                                                   <  2:288425/150‑1 (MQ=255)
  aGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTACTGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCcaca                                                                                                                                                  <  2:800694/150‑1 (MQ=255)
      tgtGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCt                                                                                                                                              <  2:1093132/148‑1 (MQ=255)
        tgtggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                                                                            <  1:880733/147‑4 (MQ=255)
        tGCGGCGGGGATCATGTTGCTCATCACATTAAGGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTGGCAACATGcggt                                                                                                                                            <  2:782507/150‑4 (MQ=255)
          cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCAta                                                                                                                                          >  1:1492598/1‑149 (MQ=255)
          cggcggGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATc                                                                                                                                          <  1:110543/150‑1 (MQ=255)
              ggggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCGACATGCTCATCCTCa                                                                                                                                      <  1:1089617/150‑1 (MQ=255)
               gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAg                                                                                                                                     <  2:270515/150‑1 (MQ=255)
               gggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGGCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAg                                                                                                                                     <  1:1062334/150‑1 (MQ=255)
                ggATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAgg                                                                                                                                    >  1:82559/1‑150 (MQ=255)
                    cATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGAttt                                                                                                                                >  2:463518/1‑150 (MQ=255)
                     aTGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTa                                                                                                                               <  1:1457763/150‑1 (MQ=255)
                       gTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                                                                             <  2:190949/150‑1 (MQ=255)
                       gTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAAt                                                                                                                             <  2:1425711/150‑1 (MQ=255)
                         tGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGc                                                                                                                           >  2:601054/1‑150 (MQ=255)
                           ctcatctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCaa                                                                                                                         >  1:621165/1‑150 (MQ=255)
                               tctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAAt                                                                                                                     >  1:798954/1‑150 (MQ=255)
                                ctcatTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATc                                                                                                                    <  1:639621/150‑1 (MQ=255)
                                     tAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGc                                                                                                               <  1:832460/150‑1 (MQ=255)
                                     tAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGc                                                                                                               <  2:1221262/150‑1 (MQ=255)
                                      aaTGGAAATGCTTCCTGCCGCACTACCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTACGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCa                                                                                                              >  1:1247319/1‑150 (MQ=255)
                                        tGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGc                                                                                                            >  2:509350/1‑150 (MQ=255)
                                                   ccTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCaaaa                                                                                                 >  1:1257286/1‑150 (MQ=255)
                                                     tGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATc                                                                                               >  1:552434/1‑150 (MQ=255)
                                                      gCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCg                                                                                              <  2:935394/150‑1 (MQ=255)
                                                         gCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGCCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                                                                           >  2:56564/1‑150 (MQ=255)
                                                         gCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                                                                           <  1:719828/150‑1 (MQ=255)
                                                         gCACTAGCAGCTGAAGGAATGTCGCCTGTGATGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATc                                                                                           <  2:954098/150‑1 (MQ=255)
                                                              agcagcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAgcg                                                                                      >  1:1134493/1‑150 (MQ=255)
                                                                 agcTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGcac                                                                                   >  1:974773/1‑150 (MQ=255)
                                                                   cTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAg                                                                                 <  1:463518/150‑1 (MQ=255)
                                                                    tGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGc                                                                                >  1:68870/1‑150 (MQ=255)
                                                                    tGAAGGAATGTCGCCGGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGc                                                                                <  2:941067/150‑1 (MQ=255)
                                                                     gAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGcc                                                                               >  2:689495/1‑150 (MQ=255)
                                                                             gTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCt                                                                       >  1:1288764/1‑150 (MQ=255)
                                                                             gTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCt                                                                       <  1:718136/150‑1 (MQ=255)
                                                                              tCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTc                                                                      <  1:509350/150‑1 (MQ=255)
                                                                                gCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCAc                                                                    <  2:798954/150‑1 (MQ=255)
                                                                                  cTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActc                                                                  <  2:621165/150‑1 (MQ=255)
                                                                                   tgtgTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCActct                                                                 <  1:190343/150‑1 (MQ=255)
                                                                                       ttGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGc                                                             <  1:1370172/150‑1 (MQ=255)
                                                                                        tGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCa                                                            <  1:152051/150‑1 (MQ=255)
                                                                                                ggAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTg                                                    <  1:912946/150‑1 (MQ=255)
                                                                                                 gAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGc                                                   >  1:239965/1‑150 (MQ=255)
                                                                                                   aTGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCAt                                                 <  2:452792/150‑1 (MQ=255)
                                                                                                     gTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATa                                                <  2:869230/149‑1 (MQ=255)
                                                                                                     gTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATa                                                >  1:869230/1‑149 (MQ=255)
                                                                                                      tttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAAc                                              <  2:1288764/150‑1 (MQ=255)
                                                                                                      tttATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAAc                                              <  1:601054/150‑1 (MQ=255)
                                                                                                           cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTccc                                         >  2:1391127/1‑150 (MQ=255)
                                                                                                           cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTccc                                         <  1:1267736/150‑1 (MQ=255)
                                                                                                           cTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTccc                                         >  1:312877/1‑150 (MQ=255)
                                                                                                            ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCATCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCg                                        <  2:408009/150‑1 (MQ=255)
                                                                                                            ttCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCg                                        >  1:1138359/1‑150 (MQ=255)
                                                                                                             tCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCgg                                       <  2:1257286/150‑1 (MQ=255)
                                                                                                             tCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCgg                                       <  2:974773/150‑1 (MQ=255)
                                                                                                                  cTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAggg                                  <  2:552434/150‑1 (MQ=255)
                                                                                                                  cTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCCGCATAACTTCCCGGAAggg                                  >  2:11564/1‑150 (MQ=255)
                                                                                                                   tGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGa                                 >  1:1367824/1‑150 (MQ=255)
                                                                                                                    gCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGAt                                >  1:1536560/1‑150 (MQ=255)
                                                                                                                    gCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGAt                                >  1:976094/1‑150 (MQ=255)
                                                                                                                    gCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGCCTGCATAACTTCCCGGAAGGGAt                                >  2:1003509/1‑150 (MQ=255)
                                                                                                                      ttGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTg                              <  1:279090/150‑1 (MQ=255)
                                                                                                                               ttGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTg                     <  2:51497/150‑1 (MQ=255)
                                                                                                                               ttGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTg                     <  2:395408/150‑1 (MQ=255)
                                                                                                                                  gCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCa                  >  1:575889/1‑150 (MQ=255)
                                                                                                                                     ggatgagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATAAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACgg               >  2:1223872/7‑150 (MQ=255)
                                                                                                                                        aCCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCCCTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCgc            >  1:1533178/1‑149 (MQ=255)
                                                                                                                                         gagCATGTGGCAACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGag           >  2:80056/3‑150 (MQ=255)
                                                                                                                                          cGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGagc          <  1:1143086/150‑1 (MQ=255)
                                                                                                                                                  gCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCt  >  1:94105/1‑150 (MQ=255)
                                                                                                                                                  gCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCt  >  1:1352692/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTAGGTTTTGCGGCGGGGATCATGTTGCTCATCTCATTAATGGAAATGCTTCCTGCCGCACTAGCAGCTGAAGGAATGTCGCCTGTGTTGGGTTATGGAATGTTTATCTTCGGTCTGCTTGGCTATTTTGGCCTGGACCGCATGTTGCCACATGCTCATCCGCAGGATTTAATGCAAAAATCGGTGCAGCCGTTGCCAAAATCGATCAAGCGCACAGCCATTCTGCTCACTCTCGGCATCAGTCTGCATAACTTCCCGGAAGGGATTGCCACCTTTGTCACGGCGAGCAGCAACCT  >  NC_000913/3182658‑3182953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑