Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,458:1 +T 5.6% intergenic (+90/+23) yncE → / ← ansP ATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,524,4581.T 5.6% 671.9 / 35.3 161intergenic (+90/+23)yncE/ansPATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base . (84/68);  new base T (3/6);  total (87/74)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 165 total aligned reads displayed.

AATCCACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAA  >  NC_000913/1524310‑1524604
                                                                                                                                                     |                                                                                                                                                  
aaTCCACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGc                                                                                                                                                     <  1:1443034/149‑1 (MQ=255)
  tCCACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCAc                                                                                                                                                  >  2:85130/1‑149 (MQ=255)
   ccACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACg                                                                                                                                                 <  3:533448/149‑1 (MQ=255)
    cACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGt                                                                                                                                                >  4:729903/1‑149 (MQ=255)
       tAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGc                                                                                                                                             >  1:826746/1‑149 (MQ=255)
         aaCAGCAGTAAGCTATCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCg                                                                                                                                           <  2:390412/149‑1 (MQ=255)
         aaCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCg                                                                                                                                           <  4:1000462/149‑1 (MQ=255)
               agGAAGCTACTCAGCCAGACGATGTGATTAGTATTGCGCTGTAATAAATATGGCGATGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTAGGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAAt                                                                                                                                     <  1:1127763/149‑1 (MQ=255)
                  aaGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCAAattc                                                                                                                                  <  4:889216/149‑2 (MQ=255)
                  aaGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAAtttt                                                                                                                                  <  3:1403421/149‑1 (MQ=255)
                   aGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCAAattca                                                                                                                                 <  2:384378/149‑3 (MQ=255)
                   aGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAAtttt                                                                                                                                  >  3:469465/1‑148 (MQ=255)
                    gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCaaattcag                                                                                                                                <  1:1341974/149‑7 (MQ=255)
                    gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGc                                                                                                                                <  4:832084/149‑1 (MQ=255)
                      tACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCaa                                                                                                                              <  1:85130/149‑1 (MQ=255)
                       aCCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAAt                                                                                                                             >  2:64027/1‑149 (MQ=255)
                        cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAAtt                                                                                                                            <  3:1066994/149‑1 (MQ=255)
                           agccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAAc                                                                                                                         >  1:728592/1‑149 (MQ=255)
                           agccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAAc                                                                                                                         >  2:1161804/1‑149 (MQ=255)
                            gccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTaa                                                                                                                          >  3:297988/1‑147 (MQ=255)
                            gccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACt                                                                                                                        >  2:1374811/1‑149 (MQ=255)
                             ccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTc                                                                                                                       <  2:150705/149‑1 (MQ=255)
                              cagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCg                                                                                                                      >  2:1048053/1‑149 (MQ=255)
                                gACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGcc                                                                                                                    <  4:16516/149‑1 (MQ=255)
                                   gATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGtt                                                                                                                 >  3:8352/1‑149 (MQ=255)
                                   gATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGtt                                                                                                                 >  2:1065166/1‑149 (MQ=255)
                                    aTGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGttt                                                                                                                >  1:1064706/1‑149 (MQ=255)
                                     tgtgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTc                                                                                                               >  2:464120/1‑149 (MQ=255)
                                       tgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAg                                                                                                             <  4:1089927/149‑1 (MQ=255)
                                       tgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAg                                                                                                             <  2:1016344/149‑1 (MQ=255)
                                        gATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                            <  2:561766/149‑1 (MQ=255)
                                        gATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                            <  1:640124/149‑1 (MQ=255)
                                              tATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaaa                                                                                                      <  1:688276/149‑1 (MQ=255)
                                              tATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaaa                                                                                                      <  1:1105698/149‑1 (MQ=255)
                                                   cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                 <  3:848942/149‑1 (MQ=255)
                                                   cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                 >  1:1315755/1‑149 (MQ=255)
                                                      tGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAttt                                                                                              >  4:769129/1‑149 (MQ=255)
                                                        tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                            >  3:472327/1‑149 (MQ=255)
                                                          aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                          >  1:682570/1‑149 (MQ=255)
                                                          aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGCAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                          >  4:1212472/1‑149 (MQ=255)
                                                             gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGttt                                                                                       >  3:64543/1‑149 (MQ=255)
                                                             gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGttt                                                                                       <  1:1255421/149‑1 (MQ=255)
                                                              atatGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtttt                                                                                      <  4:52713/149‑1 (MQ=255)
                                                                atGGCGAGGGTGAAAACGGTGAGACTTGCCTGCTTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttc                                                                                    <  3:981584/149‑1 (MQ=255)
                                                                atGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttc                                                                                    <  2:1308152/149‑1 (MQ=255)
                                                                 tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                   <  4:163583/149‑1 (MQ=255)
                                                                       gggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCt                                                                             >  1:583765/1‑149 (MQ=255)
                                                                         gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTc                                                                           >  4:749122/1‑149 (MQ=255)
                                                                             aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcg                                                                       >  4:1130217/1‑149 (MQ=255)
                                                                              aaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcga                                                                      >  3:328629/1‑149 (MQ=255)
                                                                                     agaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTgg                                                               <  4:1426783/149‑1 (MQ=255)
                                                                                           gcctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGt                                                         <  4:299448/149‑1 (MQ=255)
                                                                                             ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                                                                                >  2:1217009/1‑108 (MQ=255)
                                                                                             ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                                                                                <  1:1217009/108‑1 (MQ=255)
                                                                                              tgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCt                                                      >  2:243565/1‑149 (MQ=255)
                                                                                               gcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttctt                                                                                  <  2:684023/120‑1 (MQ=255)
                                                                                                cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCAt                                                                               >  3:66957/1‑122 (MQ=255)
                                                                                                cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgt                                                    >  1:680430/1‑149 (MQ=255)
                                                                                                 ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                          >  3:764011/1‑110 (MQ=255)
                                                                                                 ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                   <  3:948213/149‑1 (MQ=255)
                                                                                                  ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                          >  1:157657/1‑125 (MQ=255)
                                                                                                  ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                          <  2:157657/125‑1 (MQ=255)
                                                                                                  ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                  >  3:1215774/1‑149 (MQ=255)
                                                                                                  ttttATCAAGGCGTATCCGGCACTTTCGTCACTAATCTGAATTTGCCCAGCACGTAGCCGGGCAATTTTGAAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtta                                                  >  2:1329637/1‑147 (MQ=255)
                                                                                                       tCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTc                                             >  3:86209/1‑149 (MQ=255)
                                                                                                        cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                            <  3:1212472/149‑1 (MQ=255)
                                                                                                        cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                            <  2:962375/149‑1 (MQ=255)
                                                                                                         aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                           >  4:1318248/1‑149 (MQ=255)
                                                                                                            gCGTATCCGGCACTTTCGTCATTAATCTGAATTTGCCCGGCATGTAGCCGGGCAATTTTGCAATTTACTCGACGTTTCAGGCTTAAACACATTTTGCCGTCTTTGTTCATCTTCCTCCACGACTGGCGCAGTGCTGTGAAGTTCAGCag                                        <  4:250273/149‑2 (MQ=255)
                                                                                                               tATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAcgc                                     >  1:698896/1‑149 (MQ=255)
                                                                                                                   cGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTGCATCTGCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGttt                                 <  4:939247/149‑1 (MQ=255)
                                                                                                                    ggCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTg                                >  4:575021/1‑149 (MQ=255)
                                                                                                                     gCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTATTCATCTTCCTCGACGACCGGCGCAGTGCTGTGAATTTCAGCAATGCGTTTgc                               <  2:591578/149‑1 (MQ=255)
                                                                                                                        cTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCACTTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTAGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGca                            <  4:911691/149‑1 (MQ=255)
                                                                                                                         tttCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                           >  4:350773/1‑149 (MQ=255)
                                                                                                                               cAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCaaa                     <  4:817326/149‑1 (MQ=255)
                                                                                                                                 gTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGttt                                 >  1:1292852/1‑135 (MQ=255)
                                                                                                                                 gTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAcc                   >  1:1141375/1‑149 (MQ=255)
                                                                                                                                    aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                <  1:371887/149‑1 (MQ=255)
                                                                                                                                    aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                <  2:725210/149‑1 (MQ=255)
                                                                                                                                    aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                >  4:578416/1‑149 (MQ=255)
                                                                                                                                    aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                <  4:64543/149‑1 (MQ=255)
                                                                                                                                     tCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTGCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc               >  2:672948/1‑149 (MQ=255)
                                                                                                                                     tCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc               >  4:1367254/1‑149 (MQ=255)
                                                                                                                                     tCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc               >  2:468735/1‑149 (MQ=255)
                                                                                                                                     tCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc               >  1:498565/1‑149 (MQ=255)
                                                                                                                                       tGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCTGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta             >  2:1291490/1‑149 (MQ=255)
                                                                                                                                       tGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta             >  1:1406240/1‑149 (MQ=255)
                                                                                                                                        gAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCtat            <  1:464120/149‑1 (MQ=255)
                                                                                                                                          aTTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATaa          <  4:86209/149‑1 (MQ=255)
                                                                                                                                           tttGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                >  1:394339/1‑142 (MQ=255)
                                                                                                                                            ttGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                           >  3:185582/1‑130 (MQ=255)
                                                                                                                                             tGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                   <  2:683536/73‑1 (MQ=38)
                                                                                                                                             tGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                   >  1:683536/1‑73 (MQ=38)
                                                                                                                                               cccGGCACGTAGCCGTGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCGTCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc     <  1:1117022/149‑1 (MQ=255)
                                                                                                                                                ccGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                                                                    <  4:1331537/53‑1 (MQ=255)
                                                                                                                                                ccGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                                                                    >  3:1331537/1‑53 (MQ=255)
                                                                                                                                                ccGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa    >  4:1446260/1‑149 (MQ=255)
                                                                                                                                                ccGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCa    >  3:1406100/1‑149 (MQ=255)
                                                                                                                                                  ggCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaaa  >  2:1093976/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AATCCACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAA  >  NC_000913/1524310‑1524604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: