Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,386,126 C→T 5.7% R1082H (CGT→CAT)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,386,1260CT 5.7% 408.1 / 10.0 123R1082H (CGT→CAT) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base C (65/51);  new base T (4/3);  total (69/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 = 7.71e-01

Only 100 of 123 total aligned reads displayed.

TGTGGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGC  >  NC_000913/4385981‑4386272
                                                                                                                                                 |                                                                                                                                                  
tgtgGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCa                                                                                                                                                 >  4:701604/1‑149 (MQ=255)
    gCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTg                                                                                                                                             >  1:458900/1‑149 (MQ=255)
       gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGaa                                                                                                                                          <  3:542718/149‑1 (MQ=255)
       gTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGaa                                                                                                                                          <  4:550548/149‑1 (MQ=255)
        ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAg                                                                                                                                         >  1:168871/1‑149 (MQ=255)
        ttttttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAg                                                                                                                                         <  4:526253/149‑1 (MQ=255)
           tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGATGCGTCTGGaaagcc                                                                                                                                      <  3:795246/149‑5 (MQ=255)
           tttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGc                                                                                                                                      <  3:759153/149‑1 (MQ=255)
            ttCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCg                                                                                                                                     <  3:291200/149‑1 (MQ=255)
             tCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCgg                                                                                                                                    <  3:701604/149‑1 (MQ=255)
                    tATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGc                                                                                                                             <  2:73100/149‑1 (MQ=255)
                    tATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGc                                                                                                                             <  1:360295/149‑1 (MQ=255)
                    tATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGc                                                                                                                             <  2:262339/149‑1 (MQ=255)
                    tATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGc                                                                                                                             <  2:836317/149‑1 (MQ=255)
                     aTCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                                                            <  1:349045/149‑1 (MQ=255)
                     aTCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCa                                                                                                                            >  3:854589/1‑149 (MQ=255)
                       cAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCata                                                                                                                          <  1:1188/149‑1 (MQ=255)
                          ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                       >  1:105534/1‑149 (MQ=255)
                          ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                       >  3:733538/1‑149 (MQ=255)
                          ttttGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCg                                                                                                                       >  4:651674/1‑149 (MQ=255)
                           tttgtttgTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGa                                                                                                                      <  2:383482/149‑1 (MQ=255)
                              gtttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAc                                                                                                                   <  2:101856/149‑1 (MQ=255)
                               tttgtGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATAcc                                                                                                                  >  1:857393/1‑149 (MQ=255)
                                   tgAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                              >  2:264048/1‑149 (MQ=255)
                                     aGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGc                                                                                                            >  4:708953/1‑149 (MQ=255)
                                       ccGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCAt                                                                                                          >  2:330662/1‑149 (MQ=255)
                                         ggATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATgg                                                                                                        >  4:21881/1‑149 (MQ=255)
                                         ggATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATgg                                                                                                        >  1:21301/1‑149 (MQ=255)
                                           aTTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                      <  2:690734/149‑1 (MQ=255)
                                               gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                              <  3:323921/137‑1 (MQ=255)
                                               gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTgg                                                                                                              >  4:323921/1‑137 (MQ=255)
                                               gTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagcca                                                                                                  >  2:357793/1‑149 (MQ=255)
                                                ttCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAagccag                                                                                                 >  2:512800/1‑149 (MQ=255)
                                                    tCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGccagccag                                                                                             >  4:215127/1‑149 (MQ=255)
                                                     ccGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGa                                                                                            >  3:464493/1‑149 (MQ=255)
                                                          acacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGCATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAg                                                                                       >  4:850192/1‑149 (MQ=255)
                                                          acacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                               >  4:240966/1‑125 (MQ=255)
                                                          acacaAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTg                                                                                                               <  3:240966/125‑1 (MQ=255)
                                                             caaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtg                                                                                    <  1:541933/149‑1 (MQ=255)
                                                             caaacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtg                                                                                    <  1:199563/149‑1 (MQ=255)
                                                               aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                  >  4:349297/1‑149 (MQ=255)
                                                               aacttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                  <  2:843453/149‑1 (MQ=255)
                                                                  ttacaaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGa                                                                               >  1:732815/1‑149 (MQ=255)
                                                                      aaacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAtctc                                                                           >  3:63338/1‑149 (MQ=255)
                                                                       aacttCCCGCCTGACGACCTTTGCCCGTAGACGTCAGCGGTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                          <  3:779149/149‑1 (MQ=255)
                                                                       aacttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCa                                                                          <  1:330662/149‑1 (MQ=255)
                                                                        acttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCAt                                                                         <  1:800513/149‑1 (MQ=255)
                                                                        acttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCAt                                                                         >  4:25384/1‑149 (MQ=255)
                                                                         cttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATg                                                                        <  2:458900/149‑1 (MQ=255)
                                                                          ttCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTCGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATgg                                                                       >  1:688300/1‑149 (MQ=255)
                                                                                      gacCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATa                                                           <  1:512800/149‑1 (MQ=255)
                                                                                        ccTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACg                                                         >  2:59197/1‑149 (MQ=255)
                                                                                        ccTTTGCCCGCAGACGTCAGCGTTCTCCCCGCTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACg                                                         >  4:142093/1‑149 (MQ=255)
                                                                                         cTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACgg                                                        <  2:21301/149‑1 (MQ=255)
                                                                                          tttGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGt                                                       >  2:729003/1‑149 (MQ=255)
                                                                                            tGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGa                                                     <  1:264048/149‑1 (MQ=255)
                                                                                               ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggtgg                                                                                >  1:292385/1‑119 (MQ=255)
                                                                                               ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggtgg                                                                                <  2:292385/119‑1 (MQ=255)
                                                                                               ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGAccc                                                  >  3:155451/1‑149 (MQ=255)
                                                                                               ccGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGAccc                                                  >  1:346796/1‑149 (MQ=255)
                                                                                                       gTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                      <  1:278582/89‑1 (MQ=255)
                                                                                                       gTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgc                                          <  2:857393/149‑1 (MQ=255)
                                                                                                       gTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGAAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGaa                                                                                                      >  2:278582/1‑89 (MQ=255)
                                                                                                          aGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcggc                                       <  3:86161/149‑1 (MQ=255)
                                                                                                          aGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAgcggc                                       >  4:86161/1‑149 (MQ=255)
                                                                                                            cGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGt                                     <  2:492732/149‑1 (MQ=255)
                                                                                                                 tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                >  2:814693/1‑149 (MQ=255)
                                                                                                                 tcCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                >  3:457072/1‑149 (MQ=255)
                                                                                                                  ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATa                                <  1:625660/148‑1 (MQ=255)
                                                                                                                  ccccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATAc                               <  2:477367/149‑1 (MQ=255)
                                                                                                                   cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                              >  2:601948/1‑149 (MQ=255)
                                                                                                                   cccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACg                              >  1:468489/1‑149 (MQ=255)
                                                                                                                    ccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGc                             <  1:923145/149‑1 (MQ=255)
                                                                                                                    ccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGc                             >  2:45372/1‑149 (MQ=255)
                                                                                                                    ccGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGc                             >  3:829953/1‑149 (MQ=255)
                                                                                                                    ccGGTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATAGCGATACCGTGGGCATGGAAGCCAGCCAGAATCACCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACgt                             <  4:779538/149‑2 (MQ=255)
                                                                                                                          gtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCTGCGTATTTACGCCGCtga                       <  3:562596/149‑3 (MQ=255)
                                                                                                                          gtttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCg                       >  3:306202/1‑149 (MQ=255)
                                                                                                                          atttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGTATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCTGCGTATTTACGCCGCtga                       <  1:669923/148‑3 (MQ=255)
                                                                                                                           tttACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGa                      >  3:678525/1‑149 (MQ=255)
                                                                                                                              aCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaaaa                   >  4:296179/1‑149 (MQ=255)
                                                                                                                              aCCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaaaa                   >  1:347724/1‑149 (MQ=255)
                                                                                                                               ccGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAAt                  <  3:651674/149‑1 (MQ=255)
                                                                                                                                cGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTTGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATc                 <  4:63338/149‑1 (MQ=255)
                                                                                                                                 gTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCt                <  4:581009/149‑1 (MQ=255)
                                                                                                                                       gCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcc          >  4:386903/1‑149 (MQ=255)
                                                                                                                                       gCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcc          >  3:37603/1‑149 (MQ=255)
                                                                                                                                        ccttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATccc         >  4:778160/3‑149 (MQ=255)
                                                                                                                                         tttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                  >  2:688552/1‑75 (MQ=255)
                                                                                                                                         tttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc        <  3:251500/149‑1 (MQ=255)
                                                                                                                                         tttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATcccc        >  2:298584/1‑149 (MQ=255)
                                                                                                                                         tttCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGCGGGCATGGAAGCCAGCCAGAATCAGCtggt                                                                                  <  1:688552/75‑1 (MQ=255)
                                                                                                                                          ttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCt       >  4:611155/1‑149 (MQ=255)
                                                                                                                                          ttCCAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCt       >  3:537654/1‑149 (MQ=255)
                                                                                                                                            ccAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAt                                                                              <  1:415413/76‑1 (MQ=255)
                                                                                                                                            ccAGATGCGTCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGAt                                                                              >  2:415413/1‑76 (MQ=255)
                                                                                                                                              aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGGAGCCAGCCAGAAGCAGCTGGTGGATCTCATGGCGTAGCGGCATACTGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGaa                     >  4:678838/1‑131 (MQ=255)
                                                                                                                                              aGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTg   >  2:912852/1‑149 (MQ=255)
                                                                                                                                               gACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGc  <  3:215127/149‑1 (MQ=255)
                                                                                                                                               gACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGc  >  3:623202/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGTGGCAGTTTTTTCTGCGCTATCAGTTTTGTTTGTGAGCCGGATTGGTTCATCCGGCACACAAACTTACAAACTTCCCGCCTGACGACCTTTGCCCGCAGACGTCAGCGTTCTCCCCGTTTGTTTACCGTTGAGGCTTTCCAGACGCATCTGGAAGGGCGGGAATGGCATATCGATACCGTGGGCATGGAAGCCAGCCAGAATCAGCTGGTGGATCTCATGGCGTAGCGGCATACGGTGACCCATCTCAGCGGCGTAAATACGCAGCTCGAAAATCTGAATCCCCTGTTGC  >  NC_000913/4385981‑4386272

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