Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,828 G→T 5.1% intergenic (+157/+49) loiP → / ← speB metalloprotease LoiP/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,8280GT 5.1% 265.6 / 3.1 79intergenic (+157/+49)loiP/speBmetalloprotease LoiP/agmatinase
Reads supporting (aligned to +/- strand):  ref base G (42/33);  new base T (2/2);  total (44/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

TGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTT  >  NC_000913/3082679‑3082977
                                                                                                                                                     |                                                                                                                                                     
tggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccg                                                                                                                                                       <  2:551302/150‑1 (MQ=255)
   tgtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCAt                                                                                                                                                    <  1:92736/150‑1 (MQ=255)
    gtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                                                                                                                                                   <  1:999781/150‑1 (MQ=255)
      cGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                 >  2:471767/1‑150 (MQ=255)
        aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                                                                                                                                               <  1:1387730/150‑1 (MQ=255)
          cGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAACGCCTTATCCGAcct                                                                                                                                             <  2:913417/150‑3 (MQ=38)
              ggTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                                                                                                                                         >  2:442113/1‑150 (MQ=255)
                   gATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCAc                                                                                                                                    >  1:1697524/1‑150 (MQ=255)
                     tAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGATGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGt                                                                                                                                  >  2:954575/1‑150 (MQ=255)
                     tAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGATGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGt                                                                                                                                  >  2:961814/1‑150 (MQ=255)
                        ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               >  1:339161/1‑150 (MQ=255)
                        ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               <  1:955404/150‑1 (MQ=255)
                        ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               >  1:1080464/1‑150 (MQ=255)
                            tgaaCGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTAccc                                                                                                                           >  2:963308/4‑150 (MQ=255)
                             tCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCg                                                                                                                          >  1:1406980/1‑150 (MQ=255)
                              cACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGt                                                                                                                         >  2:357454/1‑150 (MQ=255)
                                 gccgcATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTgcg                                                                                                                      >  2:690032/1‑150 (MQ=255)
                                      aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 >  2:1853609/1‑150 (MQ=255)
                                      aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  1:357454/150‑1 (MQ=255)
                                      aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  2:177232/150‑1 (MQ=255)
                                      aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  1:431283/150‑1 (MQ=255)
                                      aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGAGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  2:1559327/150‑1 (MQ=255)
                                       tCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgc                                                                                                                <  1:1844002/150‑1 (MQ=255)
                                        cTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgca                                                                                                               <  1:471767/150‑1 (MQ=255)
                                        cTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgca                                                                                                               <  1:581234/150‑1 (MQ=255)
                                         tGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgcat                                                                                                              >  2:983624/1‑150 (MQ=255)
                                             aTGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTgg                                                                                                          >  2:1611551/1‑150 (MQ=255)
                                              tGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGt                                                                                                         >  1:1732238/1‑150 (MQ=255)
                                               gAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCAGTGCGCATCGCATCTGGTg                                                                                                        >  1:1808609/1‑150 (MQ=255)
                                                 aTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCt                                                                                                      <  2:1643599/150‑1 (MQ=255)
                                                 aTGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCt                                                                                                      >  2:77357/1‑150 (MQ=255)
                                                    gCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTAc                                                                                                   >  1:486128/1‑150 (MQ=255)
                                                           tGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCt                                                                                            <  1:688498/150‑1 (MQ=255)
                                                              tACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCtttt                                                                                         <  2:200532/150‑1 (MQ=255)
                                                              acaaatgtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCtttt                                                                                         <  1:1950141/145‑1 (MQ=255)
                                                                 cTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcg                                                                                      <  1:1578990/150‑1 (MQ=255)
                                                                 cTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcg                                                                                      <  2:1237046/150‑1 (MQ=255)
                                                                   tgtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgcc                                                                                    <  2:339161/150‑1 (MQ=255)
                                                                    gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgccg                                                                                   <  1:442113/150‑1 (MQ=255)
                                                                    gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgccg                                                                                   <  1:21394/150‑1 (MQ=255)
                                                                      gcgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                                 <  1:690032/150‑1 (MQ=255)
                                                                      gcgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                                 >  2:1472112/1‑150 (MQ=255)
                                                                        gcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTg                                                                               >  1:738668/1‑150 (MQ=255)
                                                                        gcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTg                                                                               >  1:745376/1‑150 (MQ=255)
                                                                        gcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTg                                                                               >  1:1608417/1‑150 (MQ=255)
                                                                         cATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGa                                                                              >  2:1086506/1‑150 (MQ=255)
                                                                             gCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAatat                                                                          >  1:1814588/1‑150 (MQ=255)
                                                                              ccGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAAtata                                                                         >  2:877952/1‑150 (MQ=255)
                                                                                ggATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACa                                                                       <  2:1697524/150‑1 (MQ=255)
                                                                                ggATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACa                                                                       >  2:817494/1‑150 (MQ=255)
                                                                                ggATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACa                                                                       >  1:1710835/1‑150 (MQ=255)
                                                                                  atgcgaggcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGc                                                                     <  1:270866/150‑1 (MQ=255)
                                                                                     cgatgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCAtt                                                                  >  1:1745641/1‑150 (MQ=255)
                                                                                       atgcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTc                                                                >  2:571472/1‑150 (MQ=255)
                                                                                         gcTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTcca                                                              >  1:1009106/1‑150 (MQ=255)
                                                                                           tGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagc                                                            <  2:1080464/150‑1 (MQ=255)
                                                                                           tGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagc                                                            >  2:1864403/1‑150 (MQ=255)
                                                                                            ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagcg                                                           >  1:1785997/1‑150 (MQ=255)
                                                                                            ggCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTAGCCCTTTTTCGCCGCCTGAATATACAGCATTTccagcg                                                           <  2:1406980/150‑1 (MQ=255)
                                                                                                     tATCCGGCCTACTGGTGGCGAACGTAGGTCGGATAAGGCGTTCAAGCCGCATACGACATTAATGGAACGTTTTACCCGTGCGCATCGAATCTGTTTCTTACTCGCCCTTTTTCGCCGCCTGAATATACATCATTTCCATCGCCAGCGTTg                                                  <  2:738668/150‑1 (MQ=255)
                                                                                                      aTCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGc                                                 >  2:1684734/1‑150 (MQ=255)
                                                                                                           gCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTg                                            <  1:40282/150‑1 (MQ=255)
                                                                                                                  ggTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGc                                     >  1:1010877/1‑150 (MQ=255)
                                                                                                                    tGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAg                                   <  2:1814588/150‑1 (MQ=255)
                                                                                                                       cGAACGTAGGGCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGa                                <  1:1382200/150‑1 (MQ=255)
                                                                                                                            gTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTcc                           <  1:77357/150‑1 (MQ=255)
                                                                                                                             tAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCg                          >  2:428940/1‑150 (MQ=255)
                                                                                                                             aatgTCGGATGCGGCGTGAACGCCTTATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCg                          >  2:1629549/4‑150 (MQ=255)
                                                                                                                                gTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACt                       >  2:495080/1‑150 (MQ=255)
                                                                                                                                   ggATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGAt                    <  2:1920741/150‑1 (MQ=255)
                                                                                                                                    gATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATc                   >  1:26777/1‑150 (MQ=255)
                                                                                                                                    gATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATc                   >  1:574114/1‑150 (MQ=255)
                                                                                                                                     aTAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCg                  <  1:817494/150‑1 (MQ=255)
                                                                                                                                      tAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGt                 >  1:527523/1‑150 (MQ=255)
                                                                                                                                        aGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTAt               >  2:27915/1‑150 (MQ=255)
                                                                                                                                             tgaaCGCCTTATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCgg          >  1:751112/4‑150 (MQ=255)
                                                                                                                                              tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGa         <  1:983624/150‑1 (MQ=255)
                                                                                                                                              tCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACACCATTTCCAGCGCCAGCGTTCCCGCTGCAAGAGCAGTGATTTCCGACTGATCGTATGCCGGa         >  1:1954701/1‑150 (MQ=255)
                                                                                                                                               cACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAg        >  2:1575282/1‑150 (MQ=255)
                                                                                                                                               aaCGCCTTATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAg        <  1:1629549/149‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGcc      <  1:1853609/150‑1 (MQ=255)
                                                                                                                                                 cgccgcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGcc      >  1:1462699/1‑150 (MQ=255)
                                                                                                                                                     gcATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCAGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGCCACtt  >  1:608090/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTT  >  NC_000913/3082679‑3082977

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: