Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,716 C→T 6.4% intergenic (+45/+161) 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,7160CT 6.4% 290.4 / 4.8 94intergenic (+45/+161)loiP/speBPhe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase
Reads supporting (aligned to +/- strand):  ref base C (43/43);  new base T (2/4);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCAT  >  NC_000913/3082573‑3082864
                                                                                                                                               |                                                                                                                                                    
gaagaaGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCt                                                                                                                                                  <  2:212749/148‑1 (MQ=255)
 aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGa                                                                                                                                                >  3:144691/1‑149 (MQ=255)
  agaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                               <  2:366033/149‑1 (MQ=255)
  agaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                               <  1:267447/149‑1 (MQ=255)
   gaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGAcc                                                                                                                                              <  2:141925/149‑2 (MQ=255)
    aagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTTGGTCAGATAAGGCGTTCACGACGCATCTGACAt                                                                                                                                             <  8:238971/149‑1 (MQ=255)
    aagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACAt                                                                                                                                             <  5:288706/149‑1 (MQ=255)
       gTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                              >  1:315615/1‑145 (MQ=255)
       gTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                              <  2:315615/145‑1 (MQ=255)
       gTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGaa                                                                                                                                          <  2:343385/149‑1 (MQ=255)
        tCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAAt                                                                                                                                         <  6:319471/149‑1 (MQ=255)
        tCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAAt                                                                                                                                         <  6:18221/149‑1 (MQ=255)
          gCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                                       >  5:83652/1‑149 (MQ=255)
              aaGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGcaca                                                                                                                                   <  2:223539/149‑1 (MQ=255)
                 cTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAttt                                                                                                                                <  7:7270/149‑1 (MQ=255)
                 cTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAttt                                                                                                                                <  1:173890/149‑1 (MQ=255)
                  tCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTg                                                                                                                               >  1:104724/1‑149 (MQ=255)
                  tCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTg                                                                                                                               >  6:163156/1‑149 (MQ=255)
                   cAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGt                                                                                                                              >  1:73131/1‑149 (MQ=255)
                       gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                              >  1:256152/1‑129 (MQ=255)
                       gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                              <  2:256152/129‑1 (MQ=255)
                        tttGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTAcc                                                                                                                          <  8:249059/148‑1 (MQ=255)
                          tGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCtt                                                                                                                        <  2:73131/148‑1 (MQ=255)
                              gacCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTgcg                                                                                                                   >  3:166484/1‑149 (MQ=255)
                                 cATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCAt                                                                                                                >  5:158112/1‑149 (MQ=255)
                                    ccTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGc                                                                                                             >  1:447170/1‑149 (MQ=255)
                                        cATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                                         <  6:158112/149‑1 (MQ=255)
                                         atgcggcgAAAGGGCCCATGATATTCGCGTTCGGATGTGCGGGGTTGGGTTTTATTTAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGga                                                                                                         <  8:425897/142‑1 (MQ=25)
                                         aTCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgat                                                                                                        >  7:145776/1‑149 (MQ=255)
                                         aTCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgat                                                                                                        >  5:30440/1‑149 (MQ=255)
                                          tCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatg                                                                                                       <  1:85787/149‑1 (MQ=255)
                                             gccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATgg                                                                                                                                       <  8:375197/114‑1 (MQ=255)
                                             accgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcga                                                                                                    <  7:210336/148‑1 (MQ=255)
                                              ccgAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGAATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcgat                                                                                                   >  7:453976/1‑149 (MQ=255)
                                                   cgcgCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTgg                                                                                              <  8:359124/149‑1 (MQ=255)
                                                    gcgcCCAGAATATTCGCGATCGCATGAGAGAGGATGGGATAAAGTAAAGCCTGGTGGGGTCGCACGTAGGTCAGATAAGGCGTGCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCGGgc                                                                                             <  5:20123/149‑1 (MQ=255)
                                                      gcCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgcgc                                                                                           <  2:82285/149‑1 (MQ=255)
                                                       cccAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCa                                                                                          <  4:144691/149‑1 (MQ=255)
                                                         cAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                                        >  6:387705/1‑149 (MQ=255)
                                                         cAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                                        >  3:424008/1‑149 (MQ=255)
                                                          aGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgtg                                                                                                                     >  4:273366/1‑119 (MQ=255)
                                                          aGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgtg                                                                                                                     <  3:273366/119‑1 (MQ=255)
                                                              tatTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATc                                                                                   >  2:125134/1‑149 (MQ=255)
                                                                ttCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCg                                                                                 >  3:120564/1‑149 (MQ=255)
                                                                  cgcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCgg                                                                                <  4:278780/148‑1 (MQ=255)
                                                                   gcgATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGcc                                                                              >  8:308522/1‑149 (MQ=255)
                                                                     gATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTa                                                                            >  7:464213/1‑149 (MQ=255)
                                                                       tCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcga                                                                                                    <  8:79478/123‑1 (MQ=255)
                                                                       tCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcga                                                                                                    >  7:79478/1‑123 (MQ=255)
                                                                         gCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACggg                                                                        >  7:63325/1‑149 (MQ=255)
                                                                         gCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACgg                                                                         <  6:83652/148‑1 (MQ=255)
                                                                          cATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGt                                                                       >  8:371678/1‑149 (MQ=255)
                                                                           aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                                        >  7:96352/1‑131 (MQ=255)
                                                                           aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                                        <  8:96352/131‑1 (MQ=255)
                                                                           aTGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTg                                                                      >  4:361435/1‑149 (MQ=255)
                                                                                  cgGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGCAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACg                                                               <  8:376156/149‑1 (MQ=255)
                                                                                   ggATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTCACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCTGATGCGATGCTGGCGCATCTTATCCTGCCTACTGGTGCCGAACg                                                               >  8:411090/1‑148 (MQ=255)
                                                                                    gATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGAGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACg                                                               >  1:116157/1‑147 (MQ=255)
                                                                                          aTTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCg                                                        >  7:82506/1‑148 (MQ=255)
                                                                                           ttAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGa                                                      <  5:237374/149‑1 (MQ=255)
                                                                                             aaGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                                     <  1:430899/132‑1 (MQ=255)
                                                                                             aaGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                                     >  2:430899/1‑132 (MQ=255)
                                                                                              aGTAAAGCCTCGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGAACATTAGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTACGTCGGAt                                                     >  5:158555/1‑147 (MQ=255)
                                                                                               gTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAg                                                  >  4:412496/1‑149 (MQ=255)
                                                                                                tAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAgg                                                 >  2:351835/1‑149 (MQ=255)
                                                                                                    gCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGtt                                             <  3:361435/149‑1 (MQ=255)
                                                                                                     ccTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGtt                                             >  2:435425/1‑148 (MQ=255)
                                                                                                     ccTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTc                                            <  6:30440/149‑1 (MQ=255)
                                                                                                         gtggtgTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCAcg                                         <  5:387705/148‑1 (MQ=255)
                                                                                                              gtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCAt                                   >  1:247624/1‑149 (MQ=255)
                                                                                                              gtCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCAt                                   >  6:320393/1‑149 (MQ=255)
                                                                                                               tCGAACGTATGTCAGATAAGGCGTTCACGCCGCATCTGACAGGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCGACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                                  <  7:411090/149‑1 (MQ=255)
                                                                                                               tCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATGCGGCGTGAAcgccttat                                   <  7:371678/148‑5 (MQ=38)
                                                                                                               tCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATc                                  >  1:184176/1‑149 (MQ=255)
                                                                                                                cGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                 <  1:125134/149‑1 (MQ=255)
                                                                                                                 gAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATcc                                 >  3:201895/1‑148 (MQ=255)
                                                                                                                  aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAtt                                                                                                                                 >  7:220676/1‑51 (MQ=255)
                                                                                                                  aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAtt                                                                                                                                 <  8:220676/51‑1 (MQ=255)
                                                                                                                   aCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                               >  1:155767/1‑148 (MQ=255)
                                                                                                                       atgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt                           >  5:377765/3‑148 (MQ=38)
                                                                                                                       atgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt                           >  3:429402/3‑148 (MQ=38)
                                                                                                                         gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                        <  2:155767/149‑1 (MQ=255)
                                                                                                                         gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt                        <  5:163156/149‑1 (MQ=255)
                                                                                                                          tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg                       >  1:324599/1‑149 (MQ=255)
                                                                                                                            aGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                                           <  2:440462/95‑1 (MQ=255)
                                                                                                                            aGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                                           >  1:440462/1‑95 (MQ=255)
                                                                                                                            aGATAAGGCGTTCACGCCGAATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc                     <  8:453976/149‑1 (MQ=255)
                                                                                                                             gATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgc                                                                                             <  3:87184/76‑1 (MQ=255)
                                                                                                                             gATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgc                                                                                             >  4:87184/1‑76 (MQ=255)
                                                                                                                             gATAAGGCGTTCACGCCGCACCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCa                    >  8:50963/1‑149 (MQ=255)
                                                                                                                                  ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGttt               >  2:420218/1‑149 (MQ=255)
                                                                                                                                   gCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt              <  2:447170/149‑1 (MQ=255)
                                                                                                                                     gtgaaCGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTAc            <  4:429402/145‑1 (MQ=255)
                                                                                                                                     gtgaaCGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTAc            <  6:377765/145‑1 (MQ=255)
                                                                                                                                           gccgcATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTgc      >  8:388044/1‑149 (MQ=255)
                                                                                                                                               cATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGcgcat  <  2:104724/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTTACCCGTGCGCAT  >  NC_000913/3082573‑3082864

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: ‑