Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,398,431 A→G 5.1% F149L (TTC→CTC)  mreD ← cell wall structural complex MreBCD transmembrane component MreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,398,4310AG 5.1% 244.9 / 2.6 78F149L (TTC→CTC) mreDcell wall structural complex MreBCD transmembrane component MreD
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base G (2/2);  total (41/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 = 9.78e-01

TCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAAT  >  NC_000913/3398284‑3398579
                                                                                                                                                   |                                                                                                                                                    
tcAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGaa                                                                                                                                                      >  1:68709/1‑148 (MQ=255)
tcAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGaa                                                                                                                                                      >  5:229169/1‑148 (MQ=255)
tcAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGaa                                                                                                                                                      >  5:210177/1‑148 (MQ=255)
      cccGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCAc                                                                                                                                               <  1:54762/149‑1 (MQ=255)
       ccGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACg                                                                                                                                              <  4:22172/149‑1 (MQ=255)
       ccGTAACAATACGTTCAAAGCTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACg                                                                                                                                              <  1:255819/149‑1 (MQ=255)
         gTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGc                                                                                                                                            <  3:299973/149‑1 (MQ=255)
         gTAAAAATACGTTCAAAGGAAACGCAAAGTTTAGCAATTAACTCATTACGACTCGGAGAACAGGAAGATAAATAAAGAGAAGTCATATAAACCTTTATTGCACTGCAAACTGCTGACGGACTATCCGCATCAGCAACAAAATCCACGGc                                                                                                                                            >  2:263989/1‑149 (MQ=255)
             cAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCaga                                                                                                                                        <  2:316204/149‑1 (MQ=255)
                    ttCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGCCGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACAGCAAACTGCTGCCGCACATTGCGCATCAGCAAGAAAATCCACgtc                                                                                                                                            >  6:362574/1‑136 (MQ=255)
                    ttCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGc                                                                                                                                            <  5:362574/138‑1 (MQ=255)
                      cAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCa                                                                                                                                <  5:251338/148‑1 (MQ=255)
                       aaaGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCAtt                                                                                                                              <  2:283435/149‑1 (MQ=255)
                       aaaGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCAtt                                                                                                                              <  6:140048/149‑1 (MQ=255)
                       aaaGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCAtt                                                                                                                              >  8:315876/1‑149 (MQ=255)
                       aaaGGTCACCCCAAGTTGCGCAAGTAACTCCTGACGACGCGCAGCACCGGAAGCTAAATACTGAGAACTGATAGAAACCTTTATTGCACTGCAAACTGCTGCCGGACTTTGCGCATCACCATGAAAATCAACGGCCACAGCACCCCAtt                                                                                                                              >  5:96074/1‑149 (MQ=255)
                        aaGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCAcctgctt                                                                                                                              >  8:180412/1‑143 (MQ=255)
                         aGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGa                                                                                                                            <  4:309522/149‑1 (MQ=255)
                               cGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGactacac                                                                                                                      >  2:256366/1‑149 (MQ=255)
                                  cAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAGAATCCACGGCCAGAGCACCCCATTGactacactac                                                                                                                   >  2:2678/1‑149 (MQ=255)
                                  cAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAGAATCCACGGCCAGAGCACCCCATTGactacactac                                                                                                                   >  2:8701/1‑149 (MQ=255)
                                     gTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTcc                                                                                                                <  2:68709/149‑1 (MQ=255)
                                         cgcAAGAAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCAAAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCCCGGCCAGAGCACCCCATTGACTACACTACTCCAGa                                                                                                             <  1:63286/148‑1 (MQ=255)
                                          gcAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGaa                                                                                                            <  2:121421/148‑1 (MQ=255)
                                          gcAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGaa                                                                                                            >  1:121421/1‑148 (MQ=255)
                                            aaGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCCCGGCCAGAGCACCCCATTGACTACACTACTCCAGAacac                                                                                                         >  4:197833/1‑149 (MQ=255)
                                             gctctcTCCTGTCGCCGCGGAGTACCGGAAGCTAAATACTGCGAAGGCATAGACTCCTTTATTGCACTGCAAACTGCTGACGGACGTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGACCACt                                                                                                        <  4:95751/144‑1 (MQ=255)
                                                 aCTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTcc                                                                                                     <  7:233067/148‑1 (MQ=255)
                                                   tCCTGACGACGCGGAGAACCGGAAGATAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAAATGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACt                                                                                                        >  3:95751/1‑143 (MQ=255)
                                                      tGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGt                                                                                                  >  1:31954/1‑146 (MQ=255)
                                                      tGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGt                                                                                                  <  2:31954/146‑1 (MQ=255)
                                                      tGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCt                                                                                                >  5:87161/1‑148 (MQ=255)
                                                         cgaagCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGaaa                                                                                            <  1:263989/149‑1 (MQ=255)
                                                              cgGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGAc                                                                                       <  6:87161/149‑1 (MQ=255)
                                                                gagaACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGAc                                                                                       <  1:290555/147‑1 (MQ=255)
                                                                gagaACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGAc                                                                                       >  2:290555/1‑147 (MQ=255)
                                                                      cGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCa                                                                               >  3:198558/1‑149 (MQ=255)
                                                                               aaTACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAAc                                                                       <  1:256366/148‑1 (MQ=255)
                                                                                 tACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAActct                                                                    <  7:315876/149‑1 (MQ=255)
                                                                                  aCAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTg                                                                   >  1:197963/1‑149 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGc                                                                  >  6:369927/1‑141 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGc                                                                  <  5:369927/141‑1 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaaaa                                                          <  6:229169/149‑1 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaaaa                                                          <  6:210177/149‑1 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaa                                                            <  1:77432/147‑1 (MQ=255)
                                                                                           tCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGaa                                                            >  2:77432/1‑147 (MQ=255)
                                                                                             aTAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAAc                                                         >  7:242632/1‑148 (MQ=255)
                                                                                                gAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTc                                                                                                      >  2:327867/1‑100 (MQ=255)
                                                                                                gAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGGTAATCACTAAAAACTCTGCCCAGAAAACaag                                                      >  6:19260/1‑147 (MQ=255)
                                                                                                     cTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATc                                                >  2:360803/1‑149 (MQ=255)
                                                                                                                tGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGc                                                                  >  4:268244/1‑120 (MQ=255)
                                                                                                                tGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGc                                                                  <  3:268244/120‑1 (MQ=255)
                                                                                                                tGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGc                                      >  6:332717/1‑148 (MQ=255)
                                                                                                                   aaaCTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGaaa                                  <  8:242632/149‑1 (MQ=255)
                                                                                                                          tgACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACAt                           >  6:102372/1‑149 (MQ=255)
                                                                                                                            aCGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATga                         >  3:59117/1‑149 (MQ=255)
                                                                                                                              ggACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgac                        >  8:995/1‑148 (MQ=255)
                                                                                                                              ggACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgac                        >  5:374195/1‑148 (MQ=255)
                                                                                                                                aCTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgacgac                     >  5:333230/1‑149 (MQ=255)
                                                                                                                                aCTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgacgac                     >  1:178492/1‑149 (MQ=255)
                                                                                                                                aCTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgacgac                     >  1:373899/1‑149 (MQ=255)
                                                                                                                                  tttGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATgacgac                     >  6:210862/1‑147 (MQ=255)
                                                                                                                                    tGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAActct                                                                    <  5:208744/98‑1 (MQ=255)
                                                                                                                                    tGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAActct                                                                    >  6:208744/1‑98 (MQ=255)
                                                                                                                                     gcgcATCAGCAACAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGATACGTTAATCACAAAAAACTCTGCCCAGAAAACAATAATATCAACCACCAGCGAAAGCAACATGACGACCaccc                >  4:216860/1‑146 (MQ=255)
                                                                                                                                         aTCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAgaga                                                                                        >  5:76138/1‑73 (MQ=255)
                                                                                                                                         aTCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAgaga                                                                                        <  6:76138/73‑1 (MQ=255)
                                                                                                                                            aGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCctgct          <  2:178492/148‑1 (MQ=255)
                                                                                                                                            aGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCctgct          <  6:333230/148‑1 (MQ=255)
                                                                                                                                             gCAAGAGAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCtgctgc        <  1:8701/149‑1 (MQ=255)
                                                                                                                                             gCAAGAGAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCtgctgc        <  1:2678/149‑1 (MQ=255)
                                                                                                                                             gCAAGAAAATCCACGGCCAGAGCATCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCtgctgc        >  5:266637/1‑149 (MQ=255)
                                                                                                                                             gCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAgaga                                                                                        <  8:29335/69‑1 (MQ=255)
                                                                                                                                             gCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCtgctgc        <  2:197963/149‑1 (MQ=255)
                                                                                                                                               aaGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATga                         <  4:273249/130‑1 (MQ=255)
                                                                                                                                               aaGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATga                         >  3:273249/1‑130 (MQ=255)
                                                                                                                                               aaGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCa      >  4:31196/1‑149 (MQ=255)
                                                                                                                                               aaGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTGCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCCACATGACGACCAGCGCCTGCTGCCa      >  1:338750/1‑149 (MQ=255)
                                                                                                                                                 gAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCGGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCAt     >  3:123087/1‑148 (MQ=255)
                                                                                                                                                 gAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGcc       <  1:126504/146‑1 (MQ=255)
                                                                                                                                                 gAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGcc       >  2:126504/1‑146 (MQ=255)
                                                                                                                                                   aaaTCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCGGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAAt  <  4:123087/149‑1 (MQ=255)
                                                                                                                                                   aaaTCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAAt  <  1:360803/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCAATGCCCGTAACAATACGTTCAAAGGTCACGCCAAGTTGCGCAAGTAACTCCTGACGACGCGGAGAACCGGAAGCTAAATACAGAGAAGTCATAGAAACCTTTATTGCACTGCAAACTGCTGACGGACTTTGCGCATCAGCAAGAAAATCCACGGCCAGAGCACCCCATTGACTACACTACTCCAGAACACTTCCGGTCTGAAAGAGACGTTAATCACTAAAAACTCTGCCCAGAAAACAATAATATCCACCACCAGCGAAAGCAACATGACGACCAGCGCCTGCTGCCATAAT  >  NC_000913/3398284‑3398579

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