Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,227 A→T 15.9% intergenic (‑75/+26) btuC ← / ← ihfA vitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,795,2270AT15.9% 229.6 / 31.9 93intergenic (‑75/+26)btuC/ihfAvitamin B12 ABC transporter permease/integration host factor (IHF), DNA‑binding protein, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (42/36);  new base T (9/6);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGC  >  NC_000913/1795128‑1795325
                                                                                                   |                                                                                                  
gttgttGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCa                                                                                                    >  1:436534/1‑100 (MQ=255)
   gttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCagag                                                                                                 >  1:1181862/1‑100 (MQ=255)
   gttgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCagag                                                                                                 >  2:2579722/1‑100 (MQ=255)
    ttgGCGGGTAAGTGTCAGCATGGTATTCGTTTTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCagatc                                                                                                <  1:2049111/100‑3 (MQ=255)
     tgGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCg                                                                                               <  1:2493835/100‑1 (MQ=255)
      ggCGGGCAAGTGTGAGAATGGTATCCGTGCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGAAAAAAAAAAAAAAAAAGGCCGCAGAGCgg                                                                                              <  2:2405154/100‑1 (MQ=255)
      ggCGGGCAAGTGTCAGCATTGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCgg                                                                                              <  1:2905128/100‑1 (MQ=255)
      ggCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCgg                                                                                              <  2:655855/100‑1 (MQ=255)
        cGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGcc                                                                                            <  1:2579754/100‑1 (MQ=255)
           gCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGAATCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCggaattt                                                                                         >  1:1002497/1‑95 (MQ=255)
             aaGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttt                                                                                       <  1:2647827/100‑1 (MQ=255)
              aGTGTCAGCATGGTATCCGTTCTTCTGAATTTTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCtttttt                                                                                      <  1:974878/100‑2 (MQ=38)
              aGTGTAACCATGGTATCCGTTCTGCCGAAGTGTCATGGAGTTGATTTTACGGTGACTCTTCGACAGTGAAAACAAAAAAGGCCGAAGAGCGGCCttttt                                                                                       <  2:160724/99‑1 (MQ=255)
               gtgtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGTGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                     <  1:1301727/100‑1 (MQ=255)
               gtgtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGAGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                     <  1:2125992/100‑1 (MQ=255)
               gtgtCAGCATGGTATCCGTTCTGCTGAAGAGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTTAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttag                                                                                     <  2:1823696/100‑1 (MQ=255)
               gtgtCAGAATGGTATCCGTTCTGCGGAAGTGTGATGGGGTTGATTTTACGGTGCCTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTtta                                                                                      <  2:2609354/99‑1 (MQ=255)
                 gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGTTTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                   <  1:1238375/100‑1 (MQ=255)
                 gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                   >  1:1828703/1‑100 (MQ=255)
                 gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                   <  1:360897/100‑1 (MQ=255)
                 gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                   <  2:1829564/100‑1 (MQ=255)
                 gtCAGCATGCTATCCGTTTTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCAGCTCTGCGGCCttttttct                                                                                    <  2:1193887/99‑4 (MQ=38)
                  tCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttcttn                                                                                  <  2:784076/100‑6 (MQ=38)
                   cAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCTCTGCGGCCttttttctttt                                                                                 <  1:2272659/100‑7 (MQ=37)
                   cAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagttag                                                                                 >  2:1378862/1‑100 (MQ=255)
                    aGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGTCGCAGAGCGGCCTTTttagtt                                                                                   <  2:1592074/97‑1 (MQ=255)
                    aGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGa                                                                                >  2:43380/1‑100 (MQ=255)
                     gCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGAt                                                                               >  1:2647137/1‑100 (MQ=255)
                      cATGGTATCCGTTCTGCTGAAGTGTCATTGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                              <  2:1443304/100‑1 (MQ=255)
                      cATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATc                                                                              <  1:1315863/100‑1 (MQ=255)
                          gTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                          <  1:2681905/100‑1 (MQ=255)
                           tATCCGTTTTGCTGAATTGTCTTGGCTTTGATTTTACGGTGACTCTTCGACATTGCAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAt                                                                          <  1:2681989/99‑1 (MQ=255)
                             tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTAc                                                                       <  2:1706692/100‑1 (MQ=255)
                             tCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTAc                                                                       >  2:2009188/1‑100 (MQ=255)
                              ccGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGACCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACt                                                                      >  1:930030/1‑100 (MQ=255)
                               cGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                     >  1:2134092/1‑100 (MQ=255)
                                gTTCTGCTGAAGTGTCATGGAGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTc                                                                     >  2:946336/1‑99 (MQ=255)
                                 ttCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCg                                                                    <  2:2861116/99‑1 (MQ=255)
                                   ctgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTc                                                                  >  1:2155659/1‑99 (MQ=255)
                                    tgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCtt                                                                >  2:703718/1‑100 (MQ=255)
                                    tgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCtt                                                                >  2:2019520/1‑100 (MQ=255)
                                    tgctgAAGTGTCATGACGTTGATTATACGGTGACTCTTAGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGAGTACACGTCtt                                                                >  2:211486/1‑100 (MQ=255)
                                     gctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTAGTCttt                                                               >  1:1440941/1‑100 (MQ=255)
                                     gctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTATAGTTAGATCAGATTACTCGTCttt                                                               >  1:282133/1‑100 (MQ=255)
                                       tgAAGTGTCATGGCGTTGATTTTACGGTGCCTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                             >  1:850966/1‑100 (MQ=255)
                                       tgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                             <  1:43380/100‑1 (MQ=255)
                                       tgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTgg                                                             <  2:2025881/100‑1 (MQ=255)
                                        gAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                            <  1:2985801/100‑1 (MQ=255)
                                        gAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                            >  1:462762/1‑100 (MQ=255)
                                             gtCCTGGCGTTGATTTTACGGTGACTCTTCGACCGTGAAAAGAAAAAAGGCCGAAGAGCGGCCCTGTTACTTAGATCAGATTACTCGTATTTAGGCGa                                                         >  2:2779582/1‑98 (MQ=255)
                                             gtCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAg                                                       >  2:2764872/1‑100 (MQ=255)
                                                aTGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGt                                                    >  2:1979170/1‑100 (MQ=255)
                                                aTGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTAGTCTTTGGGCGAAGCg                                                     >  2:1684386/1‑99 (MQ=255)
                                                 tGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGtt                                                   >  1:1583688/1‑100 (MQ=255)
                                                   gCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGtttt                                                 <  2:159321/100‑1 (MQ=255)
                                                   gCGGTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGGCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGGGtttt                                                 <  2:687615/100‑1 (MQ=255)
                                                    cGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                >  1:2574089/1‑100 (MQ=255)
                                                    cGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTc                                                >  2:946047/1‑100 (MQ=255)
                                                    cGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAATTAGATCAGATAACTCATATTTGGACGAAGCGTTTTc                                                >  1:1156550/1‑100 (MQ=255)
                                                     gTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCg                                               >  1:801812/1‑100 (MQ=255)
                                                       tGATTGTACGTGTACTCTTCGTCACTGTAAAAAAAAAAGGCCGCTGAGCTGCCTTTTTAGTTAGTTCATATTTCTCGTCTTTGCGCGAAGCGTTTTc                                                <  2:2712456/97‑1 (MQ=255)
                                                         agtttACGGTGACACTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGAccc                                           >  1:2057394/3‑100 (MQ=255)
                                                         aTTTTAGGGTTACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGTGCGGTCTTTTTAGTTCGATCAGTTTACTCGTCTTTGGGCGAAGCGTTTTCGAcct                                           <  1:2888118/100‑2 (MQ=255)
                                                          ttttACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGAccc                                           >  1:2192632/1‑99 (MQ=255)
                                                          ttttACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCg                                          >  1:2107538/1‑100 (MQ=255)
                                                          ttttACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCg                                          >  2:777822/1‑100 (MQ=255)
                                                               cGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCttt                                     >  2:253175/1‑100 (MQ=255)
                                                                 gTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCttttt                                   <  2:2174682/100‑1 (MQ=255)
                                                                  tGACTCTTCGACAGTTAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTa                                  <  1:2579761/100‑1 (MQ=255)
                                                                    aCTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTaa                                 >  2:304495/1‑99 (MQ=255)
                                                                     ctctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACt                               >  2:1889249/1‑100 (MQ=255)
                                                                     ctctTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTAGTCTTTGGGCGAAGCGTTTTCGACACGGCTTTTTACCt                               >  2:477873/1‑100 (MQ=255)
                                                                        ttCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTc                             <  1:2204298/99‑1 (MQ=255)
                                                                        ttCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCt                            <  2:2523543/100‑1 (MQ=255)
                                                                          cGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGc                          >  2:2892583/1‑100 (MQ=255)
                                                                           gACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGcc                         <  1:2110467/100‑1 (MQ=255)
                                                                            aCAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGAGCGAAGAGTTTTAGACACGGCTTTTTAACTTATGctc                        >  1:2979/1‑98 (MQ=255)
                                                                             cAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCg                       >  1:202015/1‑100 (MQ=255)
                                                                               gTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                     <  1:1979193/100‑1 (MQ=255)
                                                                               gTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                     <  2:813923/100‑1 (MQ=255)
                                                                                tGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGAGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGt                    >  1:973957/1‑100 (MQ=255)
                                                                                 gaaaagaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTc                   <  1:177472/100‑1 (MQ=255)
                                                                                     actaaaaaGGCCGCTCTTCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTGTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCt                  >  2:638977/4‑97 (MQ=38)
                                                                                      gaaaaAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGaa               >  1:2159161/1‑99 (MQ=255)
                                                                                      ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGaa               >  1:2164809/3‑99 (MQ=38)
                                                                                      ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg              >  2:157829/3‑100 (MQ=39)
                                                                                      ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg              >  2:1740419/3‑100 (MQ=38)
                                                                                      ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAACGTTATCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg              >  1:81179/3‑100 (MQ=38)
                                                                                       taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg             >  1:1575260/2‑100 (MQ=39)
                                                                                       taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg             >  1:1972292/2‑100 (MQ=39)
                                                                                       aaaaaaTGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGTCCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg              <  1:2886768/99‑1 (MQ=255)
                                                                                       aaaaaaGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg             >  1:1517449/1‑100 (MQ=255)
                                                                                        aaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGt            >  2:2317985/1‑100 (MQ=255)
                                                                                        aaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGAGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGt            >  2:75286/1‑100 (MQ=255)
                                                                                              gCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCTTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTcaccac      >  2:1228082/1‑100 (MQ=255)
                                                                                                cGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGGCTTTGGGCGAAGCGTTTTCGACCCGTCTTTTTACCTTCTGCCCCGGTCTGAAGGTcaccac      >  2:1071421/1‑98 (MQ=255)
                                                                                                  cAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACgcgc  <  2:1934720/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GTTGTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGC  >  NC_000913/1795128‑1795325

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