Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,165 (G)7→8 coding (442/630 nt) ycfK → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,1581.G85.2% 549.2 / 90.9 189coding (435/630 nt)ycfKhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (11/17);  new base G (74/87);  total (85/104)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 230 total aligned reads displayed.

ACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTA  >  NC_000913/1207066‑1207249
                                                                                             |                                                                                            
aCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctct                                                                                               >  2:6251030/1‑93 (MQ=255)
aCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctct                                                                                               >  1:1598092/1‑93 (MQ=255)
 ccGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑g                                                                                            >  1:994586/1‑93 (MQ=255)
 ccGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑g                                                                                            >  1:2514928/1‑93 (MQ=255)
  cGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gg                                                                                           >  2:5822816/1‑93 (MQ=255)
   ggTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggg                                                                                          <  2:2236034/93‑1 (MQ=255)
   ggTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggg                                                                                          >  2:1213313/1‑93 (MQ=255)
   ggTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggg                                                                                          <  2:6605543/93‑1 (MQ=255)
    gTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggg                                                                                         <  1:7136277/93‑1 (MQ=255)
    gTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGAAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggg                                                                                         <  2:1144977/93‑1 (MQ=255)
     tttCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggggg                                                                                        >  1:2290609/1‑93 (MQ=255)
      ttCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggg                                                                                       <  1:7484501/93‑1 (MQ=255)
      ttCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggg                                                                                       <  2:4371440/93‑1 (MQ=255)
      ttCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggg                                                                                       <  1:5819962/93‑1 (MQ=255)
       tCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggggggg                                                                                      >  1:7094094/1‑93 (MQ=255)
       tCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggggggg                                                                                      <  1:6262439/93‑1 (MQ=255)
        cATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggggg                                                                                     <  1:1500807/93‑2 (MQ=255)
         atatCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggggga                                                                                    <  2:4893794/93‑3 (MQ=255)
         aTAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑gggggggga                                                                                    >  2:7272906/1‑91 (MQ=255)
         aTAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGAt                                                                                    >  2:4920052/1‑93 (MQ=255)
          tAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑ggggggggat                                                                                   <  2:4094200/93‑4 (MQ=255)
             cTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGa                                                                                 <  2:6708996/93‑1 (MQ=255)
               ataCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGaaa                                                                               <  1:2894640/93‑1 (MQ=255)
               ataCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGaaa                                                                               <  1:7311535/93‑1 (MQ=255)
               ataCCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAAt                                                                              >  2:541470/1‑93 (MQ=255)
                   cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACa                                                                           <  2:1543138/93‑1 (MQ=255)
                   cGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACa                                                                           <  2:5632268/93‑1 (MQ=255)
                    ggTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAg                                                                          <  1:2582089/93‑1 (MQ=255)
                     gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAg                                                                          >  1:3954175/1‑93 (MQ=255)
                     gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGt                                                                         <  1:1064524/93‑1 (MQ=255)
                     gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGt                                                                         <  1:3238879/93‑1 (MQ=255)
                     gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGtt                                                                        <  2:5161510/93‑1 (MQ=255)
                      ttATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGtt                                                                        >  1:7474305/1‑93 (MQ=255)
                        aTGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTg                                                                       <  2:2158441/93‑1 (MQ=255)
                        aTGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGt                                                                      <  1:3773195/93‑1 (MQ=255)
                         tGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTa                                                                     >  2:3141030/1‑93 (MQ=255)
                         tGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTa                                                                     >  2:5510441/1‑93 (MQ=255)
                         tGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTa                                                                     <  1:5268762/93‑1 (MQ=255)
                          gATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTAt                                                                    >  1:5537212/1‑93 (MQ=255)
                           aTTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTAtt                                                                   >  1:1525291/1‑93 (MQ=255)
                            ttGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTc                                                                  >  1:1854898/1‑93 (MQ=255)
                              ggtggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTccc                                                                <  1:1038147/93‑1 (MQ=255)
                              ggtggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTccc                                                                >  2:1515505/1‑93 (MQ=255)
                                 ggtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGt                                                             >  2:3535314/1‑93 (MQ=255)
                                 ggtATTGAGCAATCAGTAATCATTCAGGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCgg                                                              <  1:361599/93‑1 (MQ=255)
                                  gtATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGtt                                                            <  1:597962/93‑1 (MQ=255)
                                   tATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTg                                                           <  2:5080570/93‑1 (MQ=255)
                                       gAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGtt                                                       <  1:3538680/93‑1 (MQ=255)
                                         gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTta                                                     <  1:1390410/93‑1 (MQ=255)
                                         gCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaa                                                    <  1:4991903/93‑1 (MQ=255)
                                          cAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaa                                                    <  1:5089297/93‑1 (MQ=255)
                                           aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                   <  2:5598265/93‑1 (MQ=255)
                                           aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                   >  1:5166440/1‑93 (MQ=255)
                                           aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                   >  2:5765733/1‑93 (MQ=255)
                                           aaTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                   >  2:4555305/1‑93 (MQ=255)
                                            aTCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaata                                                  >  1:1690725/1‑93 (MQ=255)
                                             ttaGTAATCATGCAGTGGGGGTGGAACGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaataa                                                 <  1:1768470/91‑1 (MQ=255)
                                             tCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaataa                                                 <  2:870932/93‑1 (MQ=255)
                                              cAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTaacaat                                                <  1:6974462/93‑1 (MQ=255)
                                               aGTAATCATTCAGTGGGGGGGGGATGCCGCAAAAACAACTGCCTCTGGGGGGGGGTGGAAATACAGATGTATTCCCGGTTGCGTTTTATAATg                                               >  2:285145/1‑93 (MQ=255)
                                                gTAATCATGCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGc                                              <  1:5509260/93‑1 (MQ=255)
                                                 tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGcc                                             >  2:3730719/1‑93 (MQ=255)
                                                 tAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGcc                                             >  2:4754224/1‑93 (MQ=255)
                                                  aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCt                                            <  2:5299659/93‑1 (MQ=255)
                                                  aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCt                                            <  1:4084774/93‑1 (MQ=255)
                                                  aaTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCt                                            >  2:7020091/1‑93 (MQ=255)
                                                   aTCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtg                                           <  2:1930450/93‑1 (MQ=255)
                                                     cATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtgtg                                         <  2:7485656/93‑1 (MQ=255)
                                                     cATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTgtgt                                        <  1:4636433/93‑1 (MQ=255)
                                                      aTTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt                                       <  2:1482401/93‑1 (MQ=255)
                                                       ttCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGtt                                       <  2:457839/93‑1 (MQ=255)
                                                         cAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGc                                     <  1:4519042/93‑1 (MQ=255)
                                                              ggggTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgttgt                               <  1:6082607/93‑1 (MQ=255)
                                                               gggTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCgttgt                               <  1:2378631/93‑1 (MQ=255)
                                                                 gTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                             >  2:385400/1‑93 (MQ=255)
                                                                    gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaa                          >  1:1264047/1‑93 (MQ=255)
                                                                    gAATGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaa                          <  1:5376845/93‑1 (MQ=255)
                                                                     aaTGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaaa                         >  2:6033829/1‑93 (MQ=255)
                                                                     aaTGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaaa                         >  2:2632504/1‑93 (MQ=255)
                                                                      aTGCCGCAAAAGCATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaaa                         >  2:7078534/1‑93 (MQ=255)
                                                                       tGCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAAtt                       <  1:875420/93‑1 (MQ=255)
                                                                       tGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTa                      <  2:5940021/93‑1 (MQ=255)
                                                                       tGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTa                      <  2:296170/93‑1 (MQ=255)
                                                                        gCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTa                      <  2:2782005/93‑1 (MQ=255)
                                                                        gCCGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTa                      <  2:5996048/93‑1 (MQ=255)
                                                                         ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTAt                     <  2:6869682/93‑1 (MQ=255)
                                                                         ccGCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTAt                     <  1:896971/93‑1 (MQ=255)
                                                                           gCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGa                   >  2:6673536/1‑93 (MQ=255)
                                                                           gCAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGa                   >  1:6851384/1‑93 (MQ=255)
                                                                            cAAAAGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGAc                  <  1:6411747/93‑1 (MQ=255)
                                                                             aaaaGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                <  1:3058427/93‑1 (MQ=255)
                                                                              aaaGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                <  2:4779995/93‑1 (MQ=255)
                                                                              aaaGCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                <  1:6114659/93‑1 (MQ=255)
                                                                               aaGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATg              >  2:2281317/1‑93 (MQ=255)
                                                                                 gCATCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGt             >  2:3947473/1‑93 (MQ=255)
                                                                                  cATCTGCCTCTGGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGt             <  1:3349265/93‑1 (MQ=255)
                                                                                   aTCTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa           >  1:3503018/1‑93 (MQ=255)
                                                                                     cTGCCTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgc         <  1:6536511/93‑1 (MQ=255)
                                                                                        ccTCTGGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa      <  2:4120309/93‑1 (MQ=255)
                                                                                           ct‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTa  >  2:7531385/1‑93 (MQ=255)
                                                                                             |                                                                                            
ACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCT‑‑GGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTA  >  NC_000913/1207066‑1207249

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