Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,101,656 T→A 9.4% intergenic (‑33/+151) yggN ← / ← yggL DUF2884 family putative periplasmic protein/DUF469 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,101,6560TA 9.4% 321.9 / 17.1 125intergenic (‑33/+151)yggN/yggLDUF2884 family putative periplasmic protein/DUF469 family protein
Reads supporting (aligned to +/- strand):  ref base T (59/54);  new base A (5/7);  total (64/61)
Fisher's exact test for biased strand distribution p-value = 5.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 125 total aligned reads displayed.

CGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGCT  >  NC_000913/3101509‑3101803
                                                                                                                                                   |                                                                                                                                                   
cGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTg                                                                                                                                                   >  2:397877/1‑150 (MQ=255)
 ggTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                                                                                                  >  1:950242/1‑150 (MQ=255)
  gTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCt                                                                                                                                                 <  2:1357681/150‑1 (MQ=255)
    ttGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                                                                                                               <  2:1554405/150‑1 (MQ=255)
    ttGCGGGCTGACAATCACAGCCTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                                                                                                               <  2:536065/150‑1 (MQ=255)
      gCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcc                                                                                                                                             <  1:28710/150‑1 (MQ=255)
        gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                                                                                                           >  1:490545/1‑150 (MQ=255)
        gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                                                                                                                                           <  2:1246061/150‑3 (MQ=255)
        gggCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctg                                                                                                                                           <  1:898543/150‑3 (MQ=255)
           cTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAAGCAGGCCAGCAAATGCTGGcctgatt                                                                                                                                        <  1:83465/150‑6 (MQ=255)
            tGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgattt                                                                                                                                       <  2:1004888/150‑7 (MQ=255)
               cAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                                                                                                    >  1:713600/1‑150 (MQ=255)
                 aTCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                                                                                                                  >  2:1435929/1‑150 (MQ=255)
                    gcagcctCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTAATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGGATCATAGTCTTCCCTCAAGAAAAAATCAGGCCACCATTTGCTGGCACCAGATTGATAACa                                                                                                                               >  2:170130/7‑150 (MQ=255)
                    gcagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATATTGCCCATCAAATAATTCCCTCAAGAAAAAAACAGGCCCGCATTTGATGGCCCCAGATGGATAACa                                                                                                                               >  2:866958/5‑150 (MQ=255)
                     cagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAATAATCAGGCCAGCATTTGCTGGCCCCAGATTGATATCaa                                                                                                                              >  1:1492540/4‑150 (MQ=255)
                     cagcGTCACGCCGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa                                                                                                                              >  1:1485709/4‑150 (MQ=255)
                      ctCGGCACGCGGCCTGACGCTGCCCTGGTAGCCGGCGTGAGCGGGCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGCCCAGCATTTGCTGGCCCCAGATTGATAACaac                                                                                                                             <  2:862389/149‑2 (MQ=255)
                       tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                            <  1:422834/150‑1 (MQ=255)
                       tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                            <  1:1362260/150‑1 (MQ=255)
                       tCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                                                                                                            <  1:1162341/150‑1 (MQ=255)
                        cGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGt                                                                                                                           <  2:1576424/150‑1 (MQ=255)
                         gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                                                                                                          >  2:246862/1‑150 (MQ=255)
                         gTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                                                                                                          <  1:330534/150‑1 (MQ=255)
                           cACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                                                                                                        <  2:691527/150‑1 (MQ=255)
                           cACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTACCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                                                                                                        <  1:80494/150‑1 (MQ=255)
                            aCGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATAGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc                                                                                                                       <  1:1095158/150‑1 (MQ=255)
                             cgcgGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgcg                                                                                                                      >  1:611649/1‑150 (MQ=255)
                              gcggcgTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGACCCAGATTGATAACAAAGTACGCGt                                                                                                                     >  1:457079/1‑150 (MQ=255)
                                     gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAAATGCTGGcctgat                                                                                                                                         >  2:1290397/1‑119 (MQ=255)
                                       cGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCACAGATTGTTAACAAAGTACGCGTTTTCCATGc                                                                                                            >  2:345250/1‑150 (MQ=255)
                                         cTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg                                                                                                          >  1:1711446/1‑150 (MQ=255)
                                            cACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt                                                                                                       >  2:383070/1‑150 (MQ=255)
                                             aCTGGTAGTCGGCGTGAGAGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATCACAAAGTCCGCGTTTTCCATGCCGGATg                                                                                                      <  2:1073385/150‑1 (MQ=255)
                                              cTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc                                                                                                     >  2:1570238/1‑150 (MQ=255)
                                              cTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgc                                                                                                     >  2:1729200/1‑150 (MQ=255)
                                                 gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc                                                                                                  <  2:478331/150‑1 (MQ=255)
                                                 gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc                                                                                                  <  1:1362239/150‑1 (MQ=255)
                                                 gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggc                                                                                                  <  2:305777/150‑1 (MQ=255)
                                                 gTAGTCGGAGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAtgctgc                                                                                                  <  1:1575479/150‑1 (MQ=255)
                                                  tAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcg                                                                                                 <  1:96446/150‑1 (MQ=255)
                                                    gTCGGCGGGAGCGGGCATTGCCGTCACGGATAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCACCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTa                                                                                               <  2:71373/150‑1 (MQ=255)
                                                      cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCCGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGAGTTTTCCATGCCGGATGCGGCGTaac                                                                                             >  2:1053229/1‑149 (MQ=255)
                                                      cGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaaa                                                                                             >  2:1635409/1‑150 (MQ=255)
                                                       ggCGTGAGCGGGCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCGGGCCCCAGATTGATAACAAAGTACGCGTTATCCATGCCGGATGCGGCGTAAac                                                                                            <  2:949323/150‑1 (MQ=255)
                                                         cGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac                                                                                          >  2:883677/1‑150 (MQ=255)
                                                           tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt                                                                                        <  2:269196/150‑1 (MQ=255)
                                                           tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt                                                                                        <  2:452739/150‑1 (MQ=255)
                                                             aGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCAATTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCAATGCCGGATGCGGCGTAAACACCTTa                                                                                      >  2:115963/1‑150 (MQ=255)
                                                               cGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATc                                                                                    >  1:1642921/1‑150 (MQ=255)
                                                                  tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCg                                                                                 <  1:246862/150‑1 (MQ=255)
                                                                    aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGcc                                                                               <  1:25962/150‑1 (MQ=255)
                                                                       gCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTCGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTAc                                                                            >  1:25406/1‑150 (MQ=255)
                                                                        ccGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACt                                                                           <  2:688911/150‑1 (MQ=255)
                                                                         cGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGAGTAAACACCTTATCCCGCCTAcac                                                                          >  2:1030967/1‑148 (MQ=255)
                                                                            cACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGAt                                                                       >  1:146400/1‑150 (MQ=255)
                                                                             aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCAATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCgctg                                                                      >  1:1443862/1‑147 (MQ=255)
                                                                              cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCg                                                                     >  1:608629/1‑150 (MQ=255)
                                                                                gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                   <  2:1436746/150‑1 (MQ=255)
                                                                                 aaaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTATTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGc                                                                  <  2:807846/150‑1 (MQ=255)
                                                                                  aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTACATGCCGGCTGCGGCGTTAACACCGTATCCCGCCTACTCGATCGTGCa                                                                 >  2:96019/1‑150 (MQ=255)
                                                                                         gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTcc                                                          >  2:203948/1‑150 (MQ=255)
                                                                                             ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGc                                                      >  2:1437665/1‑150 (MQ=255)
                                                                                                 cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGa                                                  >  2:34953/1‑150 (MQ=255)
                                                                                                  agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGAt                                                 <  2:1351830/150‑1 (MQ=255)
                                                                                                  agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGAt                                                 >  2:1484820/1‑150 (MQ=255)
                                                                                                           cgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAgc                                        >  2:736104/1‑150 (MQ=255)
                                                                                                             catcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACACGATCGTGCATCTTTCCAGGCCTGATAAGCGTAgcgc                                      >  2:455839/1‑150 (MQ=255)
                                                                                                              atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATACCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCa                                     >  2:757803/1‑150 (MQ=255)
                                                                                                                 atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCa                                  <  1:1212579/150‑1 (MQ=255)
                                                                                                                     tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCa                              <  1:1579390/150‑1 (MQ=255)
                                                                                                                         cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGttt                          >  1:1365052/1‑150 (MQ=255)
                                                                                                                           cTCAAGAAATAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTg                        <  1:314390/150‑1 (MQ=255)
                                                                                                                           cTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTg                        >  1:1642691/1‑150 (MQ=255)
                                                                                                                            tCAAGAAAAAATCAGGCCAGCATTTGCTGGCCGCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTACACCCCTTATCCCGCCGACTCGATCGTGCATCGTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGGTTTgg                       >  1:147840/1‑149 (MQ=255)
                                                                                                                             cAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCa                      <  2:1310483/150‑1 (MQ=255)
                                                                                                                               aGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAtt                    >  2:489887/1‑150 (MQ=255)
                                                                                                                                gAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCAtt                    >  2:191307/1‑149 (MQ=255)
                                                                                                                                 aaaaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCAACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTg                  <  1:275973/150‑1 (MQ=255)
                                                                                                                                   aaaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCACATCAGGCAGTTTTGCATTTGTc                <  2:650808/150‑1 (MQ=255)
                                                                                                                                    aaaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCa               <  1:397877/150‑1 (MQ=255)
                                                                                                                                     aaTCAGGCCAGCATTTGCTGGCCCCAGATTGATATCAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCATGCCTGATAAGCGTTGCGCATCAGGCAGTTTTGCATTTGTCAt              <  2:1351350/150‑1 (MQ=255)
                                                                                                                                     aaTCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCAt              >  1:115936/1‑150 (MQ=255)
                                                                                                                                        ggggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGc           >  2:1276874/3‑150 (MQ=255)
                                                                                                                                        ggggCCAGCAAATACTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTTCATTTGTAATCGc           >  2:483330/3‑150 (MQ=39)
                                                                                                                                         gggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGcc          >  1:672246/2‑150 (MQ=255)
                                                                                                                                         gggCCAGCAAATGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGcc          >  1:1450657/2‑150 (MQ=255)
                                                                                                                                         aGGCCAGCATTTGCTGGCCCCACATTGATAACAAAGTACGTGTTTTCCAAGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGGGCCTCTTTCCAGGCCTGAGAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGcc          <  2:893666/150‑1 (MQ=39)
                                                                                                                                          ggCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTACACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCCGGCCTGATAATCGTATCGCATCAGGCAGTTTTGCATTTGTCTTCTccc         >  2:1197457/1‑150 (MQ=255)
                                                                                                                                          ggCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGccc         >  1:567854/1‑150 (MQ=255)
                                                                                                                                          ggCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGccc         >  1:582483/1‑150 (MQ=255)
                                                                                                                                           gCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATGCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCTTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCg        <  1:1298264/150‑1 (MQ=255)
                                                                                                                                            ccAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt       <  2:626029/150‑1 (MQ=255)
                                                                                                                                            ccAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGt       <  2:950242/150‑1 (MQ=255)
                                                                                                                                             cAGCATTTGCTGGCCCCGGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                   <  2:1641064/89‑1 (MQ=255)
                                                                                                                                             cAGCATTTGCTGGCCCCGGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                   >  1:1641064/1‑89 (MQ=255)
                                                                                                                                             cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                   >  1:938120/1‑89 (MQ=255)
                                                                                                                                             cAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTg                                                                   <  2:938120/89‑1 (MQ=255)
                                                                                                                                               gCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATg    >  1:31873/1‑150 (MQ=255)
                                                                                                                                                 aTTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCAATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGCt  >  2:1640995/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGTTTGCGGGCTGACAATCACATCGTCACGCGGCGTGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGCCTACTCGATCGTGCATCTTTCCAGGCCTGATAAGCGTAGCGCATCAGGCAGTTTTGCATTTGTCATCGCCCGTATGCT  >  NC_000913/3101509‑3101803

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