Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,829 C→T 5.0% intergenic (+158/+48) loiP → / ← speB metalloprotease LoiP/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,8290CT 5.0% 256.4 / 2.0 79intergenic (+158/+48)loiP/speBmetalloprotease LoiP/agmatinase
Reads supporting (aligned to +/- strand):  ref base C (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 = 8.49e-01

TGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTT  >  NC_000913/3082682‑3082977
                                                                                                                                                   |                                                                                                                                                    
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:961814/1‑150 (MQ=255)
                  tAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGATGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGt                                                                                                                                  >  2:954575/1‑150 (MQ=255)
                     ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               <  1:955404/150‑1 (MQ=255)
                     ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               >  1:1080464/1‑150 (MQ=255)
                     ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt                                                                                                                               >  1:339161/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:431283/150‑1 (MQ=255)
                                   aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  2:177232/150‑1 (MQ=255)
                                   aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  1:357454/150‑1 (MQ=255)
                                   aTCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGAGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcg                                                                                                                 <  2:1559327/150‑1 (MQ=255)
                                    tCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgc                                                                                                                <  1:1844002/150‑1 (MQ=255)
                                     cTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgca                                                                                                               <  1:581234/150‑1 (MQ=255)
                                     cTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcatcgca                                                                                                               <  1:471767/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                                                                                      <  2:1237046/150‑1 (MQ=255)
                                                              cTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcg                                                                                      <  1:1578990/150‑1 (MQ=255)
                                                                tgtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgcc                                                                                    <  2:339161/150‑1 (MQ=255)
                                                                 gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgccg                                                                                   <  1:21394/150‑1 (MQ=255)
                                                                 gtgCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTcgccg                                                                                   <  1:442113/150‑1 (MQ=255)
                                                                   gcgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                                 <  1:690032/150‑1 (MQ=255)
                                                                   gcgcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCgccgcc                                                                                 >  2:1472112/1‑150 (MQ=255)
                                                                     gcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTg                                                                               >  1:745376/1‑150 (MQ=255)
                                                                     gcATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTg                                                                               >  1:738668/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                                                                       >  1:1710835/1‑150 (MQ=255)
                                                                             ggATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACa                                                                       <  2:1697524/150‑1 (MQ=255)
                                                                             ggATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACa                                                                       >  2:817494/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:574114/1‑150 (MQ=255)
                                                                                                                                 gATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATc                   >  1:26777/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)
                                                                                                                                                   |                                                                                                                                                    
TGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCATCGCATCTGGTGCTTACTCGCCCTTTTTCGCCGCCTGAATATACAGCATTTCCAGCGCCAGCGTTGCCGCTGCCAGAGCAGTGATTTCCGACTGATCGTATGCCGGAGCCACTT  >  NC_000913/3082682‑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: