Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,624 Δ1 bp 18.1% 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.18.1% 470.1 / 105.9 160coding (426/822 nt)yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (68/63);  new base . (17/12);  total (85/75)
Fisher's exact test for biased strand distribution p-value = 5.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 163 total aligned reads displayed.

CAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449525‑2449723
                                                                                                    |                                                                                                   
   ttttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTc                                                                                                  <  1:136591/100‑1 (MQ=255)
    tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                                 <  1:789232/100‑1 (MQ=255)
    tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                                                 <  1:1528732/100‑1 (MQ=255)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                               <  1:2547498/100‑1 (MQ=38)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                               <  1:1735186/100‑1 (MQ=38)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                                               <  2:395467/100‑1 (MQ=38)
       gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                                              >  2:2943417/1‑100 (MQ=255)
       gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                                              <  1:2678171/100‑1 (MQ=255)
          ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                           <  1:6712/100‑5 (MQ=25)
          ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                           <  1:3504388/100‑5 (MQ=25)
             cTGATAAAATTCTCAATCCAGGAAAGGTCGTAAGGTTGAACAGAAACGGATGGACAAATCGAGTCGGGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                                           <  2:2327055/97‑5 (MQ=16)
              tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  1:1204780/100‑1 (MQ=255)
              tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  2:2639178/100‑1 (MQ=255)
              tGACATAATGCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  2:2759161/100‑1 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:1804157/1‑100 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:2519122/1‑100 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:993415/1‑100 (MQ=255)
                aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                     <  1:1969159/100‑1 (MQ=255)
                 cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                                    >  2:837083/1‑100 (MQ=255)
                  aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                   <  2:549892/100‑1 (MQ=255)
                  aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                   <  2:1667373/100‑1 (MQ=255)
                   tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                  >  1:2881204/1‑100 (MQ=255)
                    aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                 <  1:281132/100‑1 (MQ=255)
                    aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                 >  2:344106/1‑100 (MQ=255)
                      ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                               >  1:767145/1‑100 (MQ=255)
                          cATTCCATGAAACGTCATAAGGTGGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                             <  1:3090646/98‑1 (MQ=255)
                           aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGCGAGGCTTCagag                                                                          >  1:1139014/1‑100 (MQ=255)
                            ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                         >  1:2048234/1‑100 (MQ=255)
                                aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGCGGCTTCAGAGtcatt                                                                     >  1:1459038/1‑100 (MQ=255)
                                aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                     >  1:2443249/1‑100 (MQ=255)
                                 tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                    >  2:378735/1‑100 (MQ=255)
                                   aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  1:591116/100‑1 (MQ=255)
                                       gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                                              <  1:3544301/100‑1 (MQ=255)
                                        tCATAAGGTTGAACAGAAATGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                             >  2:2399574/1‑100 (MQ=255)
                                         cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  1:790392/100‑1 (MQ=255)
                                         cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATaac                                                            <  2:2443236/100‑1 (MQ=255)
                                           tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                          >  2:2390220/1‑100 (MQ=255)
                                             agcttGAACAGAACCGGACGTACAAATCGTGTCCTGAC‑TATCAGGCCCCTGATACCTCACGCCCCCTGGAGACGCTTCAGAGTCATTCATAATaac                                                            <  2:2547191/93‑1 (MQ=255)
                                               gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaacca                                                      >  1:3903357/1‑100 (MQ=255)
                                                ttGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAGCCAt                                                     <  2:1069352/100‑1 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    <  2:3825569/100‑1 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    >  2:3061869/1‑100 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    <  2:2161353/100‑1 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   <  2:3880130/100‑1 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   >  1:2108499/1‑100 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   <  1:508995/100‑1 (MQ=255)
                                                   aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                  <  1:186034/100‑1 (MQ=255)
                                                   aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTCATAATAACCCACCATGGa                                                  >  2:1963663/1‑100 (MQ=255)
                                                       gAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                              >  1:2938214/1‑100 (MQ=255)
                                                         aaCGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGa                                             <  2:3809782/99‑1 (MQ=255)
                                                           cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                          >  1:1733855/1‑100 (MQ=255)
                                                            ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:1094206/100‑1 (MQ=255)
                                                            ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  1:922304/100‑1 (MQ=255)
                                                              aTGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                                       >  1:1602642/1‑100 (MQ=255)
                                                                ggACAAATCGTGTCGTTAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     <  2:1534791/100‑1 (MQ=255)
                                                                ggACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     >  1:2164761/1‑100 (MQ=255)
                                                                 gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTTCATAATCACCAACCATGGACACAAAAAGTagaa                                    >  2:3528040/1‑97 (MQ=255)
                                                                 gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                    >  1:2582908/1‑100 (MQ=255)
                                                                  aCAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                                   <  1:3475061/100‑1 (MQ=255)
                                                                   cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                  <  1:1337307/100‑1 (MQ=255)
                                                                    aaaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTc                                 <  2:1891849/100‑1 (MQ=255)
                                                                     aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                <  2:2582891/100‑1 (MQ=255)
                                                                     aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                >  1:102483/1‑100 (MQ=255)
                                                                      aTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTAATTCCTAATAACCACCAAGGGACAACAAAAGTAAGATCCaa                               >  2:659508/1‑100 (MQ=255)
                                                                        cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  1:287470/100‑1 (MQ=255)
                                                                        cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  2:1257433/100‑1 (MQ=255)
                                                                         gtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                            >  1:2872612/1‑100 (MQ=255)
                                                                         ggggcgtgACGTATCAGGGGCCTGATA‑GTCACACCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                            >  1:3743833/5‑100 (MQ=25)
                                                                           gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          <  2:100897/100‑1 (MQ=255)
                                                                           ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3115337/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:767622/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3613679/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3029200/3‑100 (MQ=37)
                                                                            gcgtgACGTATCAGGGGCCTGATA‑GTTACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         >  2:2087197/2‑100 (MQ=38)
                                                                             cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  2:458449/1‑100 (MQ=38)
                                                                             cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:3622871/1‑100 (MQ=38)
                                                                             cgtgAC‑TCTGAGGCCCCTGATACGTCACGACTCCTAGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACCAGCAAAGTAAGATTCAATCTCTGc                        <  2:3509094/100‑1 (MQ=255)
                                                                             cgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        <  1:1233138/100‑1 (MQ=255)
                                                                             cgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:2821126/1‑100 (MQ=255)
                                                                              gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:2390232/100‑1 (MQ=255)
                                                                               tgacgtATCAGGGGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       >  1:929948/6‑100 (MQ=255)
                                                                               tgacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      >  1:1268517/6‑100 (MQ=38)
                                                                               tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  2:911788/100‑1 (MQ=255)
                                                                                gacgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                     >  1:2302986/5‑100 (MQ=38)
                                                                                 acgtATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                    >  2:3199402/4‑100 (MQ=255)
                                                                                 aC‑TCTCCGGCCCCTGATACGTCACGACCCCTGGAGAGGCCTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  1:3144647/98‑1 (MQ=255)
                                                                                 aC‑TATCAGGACCCTGATACGTCACGCCACCTGGAGAGGCCTCATAGCCATTCATAATAATCAACCATGGACAAGAAAAGTAAGATGCAATCTCTGCagag                    <  2:416174/100‑1 (MQ=255)
                                                                                     aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 >  2:1231097/1‑100 (MQ=255)
                                                                                      tcaggggccTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                >  2:2710873/8‑99 (MQ=255)
                                                                                      tCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                >  1:3470469/1‑100 (MQ=255)
                                                                                      tCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                >  1:3348300/1‑100 (MQ=255)
                                                                                         ggCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATaa             >  1:1880833/1‑100 (MQ=255)
                                                                                           ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           >  1:2349936/1‑100 (MQ=255)
                                                                                           ccccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAAcc           <  1:387598/100‑1 (MQ=255)
                                                                                            cccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAGAAGTAAGGTTCAATCTCTGCAGAGACAATAACCt          >  2:1679282/1‑100 (MQ=255)
                                                                                             ccTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTc         >  1:3538338/1‑100 (MQ=255)
                                                                                               tGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaa       <  1:3031087/100‑1 (MQ=255)
                                                                                                gATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCaaa      <  1:2943435/100‑1 (MQ=255)
                                                                                                  tACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATg    <  1:3715621/100‑1 (MQ=255)
                                                                                                   aCGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGt   <  2:1168778/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
CAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449525‑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: ‑