Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 787,784 C→T 7.0% intergenic (‑189/+13) 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,7840CT 7.0% 275.6 / 7.8 86intergenic (‑189/+13)gpmA/galMphosphoglyceromutase 1/aldose 1‑epimerase; type‑1 mutarotase; galactose mutarotase
Reads supporting (aligned to +/- strand):  ref base C (37/43);  new base T (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.46e-01

ATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTT  >  NC_000913/787637‑787927
                                                                                                                                                   |                                                                                                                                               
aTTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                                >  1:769602/1‑149 (MQ=255)
aTTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACg                                                                                                                                                >  4:443141/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:416216/1‑149 (MQ=255)
             gacgaAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACAtt                                                                                                                                   >  4:83036/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:540817/149‑1 (MQ=255)
                                                               aGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGa                                                                                 <  3:675917/149‑1 (MQ=255)
                                                                  gtAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGc                                                                              >  1:406189/1‑149 (MQ=255)
                                                                    aaGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAg                                                                            >  1:294447/1‑149 (MQ=255)
                                                                    aaGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAg                                                                            <  1:674367/149‑1 (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)
                                                                                                                                                   |                                                                                                                                               
ATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCGGTGTTTGCTGGTGATTTGAACAATATGAGATAAAGCCCTCATGACGAGGGCGTAACATTACTCAGCAATAAACTGATATTCCGTCAGGCTGGAATACTCTTCGCCAGGACGCAGGAAGCAGTCCGGTTGCGGCCATTCAGGGTGGTTCGGGCTGTCCGGTAGAAACTCGCTTTCCAGAGCCAGCCCTT  >  NC_000913/787637‑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: