Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 25,774 A→C 5.9% intergenic (+73/‑52) lspA → / → fkpB prolipoprotein signal peptidase (signal peptidase II)/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,7740AC 5.9% 384.4 / 12.3 136intergenic (+73/‑52)lspA/fkpBprolipoprotein signal peptidase (signal peptidase II)/FKBP‑type peptidyl‑prolyl cis‑trans isomerase (rotamase)
Reads supporting (aligned to +/- strand):  ref base A (69/59);  new base C (2/6);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 1.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 137 total aligned reads displayed.

GATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTA  >  NC_000913/25627‑25922
                                                                                                                                                   |                                                                                                                                                      
gATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAg‑‑g                                                                                                                                                    <  1:991985/150‑1 (MQ=255)
  tACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCg                                                                                                                                                  >  1:1371851/1‑150 (MQ=255)
     tGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGttt                                                                                                                                               >  1:478480/1‑150 (MQ=255)
     tGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGtta                                                                                                                                               >  1:28317/1‑149 (MQ=255)
       ccATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTAc                                                                                                                                             <  1:1340986/150‑1 (MQ=255)
            tgtgTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTAcgccgc                                                                                                                                        <  1:1699978/150‑1 (MQ=255)
                tctgTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCG‑‑GCGTAAACGCCTTATcc                                                                                                                                    <  1:278219/147‑1 (MQ=38)
                tCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCG‑‑GCGTAAACGCCTTATcc                                                                                                                                    <  1:1602984/150‑1 (MQ=39)
                tCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATcc                                                                                                                                    <  2:1498243/150‑1 (MQ=255)
                   gTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGc                                                                                                                                 >  2:265649/1‑150 (MQ=255)
                    tGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCG‑‑GCGTAAACGCCTTATCCGGcc                                                                                                                                <  1:523646/150‑2 (MQ=39)
                      cGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCaaa                                                                                                                              >  2:41431/1‑150 (MQ=255)
                       ggCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGTGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCG‑‑GCGTAAACGCCTTATCCGGcctac                                                                                                                             <  2:961893/150‑5 (MQ=255)
                       ggCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCTAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCaaaa                                                                                                                             >  1:1183555/1‑150 (MQ=255)
                       ggCACTGATTGTGCTGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCaaaaa                                                                                                                            >  2:1341638/1‑150 (MQ=255)
                          aCTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGTATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATGCG‑‑GCGTAAACGCCTTATCCGGcctacgat                                                                                                                          <  2:1487386/150‑8 (MQ=38)
                          aCTGAATGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGACTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATc                                                                                                                          <  2:1117144/150‑1 (MQ=255)
                            tGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCt                                                                                                                        >  1:548711/1‑150 (MQ=255)
                              aTTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTa                                                                                                                      <  1:1554914/150‑1 (MQ=255)
                                   gCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAatat                                                                                                                 >  1:1011245/1‑150 (MQ=255)
                                   gCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAatat                                                                                                                 >  1:661948/1‑150 (MQ=255)
                                   gCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAatat                                                                                                                 >  2:1694093/1‑150 (MQ=255)
                                    cTGGAAGGTTTTTTGCCTTCTACACCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACCCATCTTATCCGGCCTACAGATTGCTTCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCATAAAtata                                                                                                                >  2:1256737/1‑150 (MQ=37)
                                      ggAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                                                                                                              >  1:1350155/1‑150 (MQ=255)
                                      ggAAGGTTTTTTGACTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAag                                                                                                              <  1:742794/150‑1 (MQ=255)
                                       gAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGTCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                    >  2:1186562/1‑143 (MQ=255)
                                       gAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGTCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                    <  1:1152459/143‑1 (MQ=255)
                                       gAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGTCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                    <  1:1186562/143‑1 (MQ=255)
                                       gAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGTCGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTaaa                                                                                                                    >  2:1152459/1‑143 (MQ=255)
                                          ggTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCa                                                                                                          <  1:1484450/150‑1 (MQ=255)
                                           gTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCaa                                                                                                         >  1:1429347/1‑150 (MQ=255)
                                            ttttttGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCaaa                                                                                                        >  2:1455582/1‑150 (MQ=255)
                                             tttttGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAAc                                                                                                       <  1:447299/150‑1 (MQ=255)
                                                  gCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCa                                                                                                  >  1:967578/1‑150 (MQ=255)
                                                    cTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATg                                                                                                <  2:463114/150‑1 (MQ=255)
                                                        tAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTg                                                                                            <  2:538230/150‑1 (MQ=255)
                                                           agCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAAt                                                                                         >  2:20467/1‑150 (MQ=255)
                                                             cGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCt                                                                                       >  2:1067159/1‑150 (MQ=255)
                                                             cGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCt                                                                                       >  1:582288/1‑150 (MQ=255)
                                                                aaaaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTa                                                                                    <  1:390743/150‑1 (MQ=255)
                                                                aaaaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTa                                                                                    <  1:384404/150‑1 (MQ=255)
                                                                aaaaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTa                                                                                    <  2:1709610/150‑1 (MQ=255)
                                                                 aaaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAc                                                                                   >  2:539084/1‑150 (MQ=255)
                                                                 aaaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTAc                                                                                   <  2:1390974/150‑1 (MQ=255)
                                                                  aaaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAATCCTGCATGTCTGAATCTGTACa                                                                                  <  2:1456178/150‑1 (MQ=255)
                                                                   aaCAATAATAAACCCTGCCGTATGCGATGCTGACGCATCTTATACGTCCTACAGATTTCTGAGAAATCGTAGTCCGGATAAG‑‑TCGTTTTCGCCGCATACGGCAAAAATCCTTAAATATATGAGCAAACCTGCATGTCTGAATCTGTACag                                                                                 <  1:424327/150‑1 (MQ=255)
                                                                   aaCAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCGGGCCGACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAACTCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACag                                                                                 <  2:865071/150‑1 (MQ=255)
                                                                     cAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAAGCTGCATGTCTGAATCTCTACagag                                                                               >  1:265793/1‑150 (MQ=255)
                                                                     cAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACagag                                                                               >  2:1621316/1‑150 (MQ=255)
                                                                     cAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATGGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAACGCCTTAAATATAAGAGCACACCTGCATGTCTGAATCTGTACagcg                                                                               >  1:1662312/1‑148 (MQ=255)
                                                                     aaataatAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACagag                                                                               <  2:1415606/149‑1 (MQ=255)
                                                                           taaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAg                                                                         <  1:174930/150‑1 (MQ=255)
                                                                           taaACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GAGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAg                                                                         >  2:779049/1‑150 (MQ=255)
                                                                                     cGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTg                                                               >  2:912787/1‑150 (MQ=255)
                                                                                      ggATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTgg                                                              >  2:1180678/1‑150 (MQ=255)
                                                                                        atgcgatgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTg                                                            <  2:310211/150‑1 (MQ=255)
                                                                                         tgcgatgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACGGAGCAATAGCGCCGTCCTGGTGc                                                           <  1:1341638/150‑1 (MQ=255)
                                                                                           cgatgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAGACCTGCATGTCTGAATCAGTACAGAGCAATAGCGCCGTCCTGGTGCag                                                         >  1:1072115/1‑149 (MQ=255)
                                                                                           cgatgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCAc                                                         >  2:1393170/1‑150 (MQ=255)
                                                                                            gatgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACt                                                        >  1:394006/1‑150 (MQ=255)
                                                                                              tgcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGGTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTc                                                      <  2:117298/150‑1 (MQ=255)
                                                                                               gcTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTcc                                                     >  1:711120/1‑149 (MQ=255)
                                                                                                cTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCAc                                                    <  1:1465071/150‑1 (MQ=255)
                                                                                                cTGACGCAGCTTATCCGGCCTACAGATTGCTGCGATATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCACACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCAc                                                    <  1:1297020/150‑1 (MQ=255)
                                                                                                 tGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACg                                                   <  2:1394189/150‑1 (MQ=255)
                                                                                                 tGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACg                                                   >  1:512388/1‑150 (MQ=255)
                                                                                                  gACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGc                                                  >  2:1400168/1‑150 (MQ=255)
                                                                                                       aTCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTaaaa                                             >  1:764315/1‑150 (MQ=255)
                                                                                                             tCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcgac                                       <  1:868930/150‑1 (MQ=255)
                                                                                                             tCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcgac                                       <  2:1587431/150‑1 (MQ=255)
                                                                                                              ccGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTcgacg                                      <  1:1021485/150‑1 (MQ=255)
                                                                                                                ggCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGAt                                    >  2:1008820/1‑150 (MQ=255)
                                                                                                                 gCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCCCTTCACGCTAAAACTCGACGATg                                   >  2:319688/1‑150 (MQ=255)
                                                                                                                 gCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATg                                   >  1:1157569/1‑150 (MQ=255)
                                                                                                                  ccTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATgg                                  <  1:1694093/150‑1 (MQ=255)
                                                                                                                    tACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCCCGCTAAAACTCGACGATGGca                                >  1:931245/1‑150 (MQ=255)
                                                                                                                       aGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGcacca                             <  2:47412/150‑1 (MQ=255)
                                                                                                                       aGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGcacca                             >  2:189221/1‑150 (MQ=255)
                                                                                                                            gctgcGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCAccgcc                        >  1:501339/1‑150 (MQ=255)
                                                                                                                             ctgcGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCTTGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCCCGCTAAAACTCGACGATGGCACCAccgccg                       >  1:1511882/1‑150 (MQ=255)
                                                                                                                               gcGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATgg                                  <  1:1658899/137‑1 (MQ=255)
                                                                                                                               gcGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATgg                                  >  2:1658899/1‑137 (MQ=255)
                                                                                                                               gcGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAg                     <  1:1173365/150‑1 (MQ=255)
                                                                                                                                 gAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGGCGGTGGTGCCATCGTCGAGtt                   >  2:1496074/1‑149 (MQ=38)
                                                                                                                                  aatctgTAGGCCGGATAAGATGCGTCAGCATCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGt                    <  1:713149/145‑1 (MQ=37)
                                                                                                                                   aaTCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTa                 <  2:1350155/150‑1 (MQ=255)
                                                                                                                                     tCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTAcc               <  2:1429347/150‑1 (MQ=255)
                                                                                                                                       gTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCg             <  1:1562467/150‑1 (MQ=255)
                                                                                                                                          tgCCGGATGCG‑‑GCGTAAACGCCTTATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTAGACGATGGCAACACCGCCGAGTCTACCCGcaa          >  2:1629373/2‑150 (MQ=255)
                                                                                                                                           gCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGcaac         <  1:176852/150‑1 (MQ=255)
                                                                                                                                            ccGGATGCG‑‑GCGTAAACGCCTTATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGcaac         >  2:669712/1‑149 (MQ=255)
                                                                                                                                             cGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGcaacaa       >  2:636707/1‑150 (MQ=255)
                                                                                                                                             cGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTATGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGcaacaa       >  2:630764/1‑150 (MQ=255)
                                                                                                                                              ggATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCaacaac      <  1:1393170/150‑1 (MQ=255)
                                                                                                                                              ggATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCaacaac      >  2:640696/1‑150 (MQ=255)
                                                                                                                                              ggATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCaacaac      >  1:129463/1‑150 (MQ=255)
                                                                                                                                                 tAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCt                  >  2:1148084/1‑135 (MQ=255)
                                                                                                                                                 tAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCt                  <  1:1148084/135‑1 (MQ=255)
                                                                                                                                                 tAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGt   <  1:1709947/150‑1 (MQ=255)
                                                                                                                                                  aaG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTa  <  1:481763/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
GATACTGCCATCTGTGTCGGTGCGGCACTGATTGTGCTGGAAGGTTTTTTGCCTTCTAGAGCGAAAAAACAATAATAAACCCTGCCGGATGCGATGCTGACGCATCTTATCCGGCCTACAGATTGCTGCGAAATCGTAGGCCGGATAAG‑‑GCGTTTACGCCGCATCCGGCAAAAATCCTTAAATATAAGAGCAAACCTGCATGTCTGAATCTGTACAGAGCAATAGCGCCGTCCTGGTGCACTTCACGCTAAAACTCGACGATGGCACCACCGCCGAGTCTACCCGCAACAACGGTA  >  NC_000913/25627‑25922

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