Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,388,508 T→A 6.2% Q288L (CAG→CTG)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,5080TA 6.2% 308.9 / 8.1 97Q288L (CAG→CTG) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base T (48/43);  new base A (4/2);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

ATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCG  >  NC_000913/4388360‑4388656
                                                                                                                                                    |                                                                                                                                                    
aTCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCt                                                                                                                                                      >  4:192251/1‑149 (MQ=255)
 tCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCt                                                                                                                                                      <  2:340499/148‑1 (MQ=255)
   acgCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                                   <  1:224553/147‑1 (MQ=255)
   aaGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGATCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                                    <  7:350599/148‑1 (MQ=255)
    aGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACg                                                                                                                                                  >  5:81662/1‑149 (MQ=255)
         ttgtGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                                             <  1:204008/149‑1 (MQ=255)
           gtgGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTg                                                                                                                                            <  5:32194/148‑1 (MQ=255)
           gtgGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGc                                                                                                                                           >  5:414263/1‑149 (MQ=255)
             ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGa                                                                                                                                         <  1:445361/149‑1 (MQ=255)
              gTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAg                                                                                                                                        >  7:433718/1‑149 (MQ=255)
              gTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAg                                                                                                                                        <  2:249679/149‑1 (MQ=255)
               tttCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAgg                                                                                                                                       >  4:454922/1‑149 (MQ=255)
                   ggCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                                    <  1:69294/148‑1 (MQ=255)
                   ggCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                                   >  8:169486/1‑149 (MQ=255)
                    gCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGACgt                                                                                                                                                 <  1:96240/134‑2 (MQ=255)
                    gCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGACgt                                                                                                                                                 >  2:96240/1‑133 (MQ=255)
                     cATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                                   >  7:58048/1‑147 (MQ=255)
                      aTTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGGGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAAc                                                                                                                                   >  1:155156/1‑146 (MQ=255)
                                 ccGGGCCCCCTGCGCCCGCAGCGCGTCGCCGAGCAGATTGGACGCGCCCCGCCATTGCGACTGTTCACGCAGCGTCTGCCACGCCTGCCGGACCTGTAACGTCTTGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCt                                                                                                                     <  7:21238/149‑1 (MQ=255)
                                   gTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                                                                                                   >  3:206456/1‑149 (MQ=255)
                                      ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCg                                                                                                                       >  3:86695/1‑142 (MQ=255)
                                      ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCg                                                                                                                       <  4:86695/142‑1 (MQ=255)
                                      ccACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGcc                                                                                                                >  1:116225/1‑149 (MQ=255)
                                        aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                                                                                             <  7:232760/134‑1 (MQ=255)
                                        aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatc                                                                                                                             >  8:232760/1‑134 (MQ=255)
                                        aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTg                                                                                                              >  6:388050/1‑149 (MQ=255)
                                        aCCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTg                                                                                                              >  5:109157/1‑149 (MQ=255)
                                         ccTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                                             <  6:81662/149‑1 (MQ=255)
                                         ccTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                                             <  5:367205/149‑1 (MQ=255)
                                          cTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGtt                                                                                                            <  5:388050/149‑1 (MQ=255)
                                             cgcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGa                                                                                                          <  8:433718/148‑1 (MQ=255)
                                             cgcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAt                                                                                                         <  8:318716/149‑1 (MQ=255)
                                              gcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                                        <  2:152758/149‑1 (MQ=255)
                                                ccGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTc                                                                                                       >  2:379415/1‑148 (MQ=255)
                                                  gCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                                    <  4:206212/149‑1 (MQ=255)
                                                       gcTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                               <  2:116225/149‑1 (MQ=255)
                                                                     aTTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGt                                                                                 <  1:39600/149‑1 (MQ=255)
                                                                         gACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                                             <  7:332815/85‑1 (MQ=255)
                                                                         gACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGtt                                                                                                                                             >  8:332815/1‑85 (MQ=255)
                                                                          aCGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAtt                                                                            <  6:414263/149‑1 (MQ=255)
                                                                                ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGaa                                                                       <  2:334045/148‑1 (MQ=255)
                                                                                ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGaa                                                                       <  2:336146/148‑1 (MQ=255)
                                                                                ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAt                                                                      >  1:207676/1‑149 (MQ=255)
                                                                                ccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAt                                                                      >  1:5685/1‑149 (MQ=255)
                                                                                 cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                                    <  4:152096/118‑1 (MQ=255)
                                                                                 cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaaa                                                                                                    >  3:152096/1‑118 (MQ=255)
                                                                                 cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                                     >  3:69864/1‑149 (MQ=255)
                                                                                 cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                                     >  7:505999/1‑149 (MQ=255)
                                                                                  aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                                    >  7:18789/1‑149 (MQ=255)
                                                                                   gCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                                                   <  2:207676/149‑1 (MQ=255)
                                                                                    ccATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgcg                                                                  >  7:356142/1‑149 (MQ=255)
                                                                                       ttGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                               <  8:58048/149‑1 (MQ=255)
                                                                                       ttGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                               >  2:183177/1‑149 (MQ=255)
                                                                                         gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                                        >  5:189894/1‑122 (MQ=255)
                                                                                         gCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                                        <  6:189894/122‑1 (MQ=255)
                                                                                          cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                                   <  1:166054/62‑1 (MQ=255)
                                                                                          cGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGAc                                                                                                                                                   >  2:166054/1‑62 (MQ=255)
                                                                                           gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATa                                                                                          <  2:497724/118‑1 (MQ=255)
                                                                                           gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATa                                                                                          >  1:497724/1‑118 (MQ=255)
                                                                                           gACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcgacga                                                            >  2:275655/1‑148 (MQ=255)
                                                                                                  cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                                     <  7:169486/148‑1 (MQ=255)
                                                                                                  cACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                                     <  4:69864/148‑1 (MQ=255)
                                                                                                   aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                               <  8:348814/105‑1 (MQ=255)
                                                                                                   aCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccgc                                                                                               >  7:348814/1‑105 (MQ=255)
                                                                                                    cGCAGCGTATTCAACGCCTGACGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCgg                                                  <  6:109157/149‑1 (MQ=255)
                                                                                                        gCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCaaa                                              >  4:404015/1‑149 (MQ=255)
                                                                                                            aTTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                                          >  5:167042/1‑149 (MQ=255)
                                                                                                             ttCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGGCGATAGTTCGCGGTTAATTATGAATTGCGCGACGATATCTTTCGGCGAATCGgc                                         <  2:155156/149‑1 (MQ=255)
                                                                                                                aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGgcgc                                       <  4:206456/148‑1 (MQ=255)
                                                                                                                aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                                      >  7:457431/1‑149 (MQ=255)
                                                                                                                aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                                      >  3:441936/1‑149 (MQ=255)
                                                                                                                aaCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                                      >  4:421597/1‑149 (MQ=255)
                                                                                                                  cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGt                                    <  2:287463/149‑1 (MQ=255)
                                                                                                                    ccggCCGGACCTGTAACGTCTGGCTGGCCGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATGCAAAGCCGCCGATAGTTAGCGCTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGttt                                  >  8:352889/4‑149 (MQ=255)
                                                                                                                        ccGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCt                               <  4:441936/148‑1 (MQ=255)
                                                                                                                        ccGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCt                               <  1:379415/148‑1 (MQ=255)
                                                                                                                          ggACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCt                                                       >  7:455888/1‑122 (MQ=255)
                                                                                                                               aacagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGAGAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCaa                       >  1:440264/6‑149 (MQ=255)
                                                                                                                                gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt                      >  4:159436/1‑149 (MQ=255)
                                                                                                                                gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt                      >  3:35240/1‑149 (MQ=255)
                                                                                                                                gTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt                      >  8:189081/1‑149 (MQ=255)
                                                                                                                                acagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                               >  3:10500/5‑108 (MQ=255)
                                                                                                                                acagcGTCAGGCTGCCAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGcga                                                               <  4:10500/104‑1 (MQ=255)
                                                                                                                                 tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCACCGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAGTTTGCATTGCGCGACGATATCGTTCGGCAAATCGGCGGTGGTTTCTGCCAGCACtt                     >  6:283711/1‑149 (MQ=255)
                                                                                                                                 tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCGGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGGAAGCAAta                     >  4:259185/1‑148 (MQ=255)
                                                                                                                                 tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAtt                     >  6:82959/1‑149 (MQ=255)
                                                                                                                                 tAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAGTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCTCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAtt                     >  6:276794/1‑149 (MQ=255)
                                                                                                                                   aCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCg                   <  3:454922/149‑1 (MQ=255)
                                                                                                                                      tCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTa                >  5:338803/1‑149 (MQ=255)
                                                                                                                                        aggctgacAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACt              >  8:306378/2‑149 (MQ=255)
                                                                                                                                         ggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACtt             >  3:120222/1‑149 (MQ=255)
                                                                                                                                            tggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTc           >  2:418682/1‑148 (MQ=255)
                                                                                                                                             ggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCg                   >  6:28001/1‑139 (MQ=255)
                                                                                                                                             ggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCg                   <  5:28001/139‑1 (MQ=255)
                                                                                                                                             ggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCAACCGATGCGCCTGTTGATTCAAAGCCGCCGATAGATCGCGCTTAATTTTGAATTGCGCGACGATATCTATCGGCAAATCGGCGCTCTATTCTGCCAGCAATTCGGTACTTTc           >  6:19839/1‑147 (MQ=255)
                                                                                                                                              gcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCa         <  7:467350/148‑1 (MQ=255)
                                                                                                                                              gcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCa         <  5:127952/148‑1 (MQ=255)
                                                                                                                                                aGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGTATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTc           <  8:177207/144‑1 (MQ=255)
                                                                                                                                                aGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGGATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTc           >  7:177207/1‑144 (MQ=255)
                                                                                                                                                    tGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCg  >  5:187137/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCG  >  NC_000913/4388360‑4388656

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