Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,702 A→C 5.1% intergenic (+31/+175) loiP → / ← speB Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,082,7020AC 5.1% 1103.1 / 35.5 303intergenic (+31/+175)loiP/speBPhe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase
Reads supporting (aligned to +/- strand):  ref base A (155/132);  new base C (6/10);  total (161/142)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 316 total aligned reads displayed.

TTGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTT  >  NC_000913/3082553‑3082851
                                                                                                                                                     |                                                                                                                                                     
 tGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGTATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAg                                                                                                                                                      >  2:3497979/1‑150 (MQ=255)
 tGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAg                                                                                                                                                      >  1:96954/1‑150 (MQ=255)
  gAAATACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCGTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCTCATGAGCGCGGATGGGATTAAGTAATGCCTGGTGGTGTCGAACGTAGGTCAGATAAgg                                                                                                                                                     >  2:2961618/1‑150 (MQ=255)
   aaaaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGc                                                                                                                                                    >  1:4090533/1‑150 (MQ=255)
    aaaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCg                                                                                                                                                   <  2:2502202/150‑1 (MQ=255)
    aaaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCg                                                                                                                                                   <  1:1854992/150‑1 (MQ=255)
     aaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATCTTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGt                                                                                                                                                  >  2:518780/1‑150 (MQ=255)
     aaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGt                                                                                                                                                  >  1:2931189/1‑150 (MQ=255)
     aaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGt                                                                                                                                                  >  1:1502896/1‑150 (MQ=255)
     aaaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGt                                                                                                                                                  <  2:1223685/150‑1 (MQ=255)
     aaaCAGGCAAAACTCGAAGAACGACGCCAAAGCACAATGTTTGACGACCATCCAGCATCCGCCGAACGCGCCCAGCATAATCGCGATCGCAAGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGt                                                                                                                                                  <  2:3830126/150‑1 (MQ=255)
      aaCTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGtt                                                                                                                                                 >  1:268770/1‑150 (MQ=255)
          ggCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcg                                                                                                                                             <  2:1753021/150‑1 (MQ=255)
           gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                            >  1:3819668/1‑150 (MQ=255)
           gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                            <  1:2128714/150‑1 (MQ=255)
           gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgc                                                                                                                                            >  2:894106/1‑150 (MQ=255)
             aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                          >  2:3175075/1‑150 (MQ=255)
             aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                          >  2:3622141/1‑150 (MQ=255)
             aaaaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                          >  1:1481520/1‑150 (MQ=255)
               aaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCa                                                                                                                                        >  1:3863064/1‑150 (MQ=255)
               aaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCa                                                                                                                                        >  1:2980146/1‑150 (MQ=255)
               aaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCa                                                                                                                                        >  1:3314638/1‑150 (MQ=255)
                aCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAAcgccttat                                                                                                                                       <  2:1928892/150‑5 (MQ=255)
                 cTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAAcgccttatc                                                                                                                                      <  1:1258039/150‑6 (MQ=255)
                 cTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATc                                                                                                                                      >  2:1154885/1‑150 (MQ=255)
                    gaagaaGGTCGCCATAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGa                                                                                                                                   <  1:347171/150‑1 (MQ=255)
                     aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                  <  1:1636552/150‑1 (MQ=255)
                      agaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                 <  1:291525/150‑1 (MQ=255)
                           gTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGCGCGCGGATGGGATTAAGTAACGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAAt                                                                                                                            <  1:3161302/150‑1 (MQ=255)
                            tCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATg                                                                                                                           <  1:623153/150‑1 (MQ=255)
                             cGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                          <  2:2592976/150‑1 (MQ=255)
                             cGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                          <  2:62030/150‑1 (MQ=255)
                              gCCAAAGCTCAATGTTTGACGACCACCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCTTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGc                                                                                                                         <  2:3433018/150‑1 (MQ=255)
                                cAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcac                                                                                                                       >  2:1117196/1‑150 (MQ=255)
                                 aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                      <  1:3680346/150‑1 (MQ=255)
                                 aaaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                      <  1:2774468/150‑1 (MQ=255)
                                           gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCt                                                                                                            >  2:1573616/1‑150 (MQ=255)
                                           gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCt                                                                                                            >  2:1970751/1‑150 (MQ=255)
                                                    ccATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGCTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCAt                                                                                                   >  2:1394875/1‑150 (MQ=255)
                                                           gCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                            <  1:2577220/150‑1 (MQ=255)
                                                            cATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgat                                                                                           <  2:2782534/150‑1 (MQ=255)
                                                             aTCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatg                                                                                          <  2:1731630/150‑1 (MQ=255)
                                                              tCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgc                                                                                         >  1:3794117/1‑150 (MQ=255)
                                                              tCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgc                                                                                         >  1:1548841/1‑150 (MQ=255)
                                                                 gccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgat                                                                                      >  1:333076/1‑150 (MQ=255)
                                                                  ccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATTTGATTTAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgatg                                                                                     <  2:3346042/150‑1 (MQ=255)
                                                                   cgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgatgc                                                                                    >  1:3754196/1‑150 (MQ=255)
                                                                            ccAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                           <  1:4946/150‑1 (MQ=255)
                                                                            ccAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                           <  1:693683/150‑1 (MQ=255)
                                                                               gCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTa                                                                        >  1:1610105/1‑150 (MQ=255)
                                                                                    ttCGCGATCGCATGAGCGCGGATGGGATGAAGTAAAGCCTGGTGGTGGCGAACGTAGGTCAGATCAGGCGTTCACGCCGCATCTGACAGGAAGGGCACAGTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCgg                                                                   >  2:2284665/1‑150 (MQ=255)
                                                                                         gATCGCATGAGCGCGGATGGTATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                              <  1:2286654/150‑1 (MQ=255)
                                                                                               aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                        >  2:3149552/1‑150 (MQ=255)
                                                                                               aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                        >  2:2645200/1‑150 (MQ=255)
                                                                                                    cgcgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAAc                                                   <  1:907344/150‑1 (MQ=255)
                                                                                                      cgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGt                                                 <  1:816061/150‑1 (MQ=255)
                                                                                                       ggATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTa                                                >  1:3571925/1‑150 (MQ=255)
                                                                                                           gggATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                            <  2:102286/150‑1 (MQ=255)
                                                                                                           gggATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                            <  2:4103844/150‑1 (MQ=255)
                                                                                                           gggATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTc                                            <  2:749587/150‑1 (MQ=255)
                                                                                                             gATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTCTCCGGCCTACGGGTGCCGAACGTAGGTCgg                                          >  2:2972835/1‑150 (MQ=255)
                                                                                                             gATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCgg                                          >  1:1918672/1‑150 (MQ=255)
                                                                                                              aTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGa                                         >  2:2334676/1‑150 (MQ=255)
                                                                                                              aTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGa                                         >  2:2464509/1‑150 (MQ=255)
                                                                                                              aTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGa                                         >  2:3765979/1‑150 (MQ=255)
                                                                                                                  aGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAg                                     <  1:3946994/150‑1 (MQ=255)
                                                                                                                  aGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAg                                     >  1:2250536/1‑150 (MQ=255)
                                                                                                                   gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                    >  1:4168116/1‑150 (MQ=255)
                                                                                                                   gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                    >  1:2449399/1‑150 (MQ=255)
                                                                                                                      aaGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                                 >  2:3090365/1‑150 (MQ=255)
                                                                                                                      aaGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGt                                 <  1:1092653/150‑1 (MQ=255)
                                                                                                                       aGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGtt                                >  2:482960/1‑150 (MQ=255)
                                                                                                                           tggtggTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGAAGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcg                            >  2:2748619/1‑150 (MQ=255)
                                                                                                                            ggtggtGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgc                           <  1:2672901/150‑1 (MQ=255)
                                                                                                                            ggtggtGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgc                           <  2:4055766/150‑1 (MQ=255)
                                                                                                                            ggtggtGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgc                           <  1:2025407/150‑1 (MQ=255)
                                                                                                                              tggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccg                         >  2:2158744/1‑150 (MQ=255)
                                                                                                                               ggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccgc                        <  1:1154885/150‑1 (MQ=255)
                                                                                                                               ggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcgccgc                        >  1:3539048/1‑150 (MQ=255)
                                                                                                                                  gtCGCACGTAGGTCAGATAAGGCGGTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGGGCGCATTGCCGGATGCGAGGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                     <  2:4169904/150‑1 (MQ=255)
                                                                                                                                     gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                  >  2:2837525/1‑150 (MQ=255)
                                                                                                                                     gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                  <  2:1496553/150‑1 (MQ=255)
                                                                                                                                     gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                  <  2:3794118/150‑1 (MQ=255)
                                                                                                                                     gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCGTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGCGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                  <  2:3640024/150‑1 (MQ=255)
                                                                                                                                      aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGAc                 >  1:767560/1‑150 (MQ=255)
                                                                                                                                         gTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt              <  1:1733398/150‑1 (MQ=255)
                                                                                                                                            tgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt           >  2:367978/2‑150 (MQ=38)
                                                                                                                                            tgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt           >  2:653646/2‑150 (MQ=38)
                                                                                                                                            tgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt           >  2:913131/2‑150 (MQ=38)
                                                                                                                                             gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg          >  2:3240535/1‑150 (MQ=255)
                                                                                                                                             gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg          >  2:1813508/1‑150 (MQ=255)
                                                                                                                                             gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg          >  2:2058292/1‑150 (MQ=255)
                                                                                                                                             gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg          >  2:814140/1‑150 (MQ=255)
                                                                                                                                              tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg         <  2:1170183/150‑1 (MQ=255)
                                                                                                                                              tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATgg         >  2:3924265/1‑150 (MQ=255)
                                                                                                                                               cAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc        <  1:1970751/150‑1 (MQ=255)
                                                                                                                                                aGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCa       >  1:3096397/1‑150 (MQ=255)
                                                                                                                                                    aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAGGGCGTTCACGCCGCATCCGACATTAATGGCACGtt   <  1:992532/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt   >  1:4046161/1‑150 (MQ=255)
                                                                                                                                                    aaGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtt   <  1:1397740/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGAAAAACTGGCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTT  >  NC_000913/3082553‑3082851

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