Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,326,411 G→C 6.2% S175R (AGC→AGG)  ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,326,4110GC 6.2% 219.2 / 6.7 65S175R (AGC→AGG) ftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (34/27);  new base C (2/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

ATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGC  >  NC_000913/3326265‑3326553
                                                                                                                                                  |                                                                                                                                              
aTTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGcggct                                                                                                                                              >  2:541543/1‑149 (MQ=255)
    cTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGcggctcggc                                                                                                                                          <  2:310969/149‑1 (MQ=255)
      aCGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCggctcggctc                                                                                                                                        <  2:346159/149‑1 (MQ=255)
           gTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCa                                                                                                                                   >  4:458243/1‑149 (MQ=255)
                aaCCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATAc                                                                                                                              <  2:464353/149‑1 (MQ=255)
                 aCCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACt                                                                                                                             >  4:58274/1‑149 (MQ=255)
                        aTAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAg                                                                                                                      <  3:51904/149‑1 (MQ=255)
                             aaaGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAg                                                                                                                 >  2:246384/1‑149 (MQ=255)
                              aaGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGc                                                                                                                >  3:108591/1‑149 (MQ=255)
                                 aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAc                                                                                                             >  3:35110/1‑149 (MQ=255)
                                 aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAc                                                                                                             >  1:220231/1‑149 (MQ=255)
                                 aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAc                                                                                                             <  4:355557/149‑1 (MQ=255)
                                 aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAc                                                                                                             >  1:443047/1‑149 (MQ=255)
                                 aaCGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAAc                                                                                                             <  4:345926/149‑1 (MQ=255)
                                    ggAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGCCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtc                                                                                                          >  1:354125/1‑149 (MQ=255)
                                    ggAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGCCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtc                                                                                                          >  1:145404/1‑149 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGGCGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtcttctt                                                                                                     >  4:516383/1‑149 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAAc                                                                                                                         >  4:386579/1‑129 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAAc                                                                                                                         <  3:386579/129‑1 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtcttctt                                                                                                     >  2:357763/1‑149 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtcttctt                                                                                                     >  1:545828/1‑149 (MQ=255)
                                         tttCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAActtcttctt                                                                                                     >  1:183482/1‑149 (MQ=255)
                                             gcttcgcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGc                                                                                                 <  1:209897/149‑1 (MQ=255)
                                                tcgcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTtc                                                                                              <  3:512868/149‑1 (MQ=255)
                                                  gcCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTtcgt                                                                                            <  3:90260/149‑1 (MQ=255)
                                                   ccTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTtcgtc                                                                                           <  4:94626/149‑1 (MQ=255)
                                                       cAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAg                                                                                       <  4:290418/149‑1 (MQ=255)
                                                         aTCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGcc                                                                                     >  1:100304/1‑149 (MQ=255)
                                                          tCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGccc                                                                                    >  2:95309/1‑149 (MQ=255)
                                                              tttCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCaa                                                                                >  3:458638/1‑149 (MQ=255)
                                                                   ccAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCa                                                                           <  4:202747/149‑1 (MQ=255)
                                                                        agcGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCaaa                                                                      <  4:253219/149‑1 (MQ=255)
                                                                         gcGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAg                                                                     >  2:56625/1‑149 (MQ=255)
                                                                           gTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGt                                                                   >  4:123005/1‑149 (MQ=255)
                                                                            ttttACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTc                                                                  >  1:568222/1‑149 (MQ=255)
                                                                               tACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGtt                                                               <  3:516429/149‑1 (MQ=255)
                                                                                aCCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAgg                                                                    >  1:173380/1‑143 (MQ=255)
                                                                                aCCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGttt                                                              >  3:396171/1‑149 (MQ=255)
                                                                                aCCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCACTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAgg                                                                    <  2:173380/143‑1 (MQ=255)
                                                                                       cccGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCACGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTtgatctg                                                       >  1:24032/1‑149 (MQ=255)
                                                                                          ggaggaCCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTtgatctgatc                                                    >  2:29788/1‑149 (MQ=255)
                                                                                            aggaCCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCtt                                                  <  4:23897/149‑1 (MQ=255)
                                                                                                  gaccATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTtgatctgat                                                     >  1:584828/1‑140 (MQ=255)
                                                                                                  gaccATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTtgatctgat                                                     <  2:584828/140‑1 (MQ=255)
                                                                                                   accATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAg                                           >  1:205306/1‑149 (MQ=255)
                                                                                                    ccATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGc                                          >  2:415586/1‑149 (MQ=255)
                                                                                                       tCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATg                                       >  4:72547/1‑149 (MQ=255)
                                                                                                        cAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATgc                                      <  4:458638/149‑1 (MQ=255)
                                                                                                        cAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTAAGCATgc                                      <  2:312234/149‑1 (MQ=255)
                                                                                                          aGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGTGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATgcgc                                    >  2:58432/1‑149 (MQ=255)
                                                                                                          aGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCAGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATgcgc                                    >  2:538283/1‑149 (MQ=255)
                                                                                                           gACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATgcgcg                                   <  1:502851/149‑1 (MQ=255)
                                                                                                             cGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCt                                 >  2:377592/1‑149 (MQ=255)
                                                                                                              gCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCtt                                <  4:38057/149‑1 (MQ=255)
                                                                                                               ccTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCCTCCGTCAGCATGCGCGCttt                               <  1:226318/149‑1 (MQ=255)
                                                                                                                 tttCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGc                             >  4:337505/1‑149 (MQ=255)
                                                                                                                     gggATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCtt                         <  4:335999/149‑1 (MQ=255)
                                                                                                                     gggATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCtt                         <  3:101600/149‑1 (MQ=255)
                                                                                                                      ggATCTTACCGCCGAGTTTCTGGAAGCGCCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTa                        <  2:354125/149‑1 (MQ=255)
                                                                                                                      ggATCTTACCGCCGAGTTTCTGGAAGCGCCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTa                        <  2:145404/149‑1 (MQ=255)
                                                                                                                         tCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCa                     >  1:585536/1‑149 (MQ=255)
                                                                                                                           ttACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCaaa                   <  1:342577/149‑1 (MQ=255)
                                                                                                                            tACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAAc                  >  4:553338/1‑149 (MQ=255)
                                                                                                                                gccgAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACa              <  3:123005/149‑1 (MQ=255)
                                                                                                                                           ggAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGCCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTg   >  2:473551/1‑149 (MQ=255)
                                                                                                                                            gAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGc  <  2:545828/149‑1 (MQ=255)
                                                                                                                                            gAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGc  <  2:585536/149‑1 (MQ=255)
                                                                                                                                            gAAGCGGCTAGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGc  >  2:134426/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ATTTCTACGAAGTCAGAACCGGAGATAGTAAAGAACGGAACTTTCGCTTCGCCTGCAATCGCTTTCGCCAGCAGCGTTTTACCGGTACCCGGAGGACCGACCATCAAGACGCCTTTCGGGATCTTACCGCCGAGTTTCTGGAAGCGGCTCGGCTCGCGCAGATACTCAACCAGTTCAGCAACTTCTTCTTTTGCTTCGTCGCAGCCCGCAACGTCAGCAAAGGTCGTTTTGATCTGATCTTCCGTCAGCATGCGCGCTTTGCTCTTACCAAACGACATGGCACCTTTGC  >  NC_000913/3326265‑3326553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑