Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,624 Δ1 bp 15.3% coding (426/822 nt) yfcO ← DUF2544 family putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,449,6240C.15.3% 468.7 / 85.3 157coding (426/822 nt)yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (58/74);  new base . (12/12);  total (70/87)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 158 total aligned reads displayed.

GTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449527‑2449723
                                                                                                  |                                                                                                    
  tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                                  <  2:1081170/100‑1 (MQ=255)
  tcttccgatcTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                                  <  2:1963881/91‑1 (MQ=255)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATGGTGTAGTGAC‑TATCAGACCCCTGATACGTCACg                                                                                                >  2:354074/1‑100 (MQ=255)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCAc                                                                                                 <  1:298679/99‑1 (MQ=38)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                                                <  2:1263405/100‑1 (MQ=255)
     gccTTACTGCTATAATTCCGAGGCCAAGAAACGTCATAAGGGAGCACAGAACCGGATGGACAAATCGCGACGTGAA‑CAGCAGGCCCCTGAGACGTCa                                                                                                  >  2:3565417/3‑97 (MQ=255)
       cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑ccc                                                                                            >  1:2720088/1‑100 (MQ=255)
       cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑ccc                                                                                            >  2:2599508/1‑100 (MQ=255)
        ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                            <  2:114219/100‑5 (MQ=37)
         tACTGACCTAATTCTCATTCCATGAAACGTCATAATGTTGAACAGACACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacacg                                                                                           <  1:1499321/100‑6 (MQ=21)
         tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑ccccc                                                                                          >  2:2567986/1‑100 (MQ=255)
         tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑ccccc                                                                                          <  1:3408156/100‑1 (MQ=255)
           cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGCAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTg                                                                                        >  2:150323/1‑100 (MQ=255)
            tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGGCGTGACGTATCAGGGGCCTGATA‑GTCACgacacgatt                                                                                        <  1:28926/100‑9 (MQ=21)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGga                                                                                      >  2:349383/1‑100 (MQ=255)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGga                                                                                      >  1:530225/1‑100 (MQ=255)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGga                                                                                      >  1:271029/1‑100 (MQ=255)
              aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑TCCCCTGgag                                                                                     >  2:2766639/1‑100 (MQ=255)
               cATAATTCTCATTCCATGACACGTCATAATGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGgag                                                                                     <  2:720289/99‑1 (MQ=255)
                aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGagag                                                                                   <  1:411752/100‑1 (MQ=255)
                 tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAgg                                                                                  >  1:2304457/1‑100 (MQ=255)
                  aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAgg                                                                                  >  1:2868589/1‑99 (MQ=255)
                   aTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TAGCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCt                                                                                <  2:527089/100‑1 (MQ=255)
                    ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCtt                                                                               <  2:1204232/100‑1 (MQ=38)
                     tctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTc                                                                              <  2:1664180/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCa                                                                             <  1:789885/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCa                                                                             <  1:150054/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCa                                                                             >  1:1534484/1‑100 (MQ=255)
                       tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCag                                                                            <  2:2954087/100‑1 (MQ=255)
                       tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCag                                                                            <  1:2630833/100‑1 (MQ=255)
                        cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCGCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCcga                                                                           >  2:1010930/1‑100 (MQ=255)
                        cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCaga                                                                           >  2:2012709/1‑100 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGTATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGGTTCagag                                                                          <  2:1907885/100‑1 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagag                                                                          <  2:2827653/100‑1 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagag                                                                          >  1:625474/1‑100 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagag                                                                          >  2:827780/1‑100 (MQ=255)
                           tCCATTAGACGTCATGAGGTGGAGCAGATACTTATGGACAAATCGTGTCGTGAC‑TAGCGGGGCCCTGATACGTGACG‑CCCCCTGGAGAGGCTTCagag                                                                          >  2:722728/1‑98 (MQ=255)
                            ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtca                                                                       <  2:848678/100‑1 (MQ=255)
                             cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCagagggcg                                                                      >  1:2232281/1‑96 (MQ=255)
                             cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcat                                                                      >  1:2169487/1‑100 (MQ=255)
                              aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                                     <  2:3201210/100‑1 (MQ=255)
                              aTGAAACGTCATAAGGTCGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATGAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcatt                                                                     >  2:2386500/1‑100 (MQ=255)
                               tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattc                                                                    <  1:1431776/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  1:3329049/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  2:1931304/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  2:530223/100‑1 (MQ=255)
                                   aCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCatan                                                                <  2:1306214/100‑2 (MQ=255)
                                    cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                               <  1:3508339/100‑1 (MQ=255)
                                    cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                               >  1:2843062/1‑100 (MQ=255)
                                       tatatGGTTGTACACAAACGGATGGACAAATCGTGCTGTGAC‑TATTGTGCCGCTTGTATGTCACG‑CCCCCTGCAGAGGCTTCAGGGGCATTTATAATaac                                                            <  2:1468563/99‑1 (MQ=255)
                                       cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  1:3046308/100‑1 (MQ=255)
                                        aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                           <  1:3029512/100‑1 (MQ=255)
                                        aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCACCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  2:3262659/99‑1 (MQ=255)
                                            ggTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                       >  1:189950/1‑100 (MQ=255)
                                             gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAA‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGGGGCTTCAGAGTCATTCATAATAaccaacca                                                      >  1:3208648/1‑100 (MQ=255)
                                                gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   >  1:2456017/1‑100 (MQ=255)
                                                 aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                  >  1:1151813/1‑100 (MQ=255)
                                                  aCAGAAACGGATGGACAAATCGTCTCGTGAC‑TATCAGGCCCCTTATACTTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                 <  2:3469036/100‑1 (MQ=255)
                                                     gAAACGGATGGACAAATCGTGTCGTGCC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGaaaa                                               >  2:865302/1‑99 (MQ=255)
                                                     gAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                              >  2:1708398/1‑100 (MQ=255)
                                                          ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:1462876/100‑1 (MQ=255)
                                                          cttccgatctaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:162496/90‑1 (MQ=255)
                                                               gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     >  2:3388504/1‑100 (MQ=255)
                                                               gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                    >  1:1714704/1‑100 (MQ=255)
                                                                aCAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                                   <  1:3575/100‑1 (MQ=255)
                                                                     tcgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                              <  1:836050/100‑1 (MQ=255)
                                                                      cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             >  1:437328/1‑100 (MQ=255)
                                                                       gtgttgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  1:1807694/99‑1 (MQ=255)
                                                                         gtcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          <  1:3536882/100‑1 (MQ=255)
                                                                         ggcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:1029726/3‑100 (MQ=37)
                                                                         ggcgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:2574354/3‑100 (MQ=37)
                                                                          tcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         <  1:1781259/100‑1 (MQ=255)
                                                                          tcgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         <  2:444730/100‑1 (MQ=255)
                                                                           cgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTgt                        >  1:3070965/1‑99 (MQ=37)
                                                                           cgtgACGTATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:1741671/1‑100 (MQ=38)
                                                                           cgtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:1642467/1‑100 (MQ=255)
                                                                            gtgacgtATCAGGGGCCTGATA‑GTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       >  1:339017/7‑100 (MQ=38)
                                                                            gtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:108665/100‑1 (MQ=255)
                                                                            gtgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:1960849/100‑1 (MQ=255)
                                                                             tgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      >  1:3255970/1‑100 (MQ=255)
                                                                             tgAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  1:419310/100‑1 (MQ=255)
                                                                              gAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                     >  1:822828/1‑100 (MQ=255)
                                                                                  tATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                  <  1:2346382/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  1:3621551/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  2:3310599/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  1:3621750/100‑1 (MQ=255)
                                                                                     cAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt               >  1:316198/1‑100 (MQ=255)
                                                                                      aGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa              <  2:2201179/100‑1 (MQ=255)
                                                                                       ggCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa             >  2:1530673/1‑100 (MQ=255)
                                                                                       ggCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa             >  1:1760075/1‑100 (MQ=255)
                                                                                         ccccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           >  1:294567/1‑100 (MQ=255)
                                                                                         ccccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           <  2:2304441/100‑1 (MQ=255)
                                                                                         ccccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGACTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATATCTACACAGACAAGAAc            >  2:3184315/1‑99 (MQ=255)
                                                                                          cccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt          <  1:2383794/100‑1 (MQ=255)
                                                                                          cccTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt          <  1:780121/100‑1 (MQ=255)
                                                                                            cTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa        >  1:2516601/1‑100 (MQ=255)
                                                                                              gATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  1:2055827/1‑100 (MQ=255)
                                                                                              gATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  1:2006102/1‑100 (MQ=255)
                                                                                              gATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  2:1196049/1‑100 (MQ=255)
                                                                                               aTACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAAt     <  1:909672/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACG‑CCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449527‑2449723

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