Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,444,419 G→T 13.2% L157L (CTC→CTA) ‡ sulA ← cell division inhibitor SulA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,444,4190GT13.2% 192.0 / 16.1 68L157L (CTC→CTA) ‡sulAcell division inhibitor SulA
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (3/6);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCG  >  NC_002947/2444277‑2444563
                                                                                                                                               |                                                                                                                                                  
ccATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCgg                                                                                                                                              <  1:247714/149‑1 (MQ=255)
      tCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGtt                                                                                                                                         >  5:42850/1‑148 (MQ=255)
        aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCa                                                                                                                                       <  5:122476/148‑1 (MQ=255)
        aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCAg                                                                                                                                      <  8:157731/149‑1 (MQ=255)
        aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTCAg                                                                                                                                      <  5:202484/149‑1 (MQ=255)
        aGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTc                                                                                                                                        <  7:29790/147‑1 (MQ=255)
          tgtgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGc                                                                                                                                    >  3:301603/1‑149 (MQ=255)
            tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTTAGTGTGCAGCAGCGTATCAGGCTATTTGCAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                   <  1:71884/149‑1 (MQ=255)
            tgGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCt                                                                                                                                   <  4:163969/148‑1 (MQ=255)
             gggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTc                                                                                                                                        >  5:213112/1‑142 (MQ=255)
             gggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCCGCCTCGGCTGATGTTc                                                                                                                                        <  6:213112/142‑1 (MQ=255)
              ggCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCttt                                                                                                                                 <  8:16692/148‑1 (MQ=255)
                cGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTtgtg                                                                                                                              <  7:63886/149‑1 (MQ=255)
                        tGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTcc                                                                                                                       <  4:137271/148‑1 (MQ=255)
                           ccGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGc                                                                                                                    >  1:45929/1‑148 (MQ=255)
                            cGGTCTCCTGATGGATGGCGATCTCGAGCGTGATGAAGGGGCGCGTTTGGTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTTAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCa                                                                                                                   <  2:45929/148‑1 (MQ=255)
                               tCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTgg                                                                                                               <  6:122567/149‑1 (MQ=255)
                                   cTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGc                                                                                                           <  3:239710/149‑1 (MQ=255)
                                    tGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCa                                                                                                          >  8:158151/1‑149 (MQ=255)
                                              cGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGAACGTAGCAgtt                                                                                                >  4:251437/1‑149 (MQ=255)
                                                aGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAgttgt                                                                                              <  5:233040/149‑1 (MQ=255)
                                                 gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                       >  6:109995/1‑107 (MQ=255)
                                                 gCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCa                                                                                                                                       <  5:109995/107‑1 (MQ=255)
                                                  cTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAgttgtt                                                                                             <  6:5147/148‑1 (MQ=255)
                                                           gatgaAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGc                                                                                   >  6:144134/1‑149 (MQ=255)
                                                                   gCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGa                                                                         >  5:176468/1‑149 (MQ=255)
                                                                    ccGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGAc                                                                        >  7:293221/1‑149 (MQ=255)
                                                                           tgctgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAACGTCA‑‑ATGCCGCTGcac                                                                 <  7:85570/149‑3 (MQ=255)
                                                                              tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagc                                                              >  2:103448/1‑149 (MQ=255)
                                                                              tgcACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagc                                                              >  3:208378/1‑149 (MQ=255)
                                                                                cACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCccagcc                                                             <  6:42850/148‑1 (MQ=255)
                                                                                   tGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCt                                                         >  7:81224/1‑149 (MQ=255)
                                                                                     aGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTca                                                       >  4:249/1‑149 (MQ=255)
                                                                                            agcagcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCATGgt                                                <  4:301603/149‑1 (MQ=255)
                                                                                              cagcGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACTTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGgtgt                                              >  5:269803/1‑149 (MQ=255)
                                                                                                tcGTATCAGGCTATTTGAAGTTACAGGCAAAGCATGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCAAGTAGCAGTTGTTGGAGTGCAGCGGAA‑‑TTGACGTTGCCCAGCCAGCTCAACACGGtgtgt                                            <  6:269803/148‑2 (MQ=255)
                                                                                                   tATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTg                                         <  8:293221/149‑1 (MQ=255)
                                                                                                      cAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                      >  6:308069/1‑149 (MQ=255)
                                                                                                      cAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAgg                                      >  2:256680/1‑149 (MQ=255)
                                                                                                       aGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑tt                                                                           >  8:258196/1‑111 (MQ=255)
                                                                                                       aGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑tt                                                                           <  7:258196/111‑1 (MQ=255)
                                                                                                        ggCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                    >  4:147943/1‑149 (MQ=255)
                                                                                                        ggCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                    >  2:234462/1‑149 (MQ=255)
                                                                                                         gCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCcca                                                                >  3:123760/1‑120 (MQ=255)
                                                                                                         gCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCcca                                                                <  4:123760/120‑1 (MQ=255)
                                                                                                         gCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGcc                                    >  2:72602/1‑148 (MQ=255)
                                                                                                                    gTGACAGGCAAAGGCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          <  7:158151/148‑1 (MQ=255)
                                                                                                                    gTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                         >  5:5664/1‑149 (MQ=255)
                                                                                                                     tGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGc                        >  4:88147/1‑149 (MQ=255)
                                                                                                                          ggCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCAGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGACCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTca                    >  5:255675/1‑148 (MQ=255)
                                                                                                                           gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCt                                          <  4:124811/125‑1 (MQ=255)
                                                                                                                           gCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCt                                          >  3:124811/1‑125 (MQ=255)
                                                                                                                            cAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                 <  5:308069/149‑1 (MQ=255)
                                                                                                                             aaaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAgg                >  7:11466/1‑149 (MQ=255)
                                                                                                                              aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCATCTTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAg                 <  1:256680/149‑1 (MQ=255)
                                                                                                                              aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          >  4:178191/1‑138 (MQ=255)
                                                                                                                              aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          <  3:178191/138‑1 (MQ=255)
                                                                                                                              aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGc               <  8:81224/149‑1 (MQ=255)
                                                                                                                              aaGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGc               >  3:286845/1‑149 (MQ=255)
                                                                                                                                  cTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGCGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                         >  6:146660/1‑135 (MQ=39)
                                                                                                                                  cTGAACATCCGCCTCGGCTGATGTTCAGGCTTTGCGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAg                         <  5:146660/135‑1 (MQ=39)
                                                                                                                                   tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          <  7:283238/133‑1 (MQ=255)
                                                                                                                                   tGAACATCCGCCTCGGCTGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCaaa                          >  8:283238/1‑133 (MQ=255)
                                                                                                                                      aCATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCa       <  6:176468/149‑1 (MQ=255)
                                                                                                                                        aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                               <  2:121824/123‑1 (MQ=255)
                                                                                                                                        aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGt                               >  1:121824/1‑123 (MQ=255)
                                                                                                                                        aTCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgc     <  3:251437/149‑1 (MQ=255)
                                                                                                                                           aGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAgcgcg  >  8:30142/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CCATGCTCAGTGTGGGCGTGCAGATGGCCGGTCGCCTGATGGATGGCGAGCTCGAGCGTGATGAAGGGCCGCGTGTGCTGCACTGAGTGTGCAGCAGCGTATCAGGCTATTTGAAGTGACAGGCAAAGCCTGAACATCAGCCGAGGCGGATGTTCAGGCTTTGTGCATTTCCGGCACTGGCTGCACGTAGCAGTTGTTGGCGTGCAGCGGCA‑‑TTGACGTTGCCCAGCCAGCTCACCACGGTGTGGCTGAGGCCCAGGTGCAAAGCCTCACAGGCCAGTTGCAGCGCG  >  NC_002947/2444277‑2444563

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: