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.

CAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  NC_000913/2449525‑2449723
                                                                                                    |                                                                                                   
   ttttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTc                                                                                                  <  1:136591/100‑1 (MQ=255)
    tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                                 <  1:789232/100‑1 (MQ=255)
    tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCa                                                                                                 <  1:1528732/100‑1 (MQ=255)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACg                                                                                               <  1:2547498/100‑1 (MQ=38)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACg                                                                                               <  1:1735186/100‑1 (MQ=38)
      tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACg                                                                                               <  2:395467/100‑1 (MQ=38)
       gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGc                                                                                              >  2:2943417/1‑100 (MQ=255)
       gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGc                                                                                              <  1:2678171/100‑1 (MQ=255)
          ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                           <  1:6712/100‑5 (MQ=25)
          ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATAGTCACgacac                                                                                           <  1:3504388/100‑5 (MQ=25)
             cTGATAAAATTCTCAATCCAGGAAAGGTCGTAAGGTTGAACAGAAACGGATGGACAAATCGAGTCGGGACGTATCAGGGGCCTGATAGTCACgacac                                                                                           <  2:2327055/97‑5 (MQ=16)
              tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  1:1204780/100‑1 (MQ=255)
              tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  2:2639178/100‑1 (MQ=255)
              tGACATAATGCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                                       <  2:2759161/100‑1 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:1804157/1‑100 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:2519122/1‑100 (MQ=255)
               gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                                      >  1:993415/1‑100 (MQ=255)
                aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                                     <  1:1969159/100‑1 (MQ=255)
                 cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                                    >  2:837083/1‑100 (MQ=255)
                  aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                   <  2:549892/100‑1 (MQ=255)
                  aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                                   <  2:1667373/100‑1 (MQ=255)
                   tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                                  >  1:2881204/1‑100 (MQ=255)
                    aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                 <  1:281132/100‑1 (MQ=255)
                    aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                                                 >  2:344106/1‑100 (MQ=255)
                      ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                                               >  1:767145/1‑100 (MQ=255)
                          cATTCCATGAAACGTCATAAGGTGGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCa                                                                             <  1:3090646/98‑1 (MQ=255)
                           aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGCGAGGCTTCagag                                                                          >  1:1139014/1‑100 (MQ=255)
                            ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                                         >  1:2048234/1‑100 (MQ=255)
                                aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGCGGCTTCAGAGtcatt                                                                     >  1:1459038/1‑100 (MQ=255)
                                aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcatt                                                                     >  1:2443249/1‑100 (MQ=255)
                                 tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                                    >  2:378735/1‑100 (MQ=255)
                                   aaaCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattcat                                                                  <  1:591116/100‑1 (MQ=255)
                                       gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                                              <  1:3544301/100‑1 (MQ=255)
                                        tCATAAGGTTGAACAGAAATGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                                             >  2:2399574/1‑100 (MQ=255)
                                         cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                                            <  1:790392/100‑1 (MQ=255)
                                         cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATaac                                                            <  2:2443236/100‑1 (MQ=255)
                                           tAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacca                                                          >  2:2390220/1‑100 (MQ=255)
                                             agcttGAACAGAACCGGACGTACAAATCGTGTCCTGACTATCAGGCCCCTGATACCTCACGCCCCCTGGAGACGCTTCAGAGTCATTCATAATaac                                                            <  2:2547191/93‑1 (MQ=255)
                                               gTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaacca                                                      >  1:3903357/1‑100 (MQ=255)
                                                ttGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAGCCAt                                                     <  2:1069352/100‑1 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    <  2:3825569/100‑1 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    >  2:3061869/1‑100 (MQ=255)
                                                 tGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                                    <  2:2161353/100‑1 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   <  2:3880130/100‑1 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   >  1:2108499/1‑100 (MQ=255)
                                                  gAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                                   <  1:508995/100‑1 (MQ=255)
                                                   aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                                  <  1:186034/100‑1 (MQ=255)
                                                   aaCAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTCATAATAACCCACCATGGa                                                  >  2:1963663/1‑100 (MQ=255)
                                                       gAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAg                                              >  1:2938214/1‑100 (MQ=255)
                                                         aaCGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGa                                             <  2:3809782/99‑1 (MQ=255)
                                                           cGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                                          >  1:1733855/1‑100 (MQ=255)
                                                            ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  2:1094206/100‑1 (MQ=255)
                                                            ggATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                                         <  1:922304/100‑1 (MQ=255)
                                                              aTGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTa                                       >  1:1602642/1‑100 (MQ=255)
                                                                ggACAAATCGTGTCGTTACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     <  2:1534791/100‑1 (MQ=255)
                                                                ggACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                                     >  1:2164761/1‑100 (MQ=255)
                                                                 gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTTCATAATCACCAACCATGGACACAAAAAGTagaa                                    >  2:3528040/1‑97 (MQ=255)
                                                                 gACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGa                                    >  1:2582908/1‑100 (MQ=255)
                                                                  aCAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAt                                   <  1:3475061/100‑1 (MQ=255)
                                                                   cAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                                  <  1:1337307/100‑1 (MQ=255)
                                                                    aaaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTc                                 <  2:1891849/100‑1 (MQ=255)
                                                                     aaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                <  2:2582891/100‑1 (MQ=255)
                                                                     aaTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa                                >  1:102483/1‑100 (MQ=255)
                                                                      aTCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTAATTCCTAATAACCACCAAGGGACAACAAAAGTAAGATCCaa                               >  2:659508/1‑100 (MQ=255)
                                                                        cgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  1:287470/100‑1 (MQ=255)
                                                                        cgtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc                             <  2:1257433/100‑1 (MQ=255)
                                                                         gtgtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                            >  1:2872612/1‑100 (MQ=255)
                                                                         ggggcgtgACGTATCAGGGGCCTGATAGTCACACCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct                            >  1:3743833/5‑100 (MQ=25)
                                                                           gtcgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          <  2:100897/100‑1 (MQ=255)
                                                                           ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3115337/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:767622/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3613679/3‑100 (MQ=37)
                                                                           ggcgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct                          >  1:3029200/3‑100 (MQ=37)
                                                                            gcgtgACGTATCAGGGGCCTGATAGTTACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg                         >  2:2087197/2‑100 (MQ=38)
                                                                             cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  2:458449/1‑100 (MQ=38)
                                                                             cgtgACGTATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:3622871/1‑100 (MQ=38)
                                                                             cgtgACTCTGAGGCCCCTGATACGTCACGACTCCTAGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACCAGCAAAGTAAGATTCAATCTCTGc                        <  2:3509094/100‑1 (MQ=255)
                                                                             cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        <  1:1233138/100‑1 (MQ=255)
                                                                             cgtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc                        >  1:2821126/1‑100 (MQ=255)
                                                                              gtgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       <  1:2390232/100‑1 (MQ=255)
                                                                               tgacgtATCAGGGGCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa                       >  1:929948/6‑100 (MQ=255)
                                                                               tgacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      >  1:1268517/6‑100 (MQ=38)
                                                                               tgACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  2:911788/100‑1 (MQ=255)
                                                                                gacgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga                     >  1:2302986/5‑100 (MQ=38)
                                                                                 acgtATCAGGGGCCTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCagag                    >  2:3199402/4‑100 (MQ=255)
                                                                                 aCTCTCCGGCCCCTGATACGTCACGACCCCTGGAGAGGCCTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag                      <  1:3144647/98‑1 (MQ=255)
                                                                                 aCTATCAGGACCCTGATACGTCACGCCACCTGGAGAGGCCTCATAGCCATTCATAATAATCAACCATGGACAAGAAAAGTAAGATGCAATCTCTGCagag                    <  2:416174/100‑1 (MQ=255)
                                                                                     aTCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACa                 >  2:1231097/1‑100 (MQ=255)
                                                                                      tcaggggccTGATAGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACaa                >  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)
                                                                                                    |                                                                                                   
CAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACTATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAGACAATAACCTCAAATGTA  >  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: