Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,102,948 T→C 5.3% K83K (AAA→AAG) ‡ csgD ← csgBAC operon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,102,9480TC 5.3% 228.2 / 3.2 75K83K (AAA→AAG) ‡csgDcsgBAC operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01

GCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCA  >  NC_000913/1102800‑1103096
                                                                                                                                                    |                                                                                                                                                    
gCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATTGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAttt                                                                                                                                                     <  2:1105545/150‑1 (MQ=255)
gCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAttt                                                                                                                                                     >  2:1215824/1‑150 (MQ=255)
gCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAttt                                                                                                                                                     <  2:85067/150‑1 (MQ=255)
 cGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATAtttt                                                                                                                                                    >  2:449226/1‑150 (MQ=255)
  gTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTg                                                                                                                                                   <  2:23151/150‑1 (MQ=255)
  gTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTg                                                                                                                                                   >  1:425260/1‑150 (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:601649/1‑150 (MQ=255)
                             ccTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTa                                                                                                                        >  1:551446/1‑150 (MQ=255)
                                 gCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCt                                                                                                                    <  2:598748/150‑1 (MQ=255)
                                    aCCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGcc                                                                                                                 <  2:1252989/150‑1 (MQ=255)
                                     cccATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCa                                                                                                                >  1:1030330/1‑150 (MQ=255)
                                     cccATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCa                                                                                                                <  2:403975/150‑1 (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                                                                      <  2:979607/150‑1 (MQ=255)
                                                                               gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  2:1025017/150‑1 (MQ=255)
                                                                               gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  2:547563/150‑1 (MQ=255)
                                                                               gCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCagaag                                                                      <  1:1234341/150‑1 (MQ=255)
                                                                                  gTTGATATGAGGCCAGTTTTCAATGTCGCGGTACTTGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaat                                                                   <  1:1063347/150‑1 (MQ=255)
                                                                                    tGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGaataa                                                                 <  1:976375/150‑1 (MQ=255)
                                                                                     gATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatca                                                                                >  2:1302369/1‑134 (MQ=255)
                                                                                     gATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCcatca                                                                                <  1:1302369/134‑1 (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)
                                                                                                                                                    |                                                                                                                                                    
GCGTAAAGTAGCATTCGCCGCGCAGGACGCCTTGCAACCCATTGACAACACGTTCTTGATCCTCCATGGAATAAAAAACGCCGTTGATATGAGGCCAGTTTTCAATGTCGCGGTACGGGTAATCTTCAGGCGTATTTAGCAACAATATTTTGATATTGTTGTTTTTCCTGCTCAAAGTATCCTGCCAATAATGGATAAGCTTTTTATCCGCTTCCATCATATCCAGAAGAATAATAGAGCCTGAAGAGATATCGTCCAGAGAACGTTGAATATTATGTAATTTTCCTGTAATTGCCA  >  NC_000913/1102800‑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: