Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,102,951 G→A 5.5% I82I (ATC→ATT csgD ← csgBAC operon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,9510GA 5.5% 220.8 / 2.0 73I82I (ATC→ATTcsgDcsgBAC operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/32);  new base A (2/2);  total (39/34)
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

GTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCA  >  NC_000913/1102802‑1103096
                                                                                                                                                     |                                                                                                                                                 
gTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTg                                                                                                                                                   >  1:425260/1‑150 (MQ=255)
gTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTg                                                                                                                                                   <  2:23151/150‑1 (MQ=255)
gTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTg                                                                                                                                                   <  2:479712/150‑1 (MQ=255)
    aGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACTCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACAGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGatat                                                                                                                                               >  2:174270/1‑150 (MQ=255)
       tgCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGTAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgt                                                                                                                                            <  1:355326/149‑1 (MQ=255)
         cATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATTGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATAttgttg                                                                                                                                          <  2:485560/150‑1 (MQ=255)
          aTTCGCCGCGCAGGACGCCTTTCAACCCATTGACAACACGTTCTTTATCCTCCAGGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgc                                                                                                                                         <  1:956720/150‑2 (MQ=255)
            tCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATCAAAATAttgttgcta                                                                                                                                       <  2:921453/150‑4 (MQ=255)
              gccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAACACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGttttt                                                                                                                                     <  1:1074450/150‑1 (MQ=255)
               ccgcGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTGTTCAATGTCGCGGTACGGGTAATCTTCCGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTc                                                                                                                                    >  1:979027/1‑150 (MQ=255)
                cgcgCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTcc                                                                                                                                   <  2:437667/150‑1 (MQ=255)
                   gcAGGACGACTTGCAACCCATTGACAACACGTTCTTGATGCTCCATGGAATCAAAAACGCCGTTGATATGAGTCCAGTTTTCAATGTCGCGGTACGGTTAAGCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGc                                                                                                                                <  1:62774/150‑1 (MQ=255)
                     aGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTc                                                                                                                              <  2:213510/150‑1 (MQ=255)
                      ggACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCa                                                                                                                             >  2:318412/1‑150 (MQ=255)
                        aCGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCATGCTCaaa                                                                                                                           >  2:849623/1‑150 (MQ=255)
                           ccTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTa                                                                                                                        >  1:551446/1‑150 (MQ=255)
                           ccTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTa                                                                                                                        >  1:601649/1‑150 (MQ=255)
                               gCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCt                                                                                                                    <  2:598748/150‑1 (MQ=255)
                                  aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGcc                                                                                                                 <  2:1252989/150‑1 (MQ=255)
                                   cccATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCa                                                                                                                <  2:403975/150‑1 (MQ=255)
                                   cccATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCa                                                                                                                >  1:1030330/1‑150 (MQ=255)
                                    ccATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGGCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCaa                                                                                                               >  2:466956/1‑150 (MQ=255)
                                            aacaCGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGAt                                                                                                       <  1:1034927/150‑1 (MQ=255)
                                             acacGTTCTTGATCCTCCATGGATTAAAACACGCCGTTGATATGCGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATa                                                                                                      >  2:244080/1‑150 (MQ=255)
                                               acGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAg                                                                                                    >  1:797175/1‑150 (MQ=255)
                                                cGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGc                                                                                                   >  1:533264/1‑150 (MQ=255)
                                                cGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGCTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTAGCCTGCCAATAATGGATATGc                                                                                                   >  2:841816/1‑150 (MQ=255)
                                                 gTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGTATAAGCt                                                                                                  >  2:196039/1‑150 (MQ=255)
                                                 gTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGAACGGGTAATCTTCAGGGGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCt                                                                                                  >  2:1214443/1‑150 (MQ=255)
                                                   tCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCttt                                                                                                >  2:1208748/1‑150 (MQ=255)
                                                      tGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTa                                                                                             <  2:516449/150‑1 (MQ=255)
                                                               aTGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTcc                                                                                    <  1:1084723/150‑1 (MQ=255)
                                                                 ggAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcat                                                                                  >  2:259607/1‑150 (MQ=255)
                                                                  gaagaaAAAACGCCGTTGATAAGAGGCGAGTTTTCAATGTCGCGGAACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatc                                                                                 >  2:781692/1‑150 (MQ=255)
                                                                  gAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatc                                                                                 >  1:90118/1‑150 (MQ=255)
                                                                   aaTAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcataa                                                                                <  1:753059/150‑3 (MQ=255)
                                                                      aaaaaaCGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCatat                                                                             >  2:851985/1‑150 (MQ=255)
                                                                           aCGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCaga                                                                        <  2:602310/150‑1 (MQ=255)
                                                                             gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  1:1234341/150‑1 (MQ=255)
                                                                             gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  2:979607/150‑1 (MQ=255)
                                                                             gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  2:1025017/150‑1 (MQ=255)
                                                                             gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  2:547563/150‑1 (MQ=255)
                                                                                gTTGATATGAGGCCAGTTTTCAATGTCGCGGTACTTGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaat                                                                   <  1:1063347/150‑1 (MQ=255)
                                                                                  tGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaataa                                                                 <  1:976375/150‑1 (MQ=255)
                                                                                   gATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatca                                                                                <  1:1302369/134‑1 (MQ=255)
                                                                                   gATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatca                                                                                >  2:1302369/1‑134 (MQ=255)
                                                                                        gAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGc                                                           >  1:822643/1‑150 (MQ=255)
                                                                                         aGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGAAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGcc                                                          >  2:931334/1‑150 (MQ=255)
                                                                                         aGGCCAGTTTTCAATGTCGCGCTACGGGAAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAAAATTGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGcc                                                          >  1:276599/1‑150 (MQ=255)
                                                                                          ggCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCt                                                         >  1:612608/1‑150 (MQ=255)
                                                                                             cAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGaa                                                      >  2:1235979/1‑150 (MQ=255)
                                                                                                  ttCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTCTCCTGCCAATAATGGATAAGCTTTTTATCAGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGat                                                 >  2:930668/1‑150 (MQ=255)
                                                                                                  ttCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGat                                                 >  1:352267/1‑150 (MQ=255)
                                                                                                  ttCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGat                                                 >  2:557914/1‑150 (MQ=255)
                                                                                                       tGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGt                                            <  1:916155/150‑1 (MQ=255)
                                                                                                        gTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTc                                           <  1:449226/150‑1 (MQ=255)
                                                                                                           gcgGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCag                                        <  1:930044/150‑1 (MQ=255)
                                                                                                           gcgGTACGGGTAATCTTCAGGCGTATTTAGCAAAAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCag                                        <  2:1146638/150‑1 (MQ=255)
                                                                                                            cgGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATTTCCTGAAGATTTATAGAGCCTGAAGAGATATCGTCCaga                                       >  1:131681/1‑150 (MQ=255)
                                                                                                             ggTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCagag                                      >  2:378854/1‑150 (MQ=255)
                                                                                                                   ggTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGtt                                <  1:196039/150‑1 (MQ=255)
                                                                                                                   ggTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGtt                                >  2:978396/1‑150 (MQ=255)
                                                                                                                      aaTCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGaa                             <  1:1061666/150‑1 (MQ=255)
                                                                                                                        tCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAata                           <  2:577180/150‑1 (MQ=255)
                                                                                                                        tCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAata                           >  2:217310/1‑150 (MQ=255)
                                                                                                                              tgCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGt                     <  1:144299/149‑1 (MQ=255)
                                                                                                                              ggCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCTCTTCCATCATATCCAGAAGAATAAAAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGt                     >  1:493176/1‑150 (MQ=255)
                                                                                                                                 gTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAAt                  >  1:1115934/1‑150 (MQ=255)
                                                                                                                                      tAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTTTCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTcc             <  2:931715/150‑1 (MQ=255)
                                                                                                                                       aacaacaaTATCAAAATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCt            >  1:517483/3‑150 (MQ=255)
                                                                                                                                        gCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTg           <  2:482172/150‑1 (MQ=255)
                                                                                                                                        acaacaaTATCAAAATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATTAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTg           >  2:256590/2‑150 (MQ=255)
                                                                                                                                         caacaaTATTTTGATTTTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGt          >  1:670511/1‑150 (MQ=255)
                                                                                                                                         caacaaTATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGt          >  1:985601/1‑150 (MQ=255)
                                                                                                                                             aaTATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTg     <  2:310515/150‑1 (MQ=255)
                                                                                                                                                 ttttGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCa  >  1:317112/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCA  >  NC_000913/1102802‑1103096

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