Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 787,780 A→G 6.8% 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 6.8% 272.4 / 5.1 86intergenic (‑185/+17)gpmA/galMphosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (37/43);  new base G (2/4);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTT  >  NC_000913/787632‑787927
                                                                                                                                                    |                                                                                                                                                   
ttCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCa                                                                                                                                                     <  3:221845/149‑1 (MQ=255)
     aTTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                                >  4:443141/1‑149 (MQ=255)
     aTTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                                >  1:769602/1‑149 (MQ=255)
      ttATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGa                                                                                                                                               <  4:565543/149‑1 (MQ=255)
        attttGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAgg                                                                                                                                             <  3:180605/146‑1 (MQ=255)
        atatTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAgg                                                                                                                                             <  1:347849/149‑1 (MQ=255)
           ttGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCGTCATGAGGGct                                                                                                                                          <  2:210452/149‑2 (MQ=255)
            tGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAGCCCTCGTCATGAGGGcttt                                                                                                                                        <  1:642486/149‑4 (MQ=255)
              ccGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTaa                                                                                                                                       <  3:223867/149‑1 (MQ=255)
               cgcgACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTCGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGCCGAGGGCGTAAc                                                                                                                                      <  4:641235/149‑1 (MQ=255)
                  gacgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAtt                                                                                                                                   >  4:83036/1‑149 (MQ=255)
                  gacgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAtt                                                                                                                                   >  4:416216/1‑149 (MQ=255)
                     gaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACt                                                                                                                                >  2:588153/1‑149 (MQ=255)
                     gaAGCAACAGCAATGCTTACGCATAACCATAGCAAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATTTGAGATAAAGCCCTCATGACGAGGGCGTAACATTACt                                                                                                                                >  1:721270/1‑149 (MQ=255)
                       aGCAACAGCAATGCTTCCGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTAGTTGTGAACATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTAATGACGAGGGCGTAACATTACTCa                                                                                                                              <  4:452217/149‑1 (MQ=255)
                       aGCAACAGCAATGCTTACGCATCACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATAAGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCa                                                                                                                              <  4:624792/149‑1 (MQ=255)
                       aGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCa                                                                                                                              >  2:64462/1‑149 (MQ=255)
                        gCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAg                                                                                                                             <  2:706439/149‑1 (MQ=255)
                        gCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAg                                                                                                                             >  4:98675/1‑149 (MQ=255)
                        gCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCAAGACGAGGGCGTAACATTACTCAg                                                                                                                             >  3:488891/1‑149 (MQ=255)
                           aCAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCaa                                                                                                                          >  4:540817/1‑149 (MQ=255)
                             aGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATa                                                                                                                        >  1:272329/1‑149 (MQ=255)
                               cAATGCTTACGCATAACCATAGCGAAACTAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATaaa                                                                                                                      <  3:522098/149‑1 (MQ=255)
                                 aTGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACt                                                                                                                    <  3:83036/149‑1 (MQ=255)
                                 aTGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACt                                                                                                                    >  1:271366/1‑149 (MQ=255)
                                  tGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTg                                                                                                                   >  2:763921/1‑149 (MQ=255)
                                   gCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGa                                                                                                                  <  2:722227/149‑1 (MQ=255)
                                     ttACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGata                                                                                                                >  4:557378/1‑149 (MQ=255)
                                        cGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTc                                                                                                             <  4:515614/149‑1 (MQ=255)
                                           aTAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGt                                                                                                          >  1:272332/1‑149 (MQ=255)
                                             aaCCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCa                                                                                                        <  4:583299/149‑1 (MQ=255)
                                                     cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                                <  4:432104/149‑1 (MQ=255)
                                                     cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                                <  1:303552/149‑1 (MQ=255)
                                                     cGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTCGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                                <  4:320582/149‑1 (MQ=255)
                                                          aTAGTGGCGCAGTGTAAGGTTGTTGCTAATATTGAGTTGCAAATATGTCGGTTTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGATGGCGTAACCTTTCTCAGCACTAATCTGATATTCCGTCAGGCTGGAATActc                                                                                           <  2:588440/149‑1 (MQ=255)
                                                           tAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATActct                                                                                          <  1:39302/149‑1 (MQ=255)
                                                                gcgcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGcc                                                                                     <  4:736772/149‑1 (MQ=255)
                                                                  gcAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAg                                                                                   <  1:30705/149‑1 (MQ=255)
                                                                    aGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGa                                                                                 <  3:675917/149‑1 (MQ=255)
                                                                    aGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGa                                                                                 <  3:540817/149‑1 (MQ=255)
                                                                       gtAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGc                                                                              >  1:406189/1‑149 (MQ=255)
                                                                         aaGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAg                                                                            <  1:674367/149‑1 (MQ=255)
                                                                         aaGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAg                                                                            >  1:294447/1‑149 (MQ=255)
                                                                                ttgtGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAg                                                                     <  2:189851/149‑1 (MQ=255)
                                                                                 tgtgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGt                                                                    <  4:697662/149‑1 (MQ=255)
                                                                                  gtgAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTc                                                                   <  4:667365/149‑1 (MQ=255)
                                                                                       tatTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCGTCATGAGGGc                                                                                                                                           >  4:383850/1‑72 (MQ=255)
                                                                                       tatTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCGTCATGAGGGc                                                                                                                                           <  3:383850/72‑1 (MQ=255)
                                                                                       tatTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGtt                                                              >  2:341872/1‑149 (MQ=255)
                                                                                             gTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGcc                                                        >  4:374305/1‑149 (MQ=255)
                                                                                                  aaaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAGTAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCa                                                   >  3:516751/1‑149 (MQ=255)
                                                                                                  aaaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                                <  1:33249/104‑1 (MQ=255)
                                                                                                  aaaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGaa                                                                                                >  2:33249/1‑104 (MQ=255)
                                                                                                   aaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTccg                                                  >  1:773986/1‑147 (MQ=255)
                                                                                                   aaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAg                                                  <  2:752987/149‑1 (MQ=255)
                                                                                                   aaTATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAg                                                  <  2:677380/149‑1 (MQ=255)
                                                                                                    atatGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAgg                                                 <  3:691976/149‑1 (MQ=255)
                                                                                                        gTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtgg                                             <  2:148392/149‑1 (MQ=255)
                                                                                                         tCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtggt                                            >  1:499774/1‑149 (MQ=255)
                                                                                                         tCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtggt                                            >  3:586772/1‑149 (MQ=255)
                                                                                                           ggTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTc                                          <  3:13502/149‑1 (MQ=255)
                                                                                                            gtgtTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCg                                         >  3:108528/1‑149 (MQ=255)
                                                                                                               tttGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGc                                      >  3:148858/1‑149 (MQ=255)
                                                                                                                    tGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTcc                                 >  3:739758/1‑149 (MQ=255)
                                                                                                                      gTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCgg                               <  3:725663/149‑1 (MQ=255)
                                                                                                                         aTTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAg                            >  1:261708/1‑149 (MQ=255)
                                                                                                                          tttGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGa                           >  3:230574/1‑149 (MQ=255)
                                                                                                                             gAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTTCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAAc                        <  4:323454/149‑1 (MQ=255)
                                                                                                                             gAACAATATGAGATAAAGCCATCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAAc                        <  2:483266/149‑1 (MQ=255)
                                                                                                                               aCAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTc                      <  1:635783/149‑1 (MQ=255)
                                                                                                                                     tGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTc                >  3:386354/1‑149 (MQ=255)
                                                                                                                                         aTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCaga            <  2:400078/149‑1 (MQ=255)
                                                                                                                                         aTAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCaga            <  2:432039/149‑1 (MQ=255)
                                                                                                                                          tAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGtt                                                              <  3:412917/98‑1 (MQ=255)
                                                                                                                                          tAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGtt                                                              >  4:412917/1‑98 (MQ=255)
                                                                                                                                           aaaGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagc          <  2:769602/149‑1 (MQ=255)
                                                                                                                                            acgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtg                                              >  3:189847/3‑112 (MQ=255)
                                                                                                                                            acgCCCTCGTCATGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGggtg                                              <  4:189847/110‑1 (MQ=255)
                                                                                                                                              gCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGagccag       <  4:108528/149‑1 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                              <  2:267610/58‑1 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATa                                                                                              >  1:267610/1‑58 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt   >  3:90902/1‑149 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt   >  2:643675/1‑149 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt   >  1:165129/1‑149 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt   >  2:343779/1‑149 (MQ=255)
                                                                                                                                                  tCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCACGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCt   >  1:606344/1‑149 (MQ=255)
                                                                                                                                                   cATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCtt  <  1:588153/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTT  >  NC_000913/787632‑787927

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