Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,449,607:1 +G 15.0% coding (443/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,6071.G15.0% 573.9 / 103.8 160coding (443/822 nt)yfcODUF2544 family putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (62/74);  new base G (12/12);  total (74/86)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

Only 100 of 162 total aligned reads displayed.

GTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAG  >  NC_000913/2449508‑2449705
                                                                                                    |                                                                                                  
 tCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑t                                                                                                   <  2:1277207/100‑1 (MQ=255)
    aCTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATc                                                                                                <  1:139625/100‑1 (MQ=255)
     cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCa                                                                                               <  2:916961/100‑1 (MQ=255)
     cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCa                                                                                               >  1:2233748/1‑100 (MQ=255)
     cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCa                                                                                               <  2:71031/100‑1 (MQ=255)
     cTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCa                                                                                               <  1:1313781/100‑1 (MQ=255)
      ttgaGGGTGAGCTGGTGTGACTTACTGATATAATTCGCAGTCTATGAAACGTCTTAAGGTTGAACAGAAACGGTTGGTCAAATCGTGTCGTGAC‑TATCAg                                                                                              <  2:3396443/97‑1 (MQ=255)
       gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAgg                                                                                             >  1:100897/1‑100 (MQ=255)
       gCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAgg                                                                                             >  2:416762/1‑100 (MQ=255)
         aGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcc                                                                                           >  2:2678154/1‑100 (MQ=255)
           ggTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGcccc                                                                                         <  1:901430/100‑1 (MQ=255)
             tGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTg                                                                                       <  2:123847/100‑1 (MQ=255)
              gAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGa                                                                                      <  1:3348116/100‑1 (MQ=255)
              gAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGa                                                                                      <  2:232881/100‑1 (MQ=255)
               aGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGAt                                                                                     <  1:1771475/100‑1 (MQ=255)
                gCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                    >  1:2203173/1‑100 (MQ=255)
                 cAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATAc                                                                                   <  1:3723939/100‑1 (MQ=255)
                  aGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATa                                                                                    <  1:1497804/98‑1 (MQ=255)
                     tttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCa                                                                               <  1:1528732/100‑1 (MQ=255)
                      ttGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCAc                                                                              >  2:236923/1‑100 (MQ=255)
                       tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                             <  2:2500783/100‑1 (MQ=38)
                       tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                             <  1:2547498/100‑1 (MQ=38)
                       tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACg                                                                             <  2:395467/100‑1 (MQ=38)
                       tGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACg                                                                             <  1:3517187/100‑1 (MQ=255)
                        gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                            <  1:2678171/100‑1 (MQ=255)
                        gACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGc                                                                            >  2:2943417/1‑100 (MQ=255)
                         aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGTATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgac                                                                           <  2:1756521/100‑3 (MQ=37)
                         aCTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgac                                                                           <  2:2616484/100‑3 (MQ=37)
                          cTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgaca                                                                          <  1:2853232/100‑4 (MQ=37)
                           ttACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGACGTATCAGGGGCCTGATA‑GTCACgacac                                                                         <  1:6712/100‑5 (MQ=25)
                              cTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTg                                                                      >  1:780538/1‑100 (MQ=255)
                               tGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTgg                                                                     <  1:1204780/100‑1 (MQ=255)
                                gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                    >  1:993415/1‑100 (MQ=255)
                                gACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGga                                                                    >  1:1239652/1‑100 (MQ=255)
                                 aCATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgag                                                                   <  1:1969159/100‑1 (MQ=255)
                                  cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                  >  2:837083/1‑100 (MQ=255)
                                  cATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                  >  2:591113/1‑100 (MQ=255)
                                  cATAATTCTCATGCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGgaga                                                                  >  1:2602716/1‑100 (MQ=255)
                                   aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                 <  2:1667373/100‑1 (MQ=255)
                                   aTAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGagag                                                                 <  1:1654804/100‑1 (MQ=255)
                                    tAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAgg                                                                >  1:2881204/1‑100 (MQ=255)
                                     aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                               >  2:344106/1‑100 (MQ=255)
                                     aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                               <  1:3461728/100‑1 (MQ=255)
                                     aaTTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGc                                                               <  1:281132/100‑1 (MQ=255)
                                       ttCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                             >  1:767145/1‑100 (MQ=255)
                                         ctcATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCtt                                                             >  2:113553/1‑98 (MQ=255)
                                            aTTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGCGAGGCTTCagag                                                        >  1:1139014/1‑100 (MQ=255)
                                             ttCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGt                                                       >  1:2048234/1‑100 (MQ=255)
                                               ccATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtca                                                     >  1:2513922/1‑100 (MQ=255)
                                                 aTGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGCGGCTTCAGAGtcatt                                                   >  1:1459038/1‑100 (MQ=255)
                                                  tGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATAAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGtcattc                                                  >  2:378735/1‑100 (MQ=255)
                                                      aCGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                              >  1:1159302/1‑100 (MQ=255)
                                                       cGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataa                                              <  2:416019/99‑1 (MQ=255)
                                                        gTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCataata                                            <  1:3544301/100‑1 (MQ=255)
                                                         tCATAAGGTTGAACAGAAATGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCAtaataa                                           >  2:2399574/1‑100 (MQ=255)
                                                         tCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTTACGCCCCCTGGAGAGGCTTCACAGTCATTCataata                                            >  1:3880160/1‑99 (MQ=255)
                                                          cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaac                                          <  1:790392/100‑1 (MQ=255)
                                                          cATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGGCATTCATAATaac                                          <  2:2443236/100‑1 (MQ=255)
                                                           aTAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATaacc                                         <  2:1045733/100‑1 (MQ=255)
                                                                gTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAaccaacca                                    >  1:3903357/1‑100 (MQ=255)
                                                                  tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                  >  2:3061869/1‑100 (MQ=255)
                                                                  tGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCACCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATg                                  <  2:2161353/100‑1 (MQ=255)
                                                                   gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                 >  1:928839/1‑100 (MQ=255)
                                                                   gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                 <  1:508995/100‑1 (MQ=255)
                                                                   gAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATgg                                 <  2:3880130/100‑1 (MQ=255)
                                                                    aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGa                                <  1:186034/100‑1 (MQ=255)
                                                                    aaCAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTGCAGAGTCATTCATAATAACCCACCATGGa                                >  2:1963663/1‑100 (MQ=255)
                                                                    aaCAGAAACGGATGGACAAATCGTGACGTGAC‑TAGCAGGCCCCTGATACGTGACGCCCCCTGGAGAGGCTTCAGCGTCATTAATAATAACTAGCCCTGGa                                >  1:662334/1‑100 (MQ=255)
                                                                      cAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACa                              <  1:2154373/100‑1 (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)
                                                                            cGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGaaaa                        >  1:3269535/1‑100 (MQ=255)
                                                                             ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                       <  2:1094206/100‑1 (MQ=255)
                                                                             ggATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAg                       <  1:922304/100‑1 (MQ=255)
                                                                                 ggACAAATCGTGTCGTTAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAg                   <  2:1534791/100‑1 (MQ=255)
                                                                                  gACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCGTTCATAATCACCAACCATGGACACAAAAAGTagaa                  >  2:3528040/1‑97 (MQ=255)
                                                                                    cAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGAtt                <  1:1337307/100‑1 (MQ=255)
                                                                                      aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa              >  1:102483/1‑100 (MQ=255)
                                                                                      aaTCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCa              <  2:2582891/100‑1 (MQ=255)
                                                                                         cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc           <  2:1257433/100‑1 (MQ=255)
                                                                                         cgtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtc           <  1:236923/100‑1 (MQ=255)
                                                                                          gtgtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct          >  1:2872612/1‑100 (MQ=255)
                                                                                          ggggcgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAAtct          >  1:311525/5‑100 (MQ=25)
                                                                                            gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct        <  2:100897/100‑1 (MQ=255)
                                                                                            gtcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATctct        >  2:3517687/1‑100 (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)
                                                                                             tcgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTg       <  2:2313486/100‑1 (MQ=255)
                                                                                              cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc      >  1:3622871/1‑100 (MQ=38)
                                                                                              cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc      >  1:3479491/1‑100 (MQ=38)
                                                                                              cgtgACGTATCAGGGGCCTGATA‑GTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc      >  2:458449/1‑100 (MQ=38)
                                                                                              cgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc      >  2:322933/1‑100 (MQ=255)
                                                                                              cgtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGc      <  1:1233138/100‑1 (MQ=255)
                                                                                               gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa     <  1:2390232/100‑1 (MQ=255)
                                                                                               gtgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCa     >  1:1746097/1‑100 (MQ=255)
                                                                                                tgAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCag    <  2:3818553/100‑1 (MQ=255)
                                                                                                 gAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCaga   >  1:252097/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
GTCGACTGCAGGGTGAGCAGTTTTGACTTACTGACATAATTCTCATTCCATGAAACGTCATAAGGTTGAACAGAAACGGATGGACAAATCGTGTCGTGAC‑TATCAGGCCCCTGATACGTCACGCCCCCTGGAGAGGCTTCAGAGTCATTCATAATAACCAACCATGGACAAGAAAAGTAAGATTCAATCTCTGCAGAG  >  NC_000913/2449508‑2449705

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