Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,139,710 C→A 6.2% intergenic (+225/+49) yegH → / ← asmA inner membrane protein YegH/putative assembly protein AsmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,139,7100CA 6.2% 323.3 / 7.2 97intergenic (+225/+49)yegH/asmAinner membrane protein YegH/putative assembly protein AsmA
Reads supporting (aligned to +/- strand):  ref base C (45/46);  new base A (4/2);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.94e-01

AGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTG  >  NC_000913/2139561‑2139854
                                                                                                                                                     |                                                                                                                                                
aGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGc                                                                                                                                                  >  2:672767/1‑150 (MQ=255)
aGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGc                                                                                                                                                  <  1:438462/150‑1 (MQ=255)
 gCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCg                                                                                                                                                 >  1:41841/1‑150 (MQ=255)
    cGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACg                                                                                                                                              >  2:1085194/1‑150 (MQ=255)
      ggTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCt                                                                                                                                            >  2:1336774/1‑150 (MQ=255)
         tGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTa                                                                                                                                         <  1:1180916/150‑1 (MQ=255)
          gCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTAc                                                                                                                                        >  1:396006/1‑150 (MQ=255)
              gggTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGt                                                                                                                                    >  2:1063564/1‑150 (MQ=255)
               ggtggtGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTc                                                                                                                                   <  2:326926/150‑1 (MQ=255)
                gtggtgGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACATCCTTCTATCTTACATATTTAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCt                                                                                                                                  <  2:361566/150‑1 (MQ=255)
                 tggtggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGACCGCCTGATGCGACGCTTTACGCGTCtt                                                                                                                                 >  2:874452/1‑150 (MQ=255)
                   gtggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTAt                                                                                                                               >  2:1437813/1‑150 (MQ=255)
                    tggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATc                                                                                                                              <  2:761750/150‑1 (MQ=255)
                     ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCa                                                                                                                             >  1:661390/1‑150 (MQ=255)
                     ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCa                                                                                                                             <  2:521158/150‑1 (MQ=255)
                     ggAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCa                                                                                                                             <  1:1463921/150‑1 (MQ=255)
                       agaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACTCGTCTTATCatg                                                                                                                           <  2:468988/150‑3 (MQ=255)
                       agaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAgg                                                                                                                           <  2:1375052/150‑1 (MQ=255)
                        gaTTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGc                                                                                                                          >  2:642722/1‑150 (MQ=255)
                             aGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACa                                                                                                                     <  1:956059/150‑1 (MQ=255)
                              gCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGTCTACAg                                                                                                                    >  2:465779/1‑150 (MQ=255)
                              gCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAg                                                                                                                    <  2:19997/150‑1 (MQ=255)
                               ccACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAgg                                                                                                                   <  2:1280770/150‑1 (MQ=255)
                               ccACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATAAGACGCGTAAAGCGTCGCATCAGGcgttcgtg                                                                                                                   <  2:1015149/150‑8 (MQ=38)
                                  cccTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACt                                                                                                                >  2:1655022/1‑150 (MQ=255)
                                  cccTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACt                                                                                                                <  2:936550/150‑1 (MQ=255)
                                  cccTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACt                                                                                                                >  1:816106/1‑150 (MQ=255)
                                   ccTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTg                                                                                                               <  1:1437813/150‑1 (MQ=255)
                                    cTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTgg                                                                                                              <  1:401486/150‑1 (MQ=255)
                                     tGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGt                                                                                                             >  2:473590/1‑150 (MQ=255)
                                      ggCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTc                                                                                                            <  2:1410059/150‑1 (MQ=255)
                                       gCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCa                                                                                                           >  2:1482653/1‑150 (MQ=255)
                                            ttCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTgg                                                                                                      <  1:672767/150‑1 (MQ=255)
                                               aCCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCt                                                                                                                                  >  1:488515/1‑119 (MQ=255)
                                               aCCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCt                                                                                                                                  <  2:488515/119‑1 (MQ=255)
                                                 cGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAAcc                                                                                                 <  2:1071847/150‑1 (MQ=255)
                                                   gCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTa                                                                                               <  1:91481/150‑1 (MQ=255)
                                                    cTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTAc                                                                                              <  2:378365/150‑1 (MQ=255)
                                                      tGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAt                                                                                            >  2:437839/1‑150 (MQ=255)
                                                       gAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCTCCTGGTAACCTACAtc                                                                                           >  2:75127/1‑150 (MQ=255)
                                                         gCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACAtctt                                                                                         >  1:271085/1‑150 (MQ=255)
                                                                aTAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGc                                                                                  <  2:396006/150‑1 (MQ=255)
                                                                   gtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAc                                                                               >  2:339925/1‑150 (MQ=255)
                                                                   gtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAc                                                                               <  2:1564929/150‑1 (MQ=255)
                                                                   gtagCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAAc                                                                               >  1:300233/1‑150 (MQ=255)
                                                                       cGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGTTAACCTATATCTTCTCCAGGAActtct                                                                           <  2:816106/150‑1 (MQ=255)
                                                                          ggAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTca                                                                        <  1:272523/150‑1 (MQ=255)
                                                                          ggAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTca                                                                        >  1:1183495/1‑150 (MQ=255)
                                                                            aaCATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTcaca                                                                      <  2:1567256/150‑1 (MQ=255)
                                                                                tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCtt                                                                  >  1:362961/1‑150 (MQ=255)
                                                                                tCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCtt                                                                  <  1:962572/150‑1 (MQ=255)
                                                                                   aGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                               >  1:556185/1‑150 (MQ=255)
                                                                                   aGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                               >  1:1506223/1‑150 (MQ=255)
                                                                                   aGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGc                                                               <  1:912601/150‑1 (MQ=255)
                                                                                     tGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCa                                                             <  2:661390/150‑1 (MQ=255)
                                                                                     tGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCa                                                             <  1:104774/150‑1 (MQ=255)
                                                                                        aCGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTa                                                          <  2:316902/150‑1 (MQ=255)
                                                                                          ggAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACg                                                        <  1:642722/150‑1 (MQ=255)
                                                                                           gAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGACTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACgg                                                       >  2:789542/1‑150 (MQ=255)
                                                                                            aaCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACggg                                                      <  2:41841/150‑1 (MQ=255)
                                                                                            aaCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACggg                                                      <  2:913842/150‑1 (MQ=255)
                                                                                            aaCAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACggg                                                      >  1:764310/1‑150 (MQ=255)
                                                                                            aaCAGCCCTCCATCATACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACggg                                                      <  2:1482255/150‑1 (MQ=255)
                                                                                                  cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                <  1:1085194/150‑1 (MQ=255)
                                                                                                  cTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCtt                                                <  1:1063564/150‑1 (MQ=255)
                                                                                                     cGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTAt                                             >  2:1437539/1‑150 (MQ=255)
                                                                                                      gATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTAtt                                            <  1:473590/150‑1 (MQ=255)
                                                                                                           tACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCt                                       >  1:730735/1‑150 (MQ=255)
                                                                                                             cATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCg                                     >  1:1048203/1‑150 (MQ=255)
                                                                                                               tatTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCTGCTCGGc                                   <  1:1249528/150‑1 (MQ=255)
                                                                                                               tatTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGc                                   >  1:1555051/1‑150 (MQ=255)
                                                                                                               tatTCAGAATATAATGAGGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACCGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGc                                   <  1:1497904/150‑1 (MQ=255)
                                                                                                                atTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGcc                                  <  2:1183495/150‑1 (MQ=255)
                                                                                                                 ttCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGccc                                 >  2:1286341/1‑150 (MQ=255)
                                                                                                                  tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCa                                <  2:362961/150‑1 (MQ=255)
                                                                                                                  tCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCa                                >  2:1572123/1‑150 (MQ=255)
                                                                                                                    aGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGt                              <  1:1482653/150‑1 (MQ=255)
                                                                                                                    aGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTATACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGt                              >  1:831160/1‑150 (MQ=255)
                                                                                                                         ataATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTc                         <  2:979645/150‑1 (MQ=255)
                                                                                                                          taATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCa                        >  2:973585/1‑150 (MQ=255)
                                                                                                                          taATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCa                        >  2:424911/1‑150 (MQ=255)
                                                                                                                           aaTGATGCCGCCACGAACGACTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTCATTCCGCTCGGCCCAGTCATTCAg                       >  2:582927/1‑150 (MQ=255)
                                                                                                                            atgatgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCTCCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAgg                      <  1:75127/150‑1 (MQ=255)
                                                                                                                                tgCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGAc                  >  2:920307/1‑150 (MQ=255)
                                                                                                                                 gccgccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCCCCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTATGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACg                 >  1:893469/1‑150 (MQ=255)
                                                                                                                                 gccgccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACg                 >  2:23210/1‑150 (MQ=255)
                                                                                                                                  ccgccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACCTCTTCTCCAGCACCTTCTTCCCATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGt                >  1:1467063/1‑150 (MQ=255)
                                                                                                                                  ccgccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGt                >  2:185207/1‑150 (MQ=255)
                                                                                                                                    gccACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGttt              <  2:1048203/150‑1 (MQ=255)
                                                                                                                                        cGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCa                            <  2:627998/132‑1 (MQ=255)
                                                                                                                                        cGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCa                            >  1:627998/1‑132 (MQ=255)
                                                                                                                                           aggCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCg       >  1:145474/3‑150 (MQ=255)
                                                                                                                                            ggCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGt      >  2:1455329/2‑150 (MQ=255)
                                                                                                                                            cGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGt      >  1:20822/1‑150 (MQ=255)
                                                                                                                                             gCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTc     >  2:821588/1‑150 (MQ=255)
                                                                                                                                             gCCTGATAAGACGCGTAAAGCGTCGCATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTc     >  2:957812/1‑150 (MQ=255)
                                                                                                                                                tGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCACCTTCTTCACATCTTTGCCATTACGGGAATCTTTGTTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTg  >  2:1015634/1‑150 (MQ=255)
                                                                                                                                                    gCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCtt         >  1:990348/1‑139 (MQ=255)
                                                                                                                                                    gCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCtt         <  2:990348/139‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCCCGGGTTGCCAGGGTGGTGGAGATTGAGCCACCCTGGCACGTTCACCGGCTTGAGCGTTCCATAGTAGCGAGGAACATCAAGTGAACGGAACAGCCTTCGATCTTACATATTCAGAATATAATGATGCCGCCACGAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACAGGACTGGTCACCTGGTAACCTACATCTTCTCCAGCAACTTCTTCACATCTTTGCCATTACGGGAATCTTTATTCCGCTCGGCCCAGTCATTCAGGCGACGTTTCGCTTCGTCCTG  >  NC_000913/2139561‑2139854

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: