Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,395,009 G→A 5.8% A305V (GCA→GTA) ‡ ynaI ← mechanosensitive channel protein, very small conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,395,0090GA 5.8% 224.1 / 3.1 69A305V (GCA→GTA) ‡ynaImechanosensitive channel protein, very small conductance
Reads supporting (aligned to +/- strand):  ref base G (30/35);  new base A (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

CCTTGTGCATAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTC  >  NC_000913/1394862‑1395157
                                                                                                                                                   |                                                                                                                                                    
ccTTGTGCATAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGCGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTgc                                                                                                                                                     >  3:501054/1‑149 (MQ=255)
  ttGTGCATAACGCCAGTGATGTCTGATGTAGCGACCCTGTTCCGGCGGTGTAATGTTACTCATATACAGCGTCTGGCTGCGGAAGGCAAAGTCCGCGCCGTGTGACCGTACAATATCGATGATCCTCAAATAAACGTCTTGCTGTgcag                                                                                                                                                   >  4:51270/1‑149 (MQ=255)
     tgCATAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCaa                                                                                                                                                <  1:326263/149‑1 (MQ=255)
       cATAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGCACAGCAAga                                                                                                                                              <  2:290559/149‑2 (MQ=255)
        aTAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGcc                                                                                                                                             >  4:35866/1‑149 (MQ=255)
           aCGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCAtt                                                                                                                                          >  2:545581/1‑149 (MQ=255)
                       tCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCg                                                                                                                              >  4:24753/1‑149 (MQ=255)
                          gATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTgg                                                                                                                           >  4:285521/1‑149 (MQ=255)
                          gATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTgg                                                                                                                           >  2:399521/1‑149 (MQ=255)
                             ttAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGttt                                                                                                                        <  4:511200/149‑1 (MQ=255)
                               aGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTg                                                                                                                      <  3:410054/149‑1 (MQ=255)
                               aGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTg                                                                                                                      <  4:513216/149‑1 (MQ=255)
                               aGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTg                                                                                                                      <  2:46875/149‑1 (MQ=255)
                               aGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGAGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTg                                                                                                                      <  2:475870/149‑1 (MQ=255)
                                gCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTgg                                                                                                                     <  3:285521/149‑1 (MQ=255)
                                 cGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGt                                                                                                                    <  1:519208/149‑1 (MQ=255)
                                      cTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGc                                                                                                               >  3:303663/1‑149 (MQ=255)
                                        gTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCaa                                                                                                             <  3:24753/149‑1 (MQ=255)
                                           ccGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATaa                                                                                                          <  2:553972/149‑1 (MQ=255)
                                                      aTGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATAtt                                                                                               <  2:276554/149‑1 (MQ=255)
                                                       tGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGCTAAAGCAATAAACCATAATATTc                                                                                              >  4:407404/1‑149 (MQ=255)
                                                         ttATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCaa                                                                                            <  4:34144/149‑1 (MQ=255)
                                                               atatACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAaga                                                                                      >  3:246693/1‑149 (MQ=255)
                                                                 ataCAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGt                                                                                    <  4:303663/149‑1 (MQ=255)
                                                                      gCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCa                                                                               >  2:30142/1‑149 (MQ=255)
                                                                        gTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCaaa                                                                             >  4:488778/1‑149 (MQ=255)
                                                                           tggctggGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTg                                                                          <  3:582112/149‑1 (MQ=255)
                                                                            ggctggGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTgg                                                                         <  4:415849/149‑1 (MQ=255)
                                                                              ctggGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGtt                                                                       <  1:30142/149‑1 (MQ=255)
                                                                               tggGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTg                                                                      <  1:82688/149‑1 (MQ=255)
                                                                                ggggAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGa                                                                     <  1:457276/149‑1 (MQ=255)
                                                                                 gggAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGaa                                                                    <  1:498981/149‑1 (MQ=255)
                                                                                   gAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAAt                                                                  >  4:106191/1‑149 (MQ=255)
                                                                                      ggCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATaaa                                                               >  1:108099/1‑149 (MQ=255)
                                                                                         aaaGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAc                                                              >  3:295399/1‑147 (MQ=255)
                                                                                         aaaGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAAc                                                              <  4:295399/147‑1 (MQ=255)
                                                                                         aaaGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCa                                                            <  2:371661/149‑1 (MQ=255)
                                                                                         aaaGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCa                                                            >  4:567259/1‑149 (MQ=255)
                                                                                              ccGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAg                                                       <  2:399429/149‑1 (MQ=255)
                                                                                              ccGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAg                                                       <  3:567259/149‑1 (MQ=255)
                                                                                                gcgcCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGt                                                     <  3:263159/149‑1 (MQ=255)
                                                                                                  gcCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGttt                                                   >  4:216339/1‑149 (MQ=255)
                                                                                                         gACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTg                                            >  4:270024/1‑149 (MQ=255)
                                                                                                           cTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGt                                          >  2:13901/1‑149 (MQ=255)
                                                                                                            tGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTgg                                                                                                                     <  3:270603/73‑1 (MQ=255)
                                                                                                            tGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTgg                                                                                                                     >  4:270603/1‑73 (MQ=255)
                                                                                                            tGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAg                                                                                                                <  1:531244/78‑1 (MQ=255)
                                                                                                            tGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAg                                                                                                                >  2:531244/1‑78 (MQ=255)
                                                                                                                  atatCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGc                                   >  1:368872/1‑149 (MQ=255)
                                                                                                                   tatCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGcc                                  >  1:130700/1‑149 (MQ=255)
                                                                                                                     tCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCgg                                <  3:509546/149‑1 (MQ=255)
                                                                                                                       gATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGt                              <  2:547402/149‑1 (MQ=255)
                                                                                                                       gATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGt                              >  2:367336/1‑149 (MQ=255)
                                                                                                                              ttCAAATNANCGNCTTGCTGTGCAGNAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAttttt                       >  2:83660/1‑149 (MQ=255)
                                                                                                                              ttCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGtt                                                                       <  2:592861/101‑1 (MQ=255)
                                                                                                                              ttCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGtt                                                                       >  1:592861/1‑101 (MQ=255)
                                                                                                                              ttCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGAttttt                       <  1:46001/149‑1 (MQ=255)
                                                                                                                               tCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATAAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTc                      <  3:407404/149‑1 (MQ=255)
                                                                                                                                cAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAGACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCa                     <  3:359806/149‑1 (MQ=255)
                                                                                                                                 aaaTAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAg                    >  4:575716/1‑149 (MQ=255)
                                                                                                                                  aaTAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGc                   >  3:71156/1‑149 (MQ=255)
                                                                                                                                   aTAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCa                  >  4:535657/1‑149 (MQ=255)
                                                                                                                                    tgaatggcTTGCTGCACAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCAt                 >  4:438009/8‑149 (MQ=255)
                                                                                                                                        tggcTTGCTGCACAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTc                      <  3:559461/137‑1 (MQ=255)
                                                                                                                                        tggcTTGCTGCACAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTc                      >  4:559461/4‑140 (MQ=255)
                                                                                                                                        cGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCa             <  4:9022/149‑1 (MQ=255)
                                                                                                                                        cGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCa             <  3:535657/149‑1 (MQ=255)
                                                                                                                                         gTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCAc            >  2:36761/1‑149 (MQ=255)
                                                                                                                                         gTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCAc            >  3:353442/1‑149 (MQ=255)
                                                                                                                                          tCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGc                                              >  3:235/1‑114 (MQ=255)
                                                                                                                                          tCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGc                                              <  4:235/114‑1 (MQ=255)
                                                                                                                                                   gcagcaAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTc  <  1:568438/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTTGTGCATAACGCCAGTGATGTCTGATTTAGCGACCCTGTTCCGGCGGTGTAATGTTATCCATATACAGCGTCTGGCTGGGGAAGGCAAAGTCCGCGCCGTGTGACTGTACAATATCGATAATCTTCAAATAAACGTCTTGCTGTGCAGCAAGCCATTCAGCCCATACCGTGGTTTTGGTAAAGCAATAAACCATAATATTCAATGAAGAGTCAGCAAACTGGTTGAAATAAACCAGTAAGGTTTGTCGCTGGTCGATGGCCGGGTGATTTTTCAGCATCTCACGTACAGCTTC  >  NC_000913/1394862‑1395157

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