Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,082,708 T→G 5.0% intergenic (+37/+169) 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,7080TG 5.0% 287.0 / 5.7 94intergenic (+37/+169)loiP/speBPhe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; Era‑binding heat shock protein/agmatinase
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base G (2/3);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTT  >  NC_000913/3082564‑3082852
                                                                                                                                                |                                                                                                                                                
gCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcg                                                                                                                                              >  6:237374/1‑149 (MQ=255)
    aaCTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCACCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCAcgccg                                                                                                                                           >  8:210336/1‑148 (MQ=255)
         gaagaaGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCt                                                                                                                                      <  2:212749/148‑1 (MQ=255)
          aagaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGa                                                                                                                                    >  3:144691/1‑149 (MQ=255)
           agaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                   <  1:267447/149‑1 (MQ=255)
           agaagGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATGCGGCGTGAACGCCTTATCTGAc                                                                                                                                   <  2:366033/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                                                                                                                                  <  2:256152/129‑1 (MQ=255)
                                gTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACa                                                                                                                                  >  1:256152/1‑129 (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                                                                                            >  5:30440/1‑149 (MQ=255)
                                                  aTCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgat                                                                                            >  7:145776/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                                                                            >  3:424008/1‑149 (MQ=255)
                                                                  cAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATc                                                                            >  6:387705/1‑149 (MQ=255)
                                                                   aGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgtg                                                                                                         <  3:273366/119‑1 (MQ=255)
                                                                   aGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTtgtg                                                                                                         >  4:273366/1‑119 (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                                                                                        >  7:79478/1‑123 (MQ=255)
                                                                                tCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGgatgcga                                                                                        <  8:79478/123‑1 (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                                                         >  2:430899/1‑132 (MQ=255)
                                                                                                      aaGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGc                                                         <  1:430899/132‑1 (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                                                                                                                     <  8:220676/51‑1 (MQ=255)
                                                                                                                           aaCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACAtt                                                                                                                     >  7:220676/1‑51 (MQ=255)
                                                                                                                            aCGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGa                   >  1:155767/1‑148 (MQ=255)
                                                                                                                                atgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt               >  5:377765/3‑148 (MQ=38)
                                                                                                                                atgTCAGATGCGGCGTGAACGCCTTATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACAtt               >  3:429402/3‑148 (MQ=38)
                                                                                                                                  gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt            <  5:163156/149‑1 (MQ=255)
                                                                                                                                  gTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAAt            <  2:155767/149‑1 (MQ=255)
                                                                                                                                   tCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATg           >  1:324599/1‑149 (MQ=255)
                                                                                                                                     aGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                               >  1:440462/1‑95 (MQ=255)
                                                                                                                                     aGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTAc                                                               <  2:440462/95‑1 (MQ=255)
                                                                                                                                     aGATAAGGCGTTCACGCCGAATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGc         <  8:453976/149‑1 (MQ=255)
                                                                                                                                      gATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgc                                                                                 >  4:87184/1‑76 (MQ=255)
                                                                                                                                      gATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGgc                                                                                 <  3:87184/76‑1 (MQ=255)
                                                                                                                                      gATAAGGCGTTCACGCCGCACCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCa        >  8:50963/1‑149 (MQ=255)
                                                                                                                                           ggCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGttt   >  2:420218/1‑149 (MQ=255)
                                                                                                                                            gCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGtttt  <  2:447170/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCAAAACTGGAAGAAGGTCGCCAAAGCTCAATGTTTGACGACCATCCTGCATCCGCCGAACGCGCCCAGCATATTCGCGATCGCATGAGCGCGGATGGGATTAAGTAAAGCCTGGTGGTGTCGAACGTAGGTCAGATAAGGCGTTCACGCCGCATCTGACATGAATGGCACATTTGTTACCTTGTGCGCATTGCCGGATGCGATGCTGGCGCATCTTATCCGGCCTACGGGTGCCGAACGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTAATGGCACGTTTT  >  NC_000913/3082564‑3082852

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