Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 787,780 A→G 5.4% intergenic (‑185/+17) gpmA ← / ← galM phosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913787,7800AG 5.4% 339.9 / 6.1 107intergenic (‑185/+17)gpmA/galMphosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (44/57);  new base G (4/2);  total (48/59)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

Only 100 of 107 total aligned reads displayed.

TCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCT  >  NC_000913/787633‑787926
                                                                                                                                                   |                                                                                                                                                  
tctcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCAt                                                                                                                                                   >  3:794077/1‑149 (MQ=255)
tctcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCAt                                                                                                                                                   <  4:360467/149‑1 (MQ=255)
 ctcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATg                                                                                                                                                  >  4:28342/1‑149 (MQ=255)
  tcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGa                                                                                                                                                 >  4:321041/1‑149 (MQ=255)
  tcATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGa                                                                                                                                                 <  4:442538/149‑1 (MQ=255)
   cATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGCAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGa                                                                                                                                                 >  1:159635/1‑148 (MQ=255)
   cATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGCGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGa                                                                                                                                                 >  2:91992/1‑148 (MQ=255)
    aTTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                               <  3:248328/149‑1 (MQ=255)
     ttATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                               >  2:717947/1‑148 (MQ=255)
     ttATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGa                                                                                                                                              >  3:75714/1‑149 (MQ=255)
     ttATATTACCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGa                                                                                                                                              <  4:691075/149‑1 (MQ=255)
     ttATATCGCCGCCACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAGAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                               >  2:111829/1‑148 (MQ=255)
       atatTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGTAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACTAgg                                                                                                                                            <  4:263206/149‑1 (MQ=255)
       atatTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAgg                                                                                                                                            <  3:786355/149‑1 (MQ=255)
          ttGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGc                                                                                                                                          >  1:727706/1‑148 (MQ=255)
           tGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCg                                                                                                                                         <  2:942541/148‑1 (MQ=255)
            gCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCGTCATGAGGGctt                                                                                                                                        <  1:541047/148‑3 (MQ=255)
                 gacgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAt                                                                                                                                   <  1:221799/148‑1 (MQ=255)
                  acgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTa                                                                                                                                 <  4:110035/149‑1 (MQ=255)
                       gCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCa                                                                                                                             <  1:941998/148‑1 (MQ=255)
                       gCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAg                                                                                                                            >  4:682853/1‑149 (MQ=255)
                        cAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAGCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAg                                                                                                                            >  1:688710/1‑148 (MQ=255)
                          aCAGCAATGCTTACGCATAACTATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAACCCCTCATGACGAGGGCGTAACATTACTCAGCa                                                                                                                          >  2:341563/1‑148 (MQ=255)
                          aCAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCa                                                                                                                          >  1:245688/1‑148 (MQ=255)
                           cAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAAt                                                                                                                        >  3:703165/1‑149 (MQ=255)
                              cAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATaaa                                                                                                                     <  4:336121/149‑1 (MQ=255)
                              cAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATaa                                                                                                                      <  2:727706/148‑1 (MQ=255)
                                aTGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACt                                                                                                                   >  3:593367/1‑149 (MQ=255)
                                 tGCTTACGCATAAGCATAGAGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTAGTGATTTGAACAATCTGAGATAGAGAACTCGTGACGAGGGCGTAGCATTACTCAGCAATGAACt                                                                                                                   >  2:693795/1‑148 (MQ=255)
                                  gCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTg                                                                                                                  <  2:699396/148‑1 (MQ=255)
                                   cTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGATGTTGTGAATATTGAGTTGCAAATATGTTGATGTTTGTTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAATATTACTCAGCAATAAACTGa                                                                                                                 <  1:693795/148‑1 (MQ=255)
                                         cATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTcc                                                                                                           <  2:159635/148‑1 (MQ=255)
                                           tAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGt                                                                                                         >  2:949557/1‑148 (MQ=255)
                                                 tAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAtt                                                                                                                                  <  2:725447/117‑1 (MQ=255)
                                                 tAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAtt                                                                                                                                  >  1:725447/1‑117 (MQ=255)
                                                   gCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTGGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTgg                                                                                                 <  1:170430/148‑1 (MQ=255)
                                                   gCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                >  3:595560/1‑149 (MQ=255)
                                                    cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                <  2:899108/148‑1 (MQ=255)
                                                    cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                <  2:921983/148‑1 (MQ=255)
                                                    cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                <  2:432867/148‑1 (MQ=255)
                                                    cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTAAGTTGCAAAGATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                               <  4:75714/149‑1 (MQ=255)
                                                    cGAAAATAGGGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGa                                                                                                <  2:921040/148‑1 (MQ=255)
                                                     gAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAGTAAACTGATATTCCGTCAGGCAGCAAt                                                                                              >  4:825454/1‑149 (MQ=255)
                                                      aaaaTAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAAt                                                                                              >  1:295241/1‑148 (MQ=255)
                                                      aaaaTAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                             <  3:471651/149‑1 (MQ=255)
                                                       aaaTAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                             <  1:717947/148‑1 (MQ=255)
                                                       aaaTAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                             <  2:480909/148‑1 (MQ=255)
                                                           aGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATActct                                                                                         <  2:501716/148‑1 (MQ=255)
                                                                cgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGcc                                                                                    <  2:246841/148‑1 (MQ=255)
                                                                   aGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAgg                                                                                 <  1:111829/148‑1 (MQ=255)
                                                                              gttgtGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGc                                                                      <  2:903050/148‑1 (MQ=255)
                                                                              gttgtGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGc                                                                      <  2:245688/148‑1 (MQ=255)
                                                                               ttgtGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCa                                                                     <  2:939980/148‑1 (MQ=255)
                                                                                 gtgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTc                                                                  >  3:688391/1‑149 (MQ=255)
                                                                                  tgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGa                                                                                                                 <  1:460426/101‑1 (MQ=255)
                                                                                   gAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTcc                                                                 >  2:20119/1‑148 (MQ=255)
                                                                                    aaTATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCgg                                                               <  4:751896/149‑1 (MQ=255)
                                                                                      tatTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGtt                                                             >  3:395332/1‑149 (MQ=255)
                                                                                       atTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTg                                                            >  4:480267/1‑149 (MQ=255)
                                                                                         tGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCg                                                          <  3:473969/149‑1 (MQ=255)
                                                                                          gAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCg                                                          <  2:917181/148‑1 (MQ=255)
                                                                                             ttGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCa                                                      <  3:72302/149‑1 (MQ=255)
                                                                                               gCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCAtt                                                    <  4:595560/149‑1 (MQ=255)
                                                                                                  aaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCa                                                  >  2:898244/1‑148 (MQ=255)
                                                                                                   atatGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGaa                                                                        <  3:188755/125‑1 (MQ=255)
                                                                                                   atatGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGaa                                                                        >  4:188755/1‑125 (MQ=255)
                                                                                                    tatGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAggg                                               <  4:810/149‑1 (MQ=255)
                                                                                                     atGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAggg                                               <  1:76214/148‑1 (MQ=255)
                                                                                                      tGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtg                                             <  3:887762/149‑1 (MQ=255)
                                                                                                               ttGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCt                                    <  3:789676/149‑1 (MQ=255)
                                                                                                                  cTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACCTTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGt                                  <  2:721905/148‑1 (MQ=255)
                                                                                                                  cTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTc                                 <  3:61615/149‑1 (MQ=255)
                                                                                                                    ggTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTcc                                >  2:680491/1‑148 (MQ=255)
                                                                                                                      tGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGCGTGGTTCGGGCTGTCCGGt                             <  3:751916/149‑1 (MQ=255)
                                                                                                                       gATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTa                            >  3:666802/1‑149 (MQ=255)
                                                                                                                         tttGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAg                           >  2:297950/1‑148 (MQ=255)
                                                                                                                            gAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAAc                       >  4:10973/1‑149 (MQ=255)
                                                                                                                             aaCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                      >  3:494045/1‑149 (MQ=255)
                                                                                                                              gcAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                      <  2:688710/147‑1 (MQ=255)
                                                                                                                              aCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                      >  2:396536/1‑148 (MQ=255)
                                                                                                                              aCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                      >  2:742834/1‑148 (MQ=255)
                                                                                                                              aCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACt                      <  2:198647/148‑1 (MQ=255)
                                                                                                                              aCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTc                     <  3:6372/149‑1 (MQ=255)
                                                                                                                               cAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTc                     >  2:725587/1‑148 (MQ=255)
                                                                                                                                   atGAGATAAAGCCCTCATGACGATGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCttt                >  3:210057/1‑149 (MQ=255)
                                                                                                                                      agaTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTcc              >  2:903142/1‑148 (MQ=255)
                                                                                                                                       gaTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCag            >  3:877601/1‑149 (MQ=255)
                                                                                                                                        gttacgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCag            >  2:943787/6‑148 (MQ=255)
                                                                                                                                        gttacgCCCTCGTCAGGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGAGGAAGGAAGCAGTACGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCAGGTAGAAACTCGCTTTCCag            >  2:588518/6‑148 (MQ=255)
                                                                                                                                        aTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCaga           <  4:437866/149‑1 (MQ=255)
                                                                                                                                          tacgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCagag          >  2:432199/4‑148 (MQ=255)
                                                                                                                                           acgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagc         <  1:943787/146‑1 (MQ=255)
                                                                                                                                             gCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagccag      >  4:832089/1‑149 (MQ=255)
                                                                                                                                                cTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGccc   <  3:429715/149‑1 (MQ=255)
                                                                                                                                                cTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGccc   >  4:734378/1‑149 (MQ=255)
                                                                                                                                                cTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACGCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagccagcc    >  2:546155/1‑148 (MQ=255)
                                                                                                                                                 tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGc                                                                             <  2:372266/74‑1 (MQ=255)
                                                                                                                                                 tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGc                                                                             >  1:372266/1‑74 (MQ=255)
                                                                                                                                                 tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt  >  4:292722/1‑149 (MQ=255)
                                                                                                                                                 tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt  >  3:680564/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCT  >  NC_000913/787633‑787926

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: