Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,074,814:1 +A 5.2% coding (36/1533 nt) yeeR → CP4‑44 prophage; inner membrane protein YeeR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,8141.A 5.2% 470.1 / 19.9 116coding (36/1533 nt)yeeRCP4‑44 prophage; inner membrane protein YeeR
Reads supporting (aligned to +/- strand):  ref base . (46/64);  new base A (2/4);  total (48/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 118 total aligned reads displayed.

ATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCG  >  NC_000913/2074666‑2074960
                                                                                                                                                      |                                                                                                                                                  
aTCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCg                                                                                                                                                   <  1:1472048/150‑1 (MQ=255)
  cGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCa                                                                                                                                                 <  1:1241010/150‑1 (MQ=255)
  cGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCa                                                                                                                                                 <  1:1246885/150‑1 (MQ=255)
   gCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAg                                                                                                                                                <  1:642906/150‑1 (MQ=255)
        tctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAttt                                                                                                                                           >  2:22640/1‑150 (MQ=255)
         ctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg                                                                                                                                          <  2:305402/150‑1 (MQ=255)
         ctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg                                                                                                                                          <  2:1668/150‑1 (MQ=255)
         ctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAtcag                                                                                                                                          <  1:446945/150‑4 (MQ=255)
         ctGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAtcag                                                                                                                                          <  2:622784/150‑4 (MQ=255)
          tgtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAtcagc                                                                                                                                         <  2:1482546/150‑5 (MQ=255)
           gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc                                                                                                                                        <  1:969640/150‑1 (MQ=255)
           gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc                                                                                                                                        <  1:963490/150‑1 (MQ=255)
           gtgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCGATCAGCAGGAtcagcg                                                                                                                                        <  2:1469570/150‑6 (MQ=255)
            tgGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCa                                                                                                                                       <  1:321897/150‑1 (MQ=255)
             ggTCCCGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAt                                                                                                                                      >  1:1352962/1‑150 (MQ=255)
             ggTCCCGGTAATAATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCAAATTCTCTGGCATTAACAATAACATGATATTCATCACGGATTGACTATTTTAAAGATAGTCTGCGCGCTGATCCTGCTTATCGCAGGATTTGCCAt                                                                                                                                      >  2:1428381/1‑150 (MQ=255)
              gTCCCGGTCATGATGACCGGGACCACAGAGAGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAAATGATATTCATCACGGAGTGAGTATGTTACAGATATTCGGCGCGATGAACATGCTGATCCCAGGATTTGCCAt                                                                                                                                      >  2:543093/1‑150 (MQ=255)
              gTCCCGGTCATGATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAtt                                                                                                                                     >  2:819976/1‑150 (MQ=255)
                cccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                                                   >  2:1020830/1‑150 (MQ=255)
                cccGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAttct                                                                                                                                   >  1:1623256/1‑150 (MQ=255)
                   ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCg                                                                                                                                <  2:1447923/150‑1 (MQ=255)
                   ggTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCg                                                                                                                                <  1:1263256/150‑1 (MQ=255)
                      catcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGAGTATGTTACAGATAGTCGGCGCGCTGAACCTGCTGATCGCAGGATTTGCCATTCTTCGCCt                                                                                                                             >  1:1506019/1‑150 (MQ=255)
                        tcatGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCttt                                                                                                                           <  1:413423/150‑1 (MQ=255)
                          atGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTg                                                                                                                         <  2:1352962/150‑1 (MQ=255)
                             aCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGAAAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTc                                                                                                                      >  2:881744/1‑150 (MQ=255)
                               cGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                                                    <  2:1066416/150‑1 (MQ=255)
                               cGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                                                    <  2:1457581/150‑1 (MQ=255)
                                gggACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCaga                                                                                                                   <  2:606144/150‑1 (MQ=255)
                                    cccGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAt                                                                                                               >  1:318085/1‑150 (MQ=255)
                                    cccGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAt                                                                                                               >  1:624340/1‑150 (MQ=255)
                                      ccttaCCGTCGAAACGGATCTTCAACTCCACAATAGATGGAATTAACAATAACATTATAATAATCACAGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCAtt                                                                                                              >  2:542320/5‑150 (MQ=255)
                                        ggACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTaa                                                                                                            <  1:1547209/150‑1 (MQ=255)
                                           ccGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCa                                                                                                         <  2:18750/150‑1 (MQ=255)
                                            cGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAg                                                                                                        <  1:703913/150‑1 (MQ=255)
                                                  aaCGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCag                                                                                                                    >  1:1438504/1‑132 (MQ=255)
                                                  aaCGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGc                                                                                                  >  1:572678/1‑150 (MQ=255)
                                                          ttCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTgcgc                                                                                          >  1:1180/1‑150 (MQ=255)
                                                           tCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCt                                                                                         <  1:1436986/150‑1 (MQ=255)
                                                            cAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTg                                                                                        <  1:786311/150‑1 (MQ=255)
                                                            cAACGCCACATTCGCTGCCATTAACAATAACATGATATTCATCACGGAGCGACTATGTTACAGCTAGTCGGCGCGCTGATCCTGCTCATCGCAGGATTTGCCATTCTTCGCCTGTTGTTCAGAGCATAAATCAGCAAGGCTTCTGCGCTg                                                                                        <  2:1157274/150‑1 (MQ=255)
                                                              aCGCCACATTCGCTGGCATTAACTATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCTGGATTTGGCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGc                                                                                      <  1:837227/150‑1 (MQ=255)
                                                               cGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCa                                                                                     <  1:1192464/150‑1 (MQ=255)
                                                                 ccACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAgg                                                                                   <  2:587109/150‑1 (MQ=255)
                                                                  cacaTTCGCTGGCATTAAGAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTTGCAggg                                                                                  <  1:200649/150‑1 (MQ=255)
                                                                   acaTTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGc                                                                                 <  2:1383824/150‑1 (MQ=255)
                                                                      ttCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCa                                                                              <  2:1180/150‑1 (MQ=255)
                                                                            ggCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                        >  1:985542/1‑150 (MQ=255)
                                                                            ggCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGc                                                                        <  2:624340/150‑1 (MQ=255)
                                                                              cATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtg                                                                      >  1:202488/1‑150 (MQ=255)
                                                                                ttAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCtgtg                                                                    <  1:819976/150‑1 (MQ=255)
                                                                                  aaCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCATGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                                                  <  2:210679/150‑1 (MQ=255)
                                                                                  aaCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCGTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                                                  >  1:721005/1‑150 (MQ=255)
                                                                                  aaCAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTc                                                                  <  2:572678/150‑1 (MQ=255)
                                                                                    cAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTg                                                                <  1:543093/150‑1 (MQ=255)
                                                                                       tatcttGATATTCATCACGGTGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCATGATTTGCCATTGTTCGCCTTTTGTTCAGAGCATTAATCAGCACTGCTTCTGCGATGGCAGGGCTGATATTGCTGTGTCTGTTc                                                             <  2:1381584/145‑1 (MQ=255)
                                                                                       tAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTc                                                             >  2:535838/1‑150 (MQ=255)
                                                                                        aaCATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcg                                                            <  1:1463754/150‑1 (MQ=255)
                                                                                          cATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggc                                                          <  2:440425/150‑1 (MQ=255)
                                                                                           aTGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcc                                                         <  1:1020830/150‑1 (MQ=255)
                                                                                                atTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccggcc                                                    <  2:274860/150‑1 (MQ=255)
                                                                                                 ttCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCt                                                   <  2:96066/150‑1 (MQ=255)
                                                                                                  tcatcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCtt                                                  >  2:1155026/1‑150 (MQ=255)
                                                                                                   catcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTa                                                 >  1:351178/1‑150 (MQ=255)
                                                                                                    atcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTAc                                                >  2:1473014/1‑150 (MQ=255)
                                                                                                     tcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGc                                                                                                                 >  1:164119/1‑84 (MQ=255)
                                                                                                     tcaCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGc                                                                                                                 <  2:164119/84‑1 (MQ=255)
                                                                                                      caCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTcggcccggc                                                     <  2:1196085/143‑1 (MQ=255)
                                                                                                      caCGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActg                                              >  1:694018/1‑150 (MQ=255)
                                                                                                        cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggc                                            <  2:223268/150‑1 (MQ=255)
                                                                                                        cGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActggc                                            <  1:632233/150‑1 (MQ=255)
                                                                                                         ggAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCt                                                                                              <  2:136167/99‑1 (MQ=255)
                                                                                                             tGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTggctggct                                       >  2:684832/1‑150 (MQ=255)
                                                                                                              gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                      <  1:22640/150‑1 (MQ=255)
                                                                                                              gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                      <  1:1068213/150‑1 (MQ=255)
                                                                                                              gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                      <  1:18316/150‑1 (MQ=255)
                                                                                                              gACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCta                                      <  1:929834/150‑1 (MQ=255)
                                                                                                                cTATGTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTatat                                   >  1:981148/1‑150 (MQ=255)
                                                                                                                    gTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGTCTATATCAc                                >  1:1494208/1‑150 (MQ=255)
                                                                                                                     ttACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCAcc                               <  1:439941/150‑1 (MQ=255)
                                                                                                                     ttACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCAcc                               >  1:110518/1‑150 (MQ=255)
                                                                                                                     ttACAGAAAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCAcc                               <  1:881744/150‑1 (MQ=255)
                                                                                                                      tACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGAATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTActg                                              <  2:642370/134‑1 (MQ=255)
                                                                                                                          gATAGTCGGGGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACg                          >  1:1449133/1‑150 (MQ=255)
                                                                                                                          gATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACg                          >  2:548593/1‑150 (MQ=255)
                                                                                                                          gATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACg                          <  2:930466/150‑1 (MQ=255)
                                                                                                                            tAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCa                        <  1:1220559/150‑1 (MQ=255)
                                                                                                                            tAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCa                        >  2:1473521/1‑150 (MQ=255)
                                                                                                                            tAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCa                        <  2:1163888/150‑1 (MQ=255)
                                                                                                                              gTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATa                      >  1:689782/1‑150 (MQ=255)
                                                                                                                                cGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAAc                    >  1:1259437/1‑150 (MQ=255)
                                                                                                                                 ggCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAAcc                   >  1:60266/1‑150 (MQ=255)
                                                                                                                                     ggcaaaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgt               >  2:1146560/6‑150 (MQ=255)
                                                                                                                                      gcaaaTCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgtt              >  1:15836/5‑150 (MQ=255)
                                                                                                                                       caaaTCCTGCGATCAGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttg             >  1:647855/4‑150 (MQ=255)
                                                                                                                                       cTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGgttg             <  1:1482154/150‑1 (MQ=255)
                                                                                                                                            ccTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCa        >  2:1609938/1‑150 (MQ=255)
                                                                                                                                              tgctgATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATACCCCGGTTGTTCCata      >  2:652646/1‑150 (MQ=255)
                                                                                                                                                ctgATCGCAGGATTTGCCATGCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATAtt    <  1:1004985/150‑1 (MQ=255)
                                                                                                                                                  gATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCg  <  2:202488/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTTGTTCAGAGCATTAATCAGCACGGCTTCTGCGCTGGCAGGGCTCATATTGCTGTGTCTGTTCGGCCCGGCCTTACTGGCTGGCTATATCACCGAACGCATAACCCGGTTGTTCCATATTCG  >  NC_000913/2074666‑2074960

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: