Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 111,564 T→G 5.4% intergenic (+131/+85) mutT → / ← yacG dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase/DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,5640TG 5.4% 321.4 / 5.6 109intergenic (+131/+85)mutT/yacGdGTP‑preferring nucleoside triphosphate pyrophosphohydrolase/DNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base T (52/51);  new base G (4/2);  total (56/53)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 111 total aligned reads displayed.

GAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATT  >  NC_000913/111415‑111709
                                                                                                                                                     |                                                                                                                                                 
gAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGtt                                                                                                                                                   <  1:1007268/150‑1 (MQ=255)
gAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGtt                                                                                                                                                   <  2:1019084/150‑1 (MQ=255)
 aaGCCTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTc                                                                                                                                                  <  1:17770/150‑1 (MQ=255)
  aGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCa                                                                                                                                                 >  1:197653/1‑150 (MQ=255)
  aGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCa                                                                                                                                                 <  1:684766/150‑1 (MQ=255)
    cTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAcg                                                                                                                                               >  2:689737/1‑150 (MQ=255)
       aaaCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAcgccg                                                                                                                                            <  1:1606709/150‑1 (MQ=255)
        aaCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAcgccgc                                                                                                                                           >  2:1399259/1‑150 (MQ=255)
          cGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                                                                                                         >  2:165276/1‑150 (MQ=255)
          cGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                                                                                                         <  1:33469/150‑1 (MQ=255)
            tCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                                                                                                                       <  1:736257/150‑1 (MQ=255)
             cTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                      >  2:1130723/1‑150 (MQ=255)
               gTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                                                                                                                                    <  1:365255/150‑1 (MQ=255)
                 aGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTTAAGGGATTTCTAACTCATTGATAAATTTTTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAt                                                                                                                                  <  2:154159/150‑1 (MQ=255)
                 aGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAcct                                                                                                                                  <  2:14977/150‑3 (MQ=37)
                  ggTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGTGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt                                                                                                                                 <  2:612376/150‑1 (MQ=255)
                  ggTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt                                                                                                                                 >  2:1322897/1‑150 (MQ=255)
                   gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAttt                                                                                                                                >  2:1037966/1‑150 (MQ=255)
                   gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAttt                                                                                                                                >  2:839117/1‑150 (MQ=255)
                   gTCAGATAAGGCGTTTTCGCCGCATCCGACATTAGCACACGCTGCCTGATGCGACTCTGGCGCGGCTTAGCAGTCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTGTGTAGGTCGGATAAGGCGTTAACGCCGCATCCGACAttt                                                                                                                                >  1:733521/1‑150 (MQ=255)
                    tCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTg                                                                                                                               <  1:1046492/150‑1 (MQ=255)
                       gATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGcac                                                                                                                            <  1:1011953/150‑1 (MQ=255)
                         aaaGGCGGTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACaa                                                                                                                          <  1:1432863/149‑1 (MQ=255)
                          aaGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAg                                                                                                                         <  1:42823/150‑1 (MQ=255)
                           aGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGa                                                                                                                        <  1:998057/150‑1 (MQ=255)
                              cGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGc                                                                                                                     >  2:1050590/1‑150 (MQ=255)
                               gTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATATGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGcc                                                                                                                    <  1:299432/150‑1 (MQ=255)
                                 tttCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTAATCAGGCCTAAAGGGATTTCTAACTCAATGATAAATTTGTTTTTATAGGTCGGATAAGGCGTTCACGCCCCATCCCACAATTGCACAAGATGCCTg                                                                                                                  >  1:131148/1‑150 (MQ=39)
                                    cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATg                                                                                                               <  1:1202008/150‑1 (MQ=255)
                                        gcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGAc                                                                                                           <  2:1031029/150‑1 (MQ=255)
                                         cATCCGACACTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACg                                                                                                          <  2:549599/150‑1 (MQ=255)
                                              gACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTc                                                                                                     <  1:1050590/150‑1 (MQ=255)
                                               aCATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTcc                                                                                                    >  1:1328951/1‑150 (MQ=255)
                                                cATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgc                                                                                                  >  1:396048/1‑150 (MQ=255)
                                                    cGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTc                                                                                               >  2:1657512/1‑150 (MQ=255)
                                                    cGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTc                                                                                               <  2:155447/150‑1 (MQ=255)
                                                     gCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCt                                                                                              >  2:137341/1‑150 (MQ=255)
                                                      cacaCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTTGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCtt                                                                                             >  2:748446/1‑150 (MQ=255)
                                                      cacaCGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGTCGTTCACGCCTCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCtt                                                                                             >  2:266182/1‑150 (MQ=255)
                                                       acacGATGCCTGATGCGACGCTGGCGCGTCTTATGAGGCCTAAAGTGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTa                                                                                            <  2:511097/150‑1 (MQ=255)
                                                          cGATGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATAAGACGCGGACAGCGTCGCATCa                                                                                         <  2:367569/150‑1 (MQ=38)
                                                           gATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATGTTTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCTTCCGACATTTGCACAAGATGCCTGATGCGACGCGGACAGCGTCGCATCAg                                                                                        <  2:898136/150‑1 (MQ=38)
                                                              gCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGAAAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGcc                                                                                     <  2:1498533/150‑1 (MQ=255)
                                                                  gATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACg                                                                                 <  2:2834/150‑1 (MQ=255)
                                                                        aCGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGc                                                                           <  1:1399259/150‑1 (MQ=255)
                                                                        aCGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACACATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATACGACATTTGCACAAGATGCCAGATGCGACGCTGTCAGCGTCTTTTCAGGCCTACGTGCGGc                                                                           >  2:684443/1‑150 (MQ=39)
                                                                         cGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCa                                                                          <  2:1623708/150‑1 (MQ=255)
                                                                         cGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCa                                                                          <  1:1638065/150‑1 (MQ=255)
                                                                          gCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCAt                                                                         <  1:1581446/150‑1 (MQ=255)
                                                                          gCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCAt                                                                         >  1:1440697/1‑150 (MQ=255)
                                                                           cTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATc                                                                        <  1:165276/150‑1 (MQ=255)
                                                                            tGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCa                                                                       >  1:465012/1‑150 (MQ=255)
                                                                                  gTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACaaa                                                                 <  1:1035751/150‑1 (MQ=255)
                                                                                   tCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTTCACAAGATGCCTGATGCTACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAAt                                                                >  2:825564/1‑150 (MQ=255)
                                                                                      tATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTc                                                             >  1:1062629/1‑150 (MQ=255)
                                                                                      tATCAGGCCTAAAGGGATTACTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATTCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAAAGTc                                                             >  1:805438/1‑150 (MQ=255)
                                                                                       aTCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTATGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAATATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGAGCCGCATCAGACAAAAGTCa                                                            >  1:323119/1‑150 (MQ=38)
                                                                                         cAGGCCTAAAGGGATATCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACt                                                          <  2:459602/150‑1 (MQ=255)
                                                                                            gCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCt                                                       >  2:1494892/1‑150 (MQ=255)
                                                                                              cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGCTGCATGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                     >  2:106584/1‑150 (MQ=255)
                                                                                              cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                     >  2:1378358/1‑150 (MQ=255)
                                                                                              cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                                     <  1:1130723/150‑1 (MQ=255)
                                                                                                aaaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTgg                                                   >  2:848875/1‑150 (MQ=255)
                                                                                                aaaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGACTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTgg                                                   >  2:1721945/1‑150 (MQ=255)
                                                                                                 aaggtctttCTCACTCATTGATAAATTTGGTTTTGTAGGTCGGATAAGGCGTGCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                                  <  1:328655/144‑1 (MQ=38)
                                                                                                 aaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                                  >  1:1696103/1‑150 (MQ=255)
                                                                                                 aaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                                  <  1:784367/150‑1 (MQ=255)
                                                                                                   gggATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttc                                                >  2:1647217/1‑150 (MQ=255)
                                                                                                       tttCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGAAGCTGTCCGCGTCTTATCAGGGCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCg                                            >  1:1083147/1‑150 (MQ=255)
                                                                                                           tAACTCTTTGATAAATTTGTTTTTGTAGGTCTGATAAGGCTTTCACGCCGCATCCGAAATTTGGACAAGATGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                        <  2:1308334/150‑1 (MQ=38)
                                                                                                           tAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                        >  2:573548/1‑150 (MQ=255)
                                                                                                           tAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                        <  2:1442117/150‑1 (MQ=255)
                                                                                                           tAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                                        <  1:1281041/150‑1 (MQ=255)
                                                                                                            aaCTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCa                                       >  2:481116/1‑150 (MQ=255)
                                                                                                             aCTCATTTATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAg                                      <  1:800853/150‑1 (MQ=255)
                                                                                                             aCTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAg                                      <  2:1505391/150‑1 (MQ=255)
                                                                                                                  ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCAtt                                 <  1:266182/150‑2 (MQ=255)
                                                                                                                  ttGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATc                                 >  2:1670370/1‑150 (MQ=255)
                                                                                                                   tGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCg                                >  1:1217920/1‑150 (MQ=255)
                                                                                                                   tGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCAGCg                                >  2:166820/1‑150 (MQ=255)
                                                                                                                    gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGc                               >  1:1389429/1‑150 (MQ=255)
                                                                                                                    gATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGc                               >  2:830925/1‑150 (MQ=255)
                                                                                                                     aTAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCt                              >  2:177356/1‑150 (MQ=255)
                                                                                                                         aTTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCg                          <  2:1312177/150‑1 (MQ=255)
                                                                                                                             gTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaaa                      <  2:124614/150‑1 (MQ=255)
                                                                                                                             gTTTCTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaaa                      <  2:663417/150‑1 (MQ=255)
                                                                                                                              tttttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAg                     <  1:252716/150‑1 (MQ=255)
                                                                                                                                tttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGAt                   >  2:1337607/1‑150 (MQ=255)
                                                                                                                                 ttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATc                  <  1:825564/150‑1 (MQ=255)
                                                                                                                                 ttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGGCATCGCTTTCGGAAAGATc                  >  1:1261463/1‑150 (MQ=255)
                                                                                                                                 ttGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGGCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATc                  <  1:212037/150‑1 (MQ=255)
                                                                                                                                   gTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGc                >  2:1547974/1‑150 (MQ=255)
                                                                                                                                    tAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGcc               <  1:541941/150‑1 (MQ=255)
                                                                                                                                    aatgTCGGATGCGGCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGcc               >  2:222008/4‑150 (MQ=39)
                                                                                                                                     atgTCGGATGCGGCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCa              >  2:1024299/3‑150 (MQ=39)
                                                                                                                                      ggTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGcc               >  2:447113/1‑148 (MQ=255)
                                                                                                                                      ggTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAc             >  2:482740/1‑150 (MQ=255)
                                                                                                                                          ggatgcggCGTGAACGCCTTATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCActgct         >  2:143818/7‑150 (MQ=39)
                                                                                                                                            aTAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTgg       <  1:137341/150‑1 (MQ=255)
                                                                                                                                                 gCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGAtt  >  1:1567915/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCACTGCTGGGGATT  >  NC_000913/111415‑111709

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