Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,130,410 T→C 5.5% S157G (AGC→GGC)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,130,4100TC 5.5% 225.5 / 2.9 74S157G (AGC→GGC) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base C (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGG  >  NC_000913/2130266‑2130557
                                                                                                                                                |                                                                                                                                                   
tCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTg                                                                                                                                                 <  1:348895/149‑1 (MQ=255)
 cATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAATTGATGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCATCACATAGCTGTTGc                                                                                                                                                >  1:76858/1‑149 (MQ=255)
  aTAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCt                                                                                                                                               >  1:192059/1‑149 (MQ=255)
          tGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTCCTGTCTTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGATGATCTAAAt                                                                                                                                         >  4:19837/1‑147 (MQ=255)
            ttttAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCa                                                                                                                                       <  8:102100/147‑1 (MQ=255)
              ttAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGata                                                                                                                                   >  4:346212/1‑149 (MQ=255)
                       ttGCTACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaaa                                                                                                                           <  5:110939/148‑1 (MQ=255)
                       ttGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaaa                                                                                                                           <  2:33295/148‑1 (MQ=255)
                        tGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAAt                                                                                                                          <  7:35519/148‑1 (MQ=255)
                          aTACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTGAGCACATAGCTGt                                                                                                                                                   >  6:163857/1‑121 (MQ=255)
                          aTACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTGAGCACATAGCTGt                                                                                                                                                   <  5:163857/121‑1 (MQ=255)
                            aCCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGAGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACt                                                                                                                     >  1:168949/1‑149 (MQ=255)
                            aCCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCAAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACt                                                                                                                     >  8:172916/1‑149 (MQ=255)
                              cAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTc                                                                                                                    <  1:171871/148‑1 (MQ=255)
                              cAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTc                                                                                                                    >  2:237082/1‑148 (MQ=255)
                              cAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTc                                                                                                                    >  2:171871/1‑148 (MQ=255)
                              cAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTc                                                                                                                    >  5:146896/1‑148 (MQ=255)
                               aaaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGGCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                                  <  4:323204/149‑1 (MQ=255)
                               aaaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaaa                                                                                                                           >  4:124483/1‑140 (MQ=255)
                               aaaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGaaa                                                                                                                           <  3:124483/140‑1 (MQ=255)
                                aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAGTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                                  >  3:369360/1‑148 (MQ=255)
                                aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTc                                                                                                                               >  3:265148/1‑135 (MQ=255)
                                aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTc                                                                                                                               <  4:265148/135‑1 (MQ=255)
                                aaCTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGACGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCat                                                                                                                 <  1:302265/149‑1 (MQ=255)
                                    gTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGCGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                                  <  4:285281/144‑1 (MQ=255)
                                    gTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGAGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCa                                                                                                                  >  3:285281/1‑144 (MQ=255)
                                       tGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAGTAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATAct                                                                                                           <  4:369360/148‑1 (MQ=255)
                                       tGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATAct                                                                                                           <  5:140498/148‑1 (MQ=255)
                                         aTGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATActgct                                                                                                        <  8:16625/149‑1 (MQ=255)
                                         aTGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGAAAATCAGATATTTCGAAATCGACTCCATAATActgct                                                                                                        <  1:237082/149‑1 (MQ=255)
                                                    aaTTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGa                                                                                             <  2:168949/149‑1 (MQ=255)
                                                      ttGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATAGACTCCATAATACTGCTGCCACCGGTGAg                                                                                            >  4:10464/1‑148 (MQ=255)
                                                         agaTTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGAt                                                                                         >  3:251693/1‑148 (MQ=255)
                                                          gaTTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATaa                                                                                       >  5:111094/1‑149 (MQ=255)
                                                               aTGCCAGAGAGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACTTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTtgt                                                                                  <  3:19837/149‑1 (MQ=255)
                                                                tGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTtgtg                                                                                 >  6:218118/1‑149 (MQ=255)
                                                                  ccAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGc                                                                                <  6:111094/148‑1 (MQ=255)
                                                                    agagCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCCGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAAttt                                                                                    <  3:242400/142‑1 (MQ=255)
                                                                    agagCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCCGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAAttt                                                                                    >  4:242400/1‑142 (MQ=255)
                                                                          aaaaaTGGCGGTTCGAAATAAAGTGCTGACGTGCGCGTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTAGAAATCGACTCCATAATACTGCTGCCACCGGTGATGATAATTTGTGCCATCTCCAc                                                                       >  2:267707/1‑149 (MQ=255)
                                                                            aaaTGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCAccg                                                                     >  5:272748/1‑149 (MQ=255)
                                                                               tGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTTCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCAccgccg                                                                  <  2:76858/149‑1 (MQ=255)
                                                                                    gTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgc                                                             <  1:182564/149‑1 (MQ=255)
                                                                                    gTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgc                                                             >  1:56880/1‑149 (MQ=255)
                                                                                    gTTCGAAATAAAGTGCTGTCGTGCGCATCCCGCCGAATTTAATGAAATTCAGCACATAGCAGTTGCTGAAAATCAGATATTTCGAAATCGACTCCATAATACAGCTGCCACCGGAGAGGATAATTTGTGCCATCTCCACCGCCGCCAgc                                                             >  6:228127/1‑149 (MQ=255)
                                                                                                 tGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGcc                                                <  1:301844/149‑1 (MQ=255)
                                                                                                 tGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGcc                                                <  1:226097/149‑1 (MQ=255)
                                                                                                   cTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGa                                                                                             >  8:209719/1‑102 (MQ=255)
                                                                                                   cTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGa                                                                                             <  7:209719/102‑1 (MQ=255)
                                                                                                        gTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                         <  3:346212/149‑1 (MQ=255)
                                                                                                        gTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCAGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                         >  1:359234/1‑149 (MQ=255)
                                                                                                           cgcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGAccacca                                       <  7:172916/148‑1 (MQ=255)
                                                                                                           cgcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGAccacca                                       <  7:126174/148‑1 (MQ=255)
                                                                                                           cgcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATAGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGAccacca                                       <  3:10464/148‑1 (MQ=255)
                                                                                                            gcTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGaccac                                         >  6:150695/1‑145 (MQ=255)
                                                                                                              ttGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                    >  2:106539/1‑148 (MQ=255)
                                                                                                                 ccgccgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATaaaa                                <  8:213854/149‑1 (MQ=255)
                                                                                                                  cgccgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAg                               <  2:56880/149‑1 (MQ=255)
                                                                                                                     cgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                    >  1:237422/1‑141 (MQ=255)
                                                                                                                     cgAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGAt                                    <  2:237422/141‑1 (MQ=255)
                                                                                                                           tAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGaa                             >  7:305256/1‑142 (MQ=255)
                                                                                                                           tAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGaa                             <  8:305256/142‑1 (MQ=255)
                                                                                                                              tGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATc                   >  3:248602/1‑149 (MQ=255)
                                                                                                                               gAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCt                  >  6:324937/1‑149 (MQ=255)
                                                                                                                                   ttCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgc                                                             <  1:312363/102‑1 (MQ=255)
                                                                                                                                   ttCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAgc                                                             >  2:312363/1‑102 (MQ=255)
                                                                                                                                       gCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATgg          <  5:218118/149‑1 (MQ=255)
                                                                                                                                       gCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATgg          <  2:192059/149‑1 (MQ=255)
                                                                                                                                        cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCCCCACCAGCCCGACCCCCAGATAAAAGAAACGTAATATCTTGCGATGGt         >  2:320804/1‑149 (MQ=255)
                                                                                                                                        cacaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGt         >  4:149634/1‑149 (MQ=255)
                                                                                                                                         acaTAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGtt        <  5:324937/149‑1 (MQ=255)
                                                                                                                                           aTAGCCGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGaaa                            <  7:82253/127‑1 (MQ=255)
                                                                                                                                           aTAGCCGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGaaa                            >  8:82253/1‑127 (MQ=255)
                                                                                                                                               cTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCAGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTgg  >  1:363221/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TCATAGCATCTGTTTTAGGCGTTTTGATACCAAACTGTTTGATGCTGAGCCAAATTGAGATTAATGCCAGAGCGAAAAATGCCGGTTCGAAATAAAGTGCTGTCGTGCGCTTGCCGCCGAATTTAATGAAATTCAGCACATAGCTGTTGCTGTAAATCAGATATTTCGAAATCGACTCCATAATACTGCTGCCACCGGTGAGGATAATTTGTGCCATCTCCACCGCCGCCAGCGCCACCACCAGCCCGACCACCAGATAAAAGAAACGTAATATCTTGCGATGGTTGTGTGG  >  NC_000913/2130266‑2130557

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