Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,387,774 G→T 5.1% P533T (CCG→ACG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,387,7740GT 5.1% 346.1 / 7.0 99P533T (CCG→ACG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (43/51);  new base T (2/3);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

CAAGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATC  >  NC_000913/4387628‑4387922
                                                                                                                                                  |                                                                                                                                                    
cAAGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTa                                                                                                                                                    >  4:124448/1‑149 (MQ=255)
  aGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAAc                                                                                                                                                  >  1:782378/1‑149 (MQ=255)
  aGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAAc                                                                                                                                                  >  3:690617/1‑149 (MQ=255)
   gCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACg                                                                                                                                                 >  2:388137/1‑149 (MQ=255)
      gTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTg                                                                                                                                              <  3:22828/149‑1 (MQ=255)
        tCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGAc                                                                                                                                            <  4:713076/149‑1 (MQ=255)
          gggcgggcAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGTTACCGGTgctgt                                                                                                                                          <  3:471514/149‑5 (MQ=255)
               ggcAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGcgacga                                                                                                                                     <  3:71791/149‑1 (MQ=255)
                gcAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACgag                                                                                                                                    >  2:670245/1‑149 (MQ=255)
                 cAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACgaga                                                                                                                                   >  3:106555/1‑149 (MQ=255)
                  aaaGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGcgacga                                                                                                                                     <  4:622229/146‑1 (MQ=255)
                  aaaGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGcgacga                                                                                                                                     >  3:622229/1‑146 (MQ=255)
                   aaGGTGGCGCAAGTCAACACCCCCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGTCCGCCAGCGGATAAGGCCACGCCTCGCGCAAGACGTAACCCATCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAAt                                                                                                                                 <  2:831267/149‑1 (MQ=255)
                   aaGGTGGCGCAAATCATCACTACCCACAGCAGCGGTACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCCTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAAt                                                                                                                                 <  3:602921/149‑1 (MQ=255)
                   aaGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAAt                                                                                                                                 <  2:848077/149‑1 (MQ=255)
                   aaGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAAt                                                                                                                                 <  1:388137/149‑1 (MQ=255)
                      gTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGa                                                                                                                              <  4:44395/149‑1 (MQ=255)
                       tGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGAc                                                                                                                             <  2:351258/149‑1 (MQ=255)
                                tcatcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAaga                                                                                                                    <  3:544006/149‑1 (MQ=255)
                                tcatcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAaga                                                                                                                    <  1:873919/149‑1 (MQ=255)
                                tcatcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAaga                                                                                                                    <  3:324032/149‑1 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGTTACCGGTgctgt                                                                                                                                          <  4:411224/125‑5 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGTTACCGGTgctgt                                                                                                                                          >  3:411224/1‑121 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGt                                                                                                                  <  2:782378/149‑1 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGt                                                                                                                  >  2:880981/1‑149 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGt                                                                                                                  >  3:731101/1‑149 (MQ=255)
                                  atcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGt                                                                                                                  >  4:858716/1‑149 (MQ=255)
                                       tACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCa                                                                                                             >  1:288305/1‑149 (MQ=255)
                                          ccACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACg                                                                                                          >  3:102433/1‑149 (MQ=255)
                                          ccACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACg                                                                                                          <  3:166619/149‑1 (MQ=255)
                                             cagcagCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCa                                                                                                       <  2:216861/149‑1 (MQ=255)
                                                    ggCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGaa                                                                                                <  4:106555/149‑1 (MQ=255)
                                                     gCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAg                                                                                               >  3:288936/1‑149 (MQ=255)
                                                       aCGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTg                                                                                             >  4:42125/1‑149 (MQ=255)
                                                        cGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                                                                            <  1:700767/149‑1 (MQ=255)
                                                           tGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATcc                                                                                         <  1:328441/149‑1 (MQ=255)
                                                            ggCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCt                                                                                        <  4:827842/149‑1 (MQ=255)
                                                            ggCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCt                                                                                        <  4:258882/149‑1 (MQ=255)
                                                             gCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTg                                                                                       >  4:226759/1‑149 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  2:31703/149‑1 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  4:690617/149‑1 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      >  2:737080/1‑149 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  1:66798/149‑1 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  3:456322/149‑1 (MQ=255)
                                                              ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  4:282189/149‑1 (MQ=255)
                                                                   acacCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCAc                                                                                 <  2:406444/149‑1 (MQ=255)
                                                                     acCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACtt                                                                               <  2:42120/149‑1 (MQ=255)
                                                                       cATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTg                                                                             >  3:793764/1‑149 (MQ=255)
                                                                        aTCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCTCAACGTCAGCCAGAAGTGATACTGAGTCACTTTGc                                                                            >  1:173397/1‑149 (MQ=255)
                                                                        aTCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      >  4:701039/1‑139 (MQ=255)
                                                                        aTCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  3:701039/139‑1 (MQ=255)
                                                                        aTCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGc                                                                            >  2:795138/1‑149 (MQ=255)
                                                                            ccAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACTCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAAc                                                                        <  3:553224/149‑1 (MQ=255)
                                                                            ccAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAAc                                                                        >  2:527221/1‑149 (MQ=255)
                                                                               aTCGCGACCGCCAGCGGATAAGGCCACGACTCGCGCAAGCCGTGACCCAGCGTCATCCACAGCACCGGTGACGGTGATGCGACGAGAATCGACCAGACAAGAGTGCGCAACGTCAGCCAGAAGCGATCCTGAGTCACTTTGCCAACttt                                                                     <  4:548698/149‑1 (MQ=255)
                                                                                tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTg                                                                                                                 >  4:422052/1‑104 (MQ=255)
                                                                                tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTg                                                                                                                 <  3:422052/104‑1 (MQ=255)
                                                                                tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTc                                                                    <  4:753286/149‑1 (MQ=255)
                                                                                tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTc                                                                    <  2:838406/149‑1 (MQ=255)
                                                                                 cgcgACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCTTCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTcg                                                                   >  3:76888/1‑149 (MQ=255)
                                                                                   cgACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTcgcc                                                                 >  3:303731/1‑149 (MQ=255)
                                                                                   cgACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTcgcc                                                                 <  4:731101/149‑1 (MQ=255)
                                                                                                aTAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCaa                                                    >  3:810036/1‑149 (MQ=255)
                                                                                                    ggCCACGCCTCGCGCAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACg                                                                                                                                                 >  2:211582/1‑51 (MQ=255)
                                                                                                    ggCCACGCCTCGCGCAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACg                                                                                                                                                 <  1:211582/51‑1 (MQ=255)
                                                                                                     gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      >  3:501316/1‑110 (MQ=255)
                                                                                                     gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  4:501316/110‑1 (MQ=255)
                                                                                                      ccACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAAc                                              >  1:161792/1‑149 (MQ=255)
                                                                                                      ccACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAAc                                              <  2:813624/149‑1 (MQ=255)
                                                                                                      ccACGCCTCGCGCAAGCCGTAACCCAGCGACAGCCACAGCAACGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTACGTCCCTTTGCCAACTTTCGCGGCCGAACGTTCAACAGAAc                                              >  2:63264/1‑149 (MQ=255)
                                                                                                            cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                                                                            <  2:466864/97‑1 (MQ=255)
                                                                                                            cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                                                                            >  1:466864/1‑97 (MQ=255)
                                                                                                            cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGAt                                                                                           <  2:91563/98‑1 (MQ=255)
                                                                                                            cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGAt                                                                                           >  1:91563/1‑98 (MQ=255)
                                                                                                            cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTAACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGa                                        >  2:760603/1‑149 (MQ=255)
                                                                                                             tCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGaa                                       >  2:688121/1‑149 (MQ=255)
                                                                                                             tCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGaa                                       >  1:757843/1‑149 (MQ=255)
                                                                                                                  cAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgc                                  <  2:496330/149‑1 (MQ=255)
                                                                                                                  cAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgc                                  <  2:281741/149‑1 (MQ=255)
                                                                                                                    aGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgcgg                                <  2:846320/149‑1 (MQ=255)
                                                                                                                     gCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgcggc                               <  2:23558/149‑1 (MQ=255)
                                                                                                                        gTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      >  2:696992/1‑91 (MQ=255)
                                                                                                                        gTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  1:696992/91‑1 (MQ=255)
                                                                                                                        gTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCAGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGcgcg                            >  4:607697/1‑149 (MQ=255)
                                                                                                                             ccAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTaa                       >  2:161977/1‑149 (MQ=255)
                                                                                                                              cAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTaaa                      <  1:561643/149‑1 (MQ=255)
                                                                                                                              cAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTaaa                      >  2:561643/1‑149 (MQ=255)
                                                                                                                               aGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      <  3:622711/84‑1 (MQ=255)
                                                                                                                               aGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                      >  4:622711/1‑84 (MQ=255)
                                                                                                                               aGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAAt                     >  2:845754/1‑149 (MQ=255)
                                                                                                                               aGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAAt                     <  2:866325/149‑1 (MQ=255)
                                                                                                                                gCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATa                    <  3:222751/149‑1 (MQ=255)
                                                                                                                                  gTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACt                  >  1:41693/1‑149 (MQ=255)
                                                                                                                                  gTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACt                  <  1:29937/149‑1 (MQ=255)
                                                                                                                                      tCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCaa              >  2:370875/1‑149 (MQ=255)
                                                                                                                                          cgtcACCGTTACCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGt                         >  3:771621/4‑134 (MQ=255)
                                                                                                                                          cgtcACCGTTACCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGt                         <  4:771621/131‑1 (MQ=255)
                                                                                                                                               ccGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGa     >  1:363856/1‑149 (MQ=255)
                                                                                                                                               ccGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGa     <  2:402670/149‑1 (MQ=255)
                                                                                                                                                  gTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATc  >  3:104196/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAAGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATC  >  NC_000913/4387628‑4387922

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