Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,387,777 A→C 5.1% L532V (TTA→GTA)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,387,7770AC 5.1% 343.2 / 9.6 99L532V (TTA→GTA) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base A (43/51);  new base C (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 = 1.92e-01

Only 100 of 101 total aligned reads displayed.

AGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAAC  >  NC_000913/4387630‑4387925
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                       <  1:873919/149‑1 (MQ=255)
                              tcatcaCTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAaga                                                                                                                       <  3:544006/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)
                                     tACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCa                                                                                                                >  1:288305/1‑149 (MQ=255)
                                        ccACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACg                                                                                                             <  3:166619/149‑1 (MQ=255)
                                        ccACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACg                                                                                                             >  3:102433/1‑149 (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                                                                                         <  1:66798/149‑1 (MQ=255)
                                                            ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         >  2:737080/1‑149 (MQ=255)
                                                            ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  4:282189/149‑1 (MQ=255)
                                                            ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  3:456322/149‑1 (MQ=255)
                                                            ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  2:31703/149‑1 (MQ=255)
                                                            ccGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  4:690617/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                                                                                         <  3:701039/139‑1 (MQ=255)
                                                                      aTCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         >  4:701039/1‑139 (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                                                                                                                    <  3:422052/104‑1 (MQ=255)
                                                                              tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTg                                                                                                                    >  4:422052/1‑104 (MQ=255)
                                                                              tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTc                                                                       <  2:838406/149‑1 (MQ=255)
                                                                              tCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTc                                                                       <  4:753286/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                                                                                                                                                    <  1:211582/51‑1 (MQ=255)
                                                                                                  ggCCACGCCTCGCGCAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACg                                                                                                                                                    >  2:211582/1‑51 (MQ=255)
                                                                                                   gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  4:501316/110‑1 (MQ=255)
                                                                                                   gCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         >  3:501316/1‑110 (MQ=255)
                                                                                                    ccACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAAc                                                 <  2:813624/149‑1 (MQ=255)
                                                                                                    ccACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAAc                                                 >  1:161792/1‑149 (MQ=255)
                                                                                                    ccACGCCTCGCGCAAGCCGTAACCCAGCGACAGCCACAGCAACGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTACGTCCCTTTGCCAACTTTCGCGGCCGAACGTTCAACAGAAc                                                 >  2:63264/1‑149 (MQ=255)
                                                                                                          cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                                                                               >  1:466864/1‑97 (MQ=255)
                                                                                                          cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGa                                                                                               <  2:466864/97‑1 (MQ=255)
                                                                                                          cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGAt                                                                                              >  1:91563/1‑98 (MQ=255)
                                                                                                          cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGAt                                                                                              <  2:91563/98‑1 (MQ=255)
                                                                                                          cTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTAACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGa                                           >  2:760603/1‑149 (MQ=255)
                                                                                                           tCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGaa                                          >  1:757843/1‑149 (MQ=255)
                                                                                                           tCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGaa                                          >  2:688121/1‑149 (MQ=255)
                                                                                                                cAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgc                                     <  2:281741/149‑1 (MQ=255)
                                                                                                                cAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgc                                     <  2:496330/149‑1 (MQ=255)
                                                                                                                  aGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgcgg                                   <  2:846320/149‑1 (MQ=255)
                                                                                                                   gCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAgcggc                                  <  2:23558/149‑1 (MQ=255)
                                                                                                                      gTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         <  1:696992/91‑1 (MQ=255)
                                                                                                                      gTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGa                                                                                         >  2:696992/1‑91 (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:866325/149‑1 (MQ=255)
                                                                                                                             aGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAAt                        >  2:845754/1‑149 (MQ=255)
                                                                                                                              gCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATa                       <  3:222751/149‑1 (MQ=255)
                                                                                                                                gTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACt                     <  1:29937/149‑1 (MQ=255)
                                                                                                                                gTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACt                     >  1:41693/1‑149 (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)
                                                                                                                                                  aaCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCaa   <  4:236999/149‑1 (MQ=255)
                                                                                                                                                   aCGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAAc  <  1:894679/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCCGTTCGGGCGGGCAAAGGTGGCGCAAATCATCACTACCCACAGCAGCGGCACGGTGGCCGTTACACCATCGCCAATCGCGACCGCCAGCGGATAAGGCCACGCCTCGCGCAAGCCGTAACCCAGCGTCATCCACAGCACCGGTAACGGTGACGCGACGAGAATCGACCAGAAAAGAGTGCGCAACGTCAGCCAGAAGTGATCCTGAGTCACTTTGCCAACTTTCGCCGCCGAACGTTCAAGAAAACGGGTGAAATAGCGGCGCGAGTAAATACTGCAACCGACCAGAATCAAC  >  NC_000913/4387630‑4387925

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