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.

GTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449527‑2449723
                                                                                                  |                                                                                                    
  tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                                  <  2:1081170/100‑1 (MQ=255)
  tcttccgatcTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                                  <  2:1963881/91‑1 (MQ=255)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATGGTGTAGTGACTATCAGACCCCTGATACGTCACg                                                                                                >  2:354074/1‑100 (MQ=255)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCAc                                                                                                 <  1:298679/99‑1 (MQ=38)
    tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACg                                                                                                <  2:1263405/100‑1 (MQ=255)
     gccTTACTGCTATAATTCCGAGGCCAAGAAACGTCATAAGGGAGCACAGAACCGGATGGACAAATCGCGACGTGAACAGCAGGCCCCTGAGACGTCa                                                                                                  >  2:3565417/3‑97 (MQ=255)
       cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                            >  1:2720088/1‑100 (MQ=255)
       cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccc                                                                                            >  2:2599508/1‑100 (MQ=255)
        ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                            <  2:114219/100‑5 (MQ=37)
         tACTGACCTAATTCTCATTCCATGAAACGTCATAATGTTGAACAGACACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacacg                                                                                           <  1:1499321/100‑6 (MQ=21)
         tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                          >  2:2567986/1‑100 (MQ=255)
         tACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGccccc                                                                                          <  1:3408156/100‑1 (MQ=255)
           cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGCAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTg                                                                                        >  2:150323/1‑100 (MQ=255)
            tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGGCGTGACGTATCAGGGGCCTGATAGTCACgacacgatt                                                                                        <  1:28926/100‑9 (MQ=21)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  2:349383/1‑100 (MQ=255)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:530225/1‑100 (MQ=255)
             gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:271029/1‑100 (MQ=255)
              aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGTCCCCTGgag                                                                                     >  2:2766639/1‑100 (MQ=255)
               cATAATTCTCATTCCATGACACGTCATAATGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                     <  2:720289/99‑1 (MQ=255)
                aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                   <  1:411752/100‑1 (MQ=255)
                 tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                  >  1:2304457/1‑100 (MQ=255)
                  aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                  >  1:2868589/1‑99 (MQ=255)
                   aTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTAGCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCt                                                                                <  2:527089/100‑1 (MQ=255)
                    ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCtt                                                                               <  2:1204232/100‑1 (MQ=38)
                     tctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTc                                                                              <  2:1664180/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                             <  1:789885/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                             <  1:150054/100‑1 (MQ=255)
                      ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                             >  1:1534484/1‑100 (MQ=255)
                       tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                            <  2:2954087/100‑1 (MQ=255)
                       tcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCag                                                                            <  1:2630833/100‑1 (MQ=255)
                        cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCcga                                                                           >  2:1010930/1‑100 (MQ=255)
                        cATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCaga                                                                           >  2:2012709/1‑100 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGTATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGGTTCagag                                                                          <  2:1907885/100‑1 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                          <  2:2827653/100‑1 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                          >  1:625474/1‑100 (MQ=255)
                         aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagag                                                                          >  2:827780/1‑100 (MQ=255)
                           tCCATTAGACGTCATGAGGTGGAGCAGATACTTATGGACAAATCGTGTCGTGACTAGCGGGGCCCTGATACGTGACGCCCCCTGGAGAGGCTTCagag                                                                          >  2:722728/1‑98 (MQ=255)
                            ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                                       <  2:848678/100‑1 (MQ=255)
                             cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCagagggcg                                                                      >  1:2232281/1‑96 (MQ=255)
                             cATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcat                                                                      >  1:2169487/1‑100 (MQ=255)
                              aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                     <  2:3201210/100‑1 (MQ=255)
                              aTGAAACGTCATAAGGTCGAACAGAAACGGATGGACAAATCGTGTCGTGACTATGAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                     >  2:2386500/1‑100 (MQ=255)
                               tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                    <  1:1431776/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  1:3329049/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  2:1931304/100‑1 (MQ=255)
                                 aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  2:530223/100‑1 (MQ=255)
                                   aCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCatan                                                                <  2:1306214/100‑2 (MQ=255)
                                    cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                               <  1:3508339/100‑1 (MQ=255)
                                    cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataat                                                               >  1:2843062/1‑100 (MQ=255)
                                       tatatGGTTGTACACAAACGGATGGACAAATCGTGCTGTGACTATTGTGCCGCTTGTATGTCACGCCCCCTGCAGAGGCTTCAGGGGCATTTATAATaac                                                            <  2:1468563/99‑1 (MQ=255)
                                       cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  1:3046308/100‑1 (MQ=255)
                                        aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                                           <  1:3029512/100‑1 (MQ=255)
                                        aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCACCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  2:3262659/99‑1 (MQ=255)
                                            ggTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaaccaacc                                                       >  1:189950/1‑100 (MQ=255)
                                             gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAATATCAGGCCCCTGATACGTCACGCCCCCTGGAGGGGCTTCAGAGTCATTCATAATAaccaacca                                                      >  1:3208648/1‑100 (MQ=255)
                                                gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   >  1:2456017/1‑100 (MQ=255)
                                                 aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                  >  1:1151813/1‑100 (MQ=255)
                                                  aCAGAAACGGATGGACAAATCGTCTCGTGACTATCAGGCCCCTTATACTTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGAc                                                 <  2:3469036/100‑1 (MQ=255)
                                                     gAAACGGATGGACAAATCGTGTCGTGCCTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGaaaa                                               >  2:865302/1‑99 (MQ=255)
                                                     gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                              >  2:1708398/1‑100 (MQ=255)
                                                          ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:1462876/100‑1 (MQ=255)
                                                          cttccgatctaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:162496/90‑1 (MQ=255)
                                                               gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     >  2:3388504/1‑100 (MQ=255)
                                                               gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                    >  1:1714704/1‑100 (MQ=255)
                                                                aCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                                   <  1:3575/100‑1 (MQ=255)
                                                                     tcgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAt                              <  1:836050/100‑1 (MQ=255)
                                                                      cgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             >  1:437328/1‑100 (MQ=255)
                                                                       gtgttgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  1:1807694/99‑1 (MQ=255)
                                                                         gtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          <  1:3536882/100‑1 (MQ=255)
                                                                         ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:1029726/3‑100 (MQ=37)
                                                                         ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:2574354/3‑100 (MQ=37)
                                                                          tcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         <  1:1781259/100‑1 (MQ=255)
                                                                          tcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         <  2:444730/100‑1 (MQ=255)
                                                                           cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTgt                        >  1:3070965/1‑99 (MQ=37)
                                                                           cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:1741671/1‑100 (MQ=38)
                                                                           cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:1642467/1‑100 (MQ=255)
                                                                            gtgacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       >  1:339017/7‑100 (MQ=38)
                                                                            gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:108665/100‑1 (MQ=255)
                                                                            gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:1960849/100‑1 (MQ=255)
                                                                             tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      >  1:3255970/1‑100 (MQ=255)
                                                                             tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  1:419310/100‑1 (MQ=255)
                                                                              gACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                     >  1:822828/1‑100 (MQ=255)
                                                                                  tATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGAc                  <  1:2346382/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  1:3621551/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  2:3310599/100‑1 (MQ=255)
                                                                                   aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 <  1:3621750/100‑1 (MQ=255)
                                                                                     cAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAAt               >  1:316198/1‑100 (MQ=255)
                                                                                      aGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATa              <  2:2201179/100‑1 (MQ=255)
                                                                                       ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa             >  2:1530673/1‑100 (MQ=255)
                                                                                       ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa             >  1:1760075/1‑100 (MQ=255)
                                                                                         ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           >  1:294567/1‑100 (MQ=255)
                                                                                         ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           <  2:2304441/100‑1 (MQ=255)
                                                                                         ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGACTCATTCATAATAAACAACCATGGACAAGAAAAGTAAGATTCAATATCTACACAGACAAGAAc            >  2:3184315/1‑99 (MQ=255)
                                                                                          cccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt          <  1:2383794/100‑1 (MQ=255)
                                                                                          cccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCt          <  1:780121/100‑1 (MQ=255)
                                                                                            cTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCa        >  1:2516601/1‑100 (MQ=255)
                                                                                              gATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  1:2055827/1‑100 (MQ=255)
                                                                                              gATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  1:2006102/1‑100 (MQ=255)
                                                                                              gATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      >  2:1196049/1‑100 (MQ=255)
                                                                                               aTACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAAt     <  1:909672/100‑1 (MQ=255)
                                                                                                  |                                                                                                    
GTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  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: