Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,795,229 A→T 15.8% intergenic (‑77/+24) 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,2290AT15.8% 227.9 / 31.6 93intergenic (‑77/+24)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 = 1.04e-01

GTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCG  >  NC_000913/1795131‑1795327
                                                                                                  |                                                                                                  
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                                                                                     <  2:1829564/100‑1 (MQ=255)
              gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                     <  1:360897/100‑1 (MQ=255)
              gtCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTttagtt                                                                                     >  1:1828703/1‑100 (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:2019520/1‑100 (MQ=255)
                                 tgctgAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCtt                                                                  >  2:703718/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:462762/1‑100 (MQ=255)
                                     gAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTggg                                                              <  1:2985801/100‑1 (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                       <  2:813923/100‑1 (MQ=255)
                                                                            gTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCggg                       <  1:1979193/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:1740419/3‑100 (MQ=38)
                                                                                   ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                >  2:157829/3‑100 (MQ=39)
                                                                                   ctaaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAACGTTATCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAg                >  1:81179/3‑100 (MQ=38)
                                                                                    taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg               >  1:1972292/2‑100 (MQ=39)
                                                                                    taaaaaGGCCGCTCTGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAgg               >  1:1575260/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)
                                                                                                 gagCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGTAGCTTTTTTGACCCGGTGTTTTAACTTATGCCCGGGTCTGAAGGTCACCACGCGCCg  <  2:898118/100‑1 (MQ=255)
                                                                                                  agCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGTGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCACGGGTCTGAAGGTAACAACGCGCCg  >  2:1548970/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
GTTGGCGGGCAAGTGTCAGCATGGTATCCGTTCTGCTGAAGTGTCATGGCGTTGATTTTACGGTGACTCTTCGACAGTGAAAAGAAAAAAGGCCGCAGAGCGGCCTTTTTAGTTAGATCAGATTACTCGTCTTTGGGCGAAGCGTTTTCGACCCGGCTTTTTAACTTCTGCCCGGGTCTGAAGGTCACCACGCGCCG  >  NC_000913/1795131‑1795327

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: