Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,130,405 A→T 5.5% Y158* (TAT→TAA wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,130,4050AT 5.5% 242.9 / 3.7 72Y158* (TAT→TAAwcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base T (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTT  >  NC_000913/2130263‑2130551
                                                                                                                                              |                                                                                                                                                  
gaATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCt                                                                                                                                               >  5:115621/1‑148 (MQ=255)
 aaTCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTg                                                                                                                                              <  7:88433/148‑1 (MQ=255)
  aTCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGt                                                                                                                                             >  3:18166/1‑148 (MQ=255)
  aTCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCCGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGt                                                                                                                                             >  7:1300/1‑148 (MQ=255)
    cATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTg                                                                                                                                           >  1:100029/1‑148 (MQ=255)
           tCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTaaa                                                                                                                                    >  3:125017/1‑148 (MQ=255)
               ttttAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATc                                                                                                                                  <  6:115621/146‑1 (MQ=255)
               ttttAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAg                                                                                                                                <  1:102860/148‑1 (MQ=255)
                tttAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAg                                                                                                                                <  8:272041/147‑1 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACTTAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  1:171271/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACTTAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  1:166227/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAg                                                                                                                                >  4:194285/1‑138 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAg                                                                                                                                <  3:194285/138‑1 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  6:125383/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  4:122665/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  5:207423/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGa                                                                                                                       >  4:184632/1‑147 (MQ=255)
                         tttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGAATAATGCCAGAGCGAAAAATCCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTTCCGCCGAATATAATGAAATTCAGCACAAATCTGTTGCAGTAAATTAGATATATCGa                                                                                                                       >  1:208004/1‑147 (MQ=255)
                             aTACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACTTAGCTGTTGCTGTAAATCAGATATTTCGAAATCg                                                                                                                  <  2:171271/148‑1 (MQ=255)
                             aTACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCg                                                                                                                  <  2:16331/148‑1 (MQ=255)
                              tACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACTTAGCTGTTGCTGTAAATCAGATATTTCGAAATCg                                                                                                                  <  2:166227/147‑1 (MQ=255)
                              tACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTaa                                                                                                                                     <  2:174715/128‑1 (MQ=255)
                              tACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTaa                                                                                                                                     >  1:174715/1‑128 (MQ=255)
                                ccAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACt                                                                                                               >  7:220996/1‑148 (MQ=255)
                                ccAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACt                                                                                                               >  1:285974/1‑148 (MQ=255)
                                  aaaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCTAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTcc                                                                                                             <  8:220996/148‑1 (MQ=255)
                                   aaCTGTTTGATGCTTAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                            <  1:294483/148‑1 (MQ=255)
                                   aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                            <  6:191261/148‑1 (MQ=255)
                                    aCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCat                                                                                                           <  5:125383/148‑1 (MQ=255)
                                        tttGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCataat                                                                                                        <  3:298789/147‑1 (MQ=255)
                                        tttGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCataat                                                                                                        <  4:125017/147‑1 (MQ=255)
                                           gATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATCct                                                                                                     <  5:14917/147‑1 (MQ=255)
                                            aTGCTGAGCCTAATTGAGATTAATGGCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGGCGTATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCGATAATActgc                                                                                                   <  2:57549/148‑1 (MQ=255)
                                                     cAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTAGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCgg                                                                                          <  2:251362/148‑1 (MQ=255)
                                                        aTTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCataat                                                                                                        >  4:229135/1‑131 (MQ=255)
                                                        aTTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCataat                                                                                                        <  3:229135/131‑1 (MQ=255)
                                                         ttGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTcc                                                                                                             <  3:189878/125‑1 (MQ=255)
                                                         ttGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTcc                                                                                                             >  4:189878/1‑125 (MQ=255)
                                                          tGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAg                                                                                      >  1:281520/1‑147 (MQ=255)
                                                                  aTGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTTTTGCTCTAAATCAGATATTTCGAAATCGACTCAATAATACAGCTGCCACCGGTGAGGATAATAtg                                                                             >  8:69340/1‑148 (MQ=255)
                                                                  aTGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTtg                                                                             <  2:281520/148‑1 (MQ=255)
                                                                   tGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTtgt                                                                            >  2:160095/1‑148 (MQ=255)
                                                                    gCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTtgtg                                                                           >  7:214810/1‑148 (MQ=255)
                                                                      cAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGcc                                                                         <  6:207423/148‑1 (MQ=255)
                                                                      cAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGcc                                                                         <  4:18166/148‑1 (MQ=255)
                                                                        gatcGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGAGAGGATAATTTGTGCCAt                                                                       >  4:56933/4‑148 (MQ=255)
                                                                        gagCGAAAAATTCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCAt                                                                       >  1:72272/1‑148 (MQ=255)
                                                                        gagCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTg                                                                                        <  8:194034/131‑1 (MQ=255)
                                                                        gagCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTg                                                                                        >  7:194034/1‑131 (MQ=255)
                                                                          gCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCAtct                                                                     >  8:314192/1‑148 (MQ=255)
                                                                          gCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCAtct                                                                     <  3:122665/148‑1 (MQ=255)
                                                                              aaaaTGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCAAc                                                                 >  3:168103/1‑148 (MQ=255)
                                                                                         tCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgcg                                                      >  7:198960/1‑148 (MQ=255)
                                                                                           gAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgcgc                                                     <  2:72272/147‑1 (MQ=255)
                                                                                            aaaTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaa                                                                                                                      >  7:304305/1‑81 (MQ=255)
                                                                                            aaaTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaa                                                                                                                      <  8:304305/81‑1 (MQ=255)
                                                                                            aaaTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGcca                                                   <  3:184632/148‑1 (MQ=255)
                                                                                                 aaGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCTAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGccaccacc                                              <  2:100029/148‑1 (MQ=255)
                                                                                                    tGCTGTCGTGCGCTTGCCGCCGATTTTAATTAAATTCAGCACATAGCGGTTGCTGTAATTGATATATTTATAAATGGACTCCATAATACTGCTTCCAAATTTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCAGCACCAGc                                           <  4:168103/148‑1 (MQ=255)
                                                                                                            tGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                   <  7:314192/148‑1 (MQ=255)
                                                                                                                      gccgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGCCCTCCAGATaa                            >  7:277781/1‑145 (MQ=255)
                                                                                                                        cgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGGGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGa                        >  4:47334/1‑147 (MQ=255)
                                                                                                                         gAATTTAATGAAATTCAGCACATAGCTTTTGCTGTAAATCAGATATTTCGAAATCGACTCCATACTACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTGCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGaaa                      <  3:177615/148‑1 (MQ=255)
                                                                                                                         gAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCa                                                                        >  3:84147/1‑98 (MQ=255)
                                                                                                                         gAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCa                                                                        <  4:84147/98‑1 (MQ=255)
                                                                                                                               aaTGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTaa                 <  7:69340/147‑1 (MQ=255)
                                                                                                                                  gAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATc             >  4:62584/1‑148 (MQ=255)
                                                                                                                                  gAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATc             >  1:188839/1‑148 (MQ=255)
                                                                                                                                  gAAATTCAGCAAATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGATGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATc             <  2:285974/148‑1 (MQ=255)
                                                                                                                                   aaaTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGGGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCt            <  3:56933/148‑1 (MQ=255)
                                                                                                                                     aTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCtt           >  7:26536/1‑147 (MQ=255)
                                                                                                                                     aTTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCtt           <  8:26536/147‑1 (MQ=255)
                                                                                                                                           cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATgg    <  1:160095/148‑1 (MQ=255)
                                                                                                                                           cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATgg    <  8:1300/148‑1 (MQ=255)
                                                                                                                                           cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATgg    >  1:130459/1‑148 (MQ=255)
                                                                                                                                            acaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGt   >  1:14267/1‑148 (MQ=255)
                                                                                                                                            acaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGt   >  1:9160/1‑148 (MQ=255)
                                                                                                                                             caTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGtt  <  3:62584/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GAATCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTT  >  NC_000913/2130263‑2130551

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