Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,903 T→A 5.4% N326Y (AAT→TAT)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,9030TA 5.4% 299.9 / 2.6 87N326Y (AAT→TAT) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (40/42);  new base A (2/3);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGA  >  NC_000913/2129761‑2130050
                                                                                                                                              |                                                                                                                                                   
ggAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTa                                                                                                                                                <  1:153353/148‑1 (MQ=255)
      aaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                              >  1:13566/1‑144 (MQ=255)
      aaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                              <  2:13566/144‑1 (MQ=255)
      aaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                                                                                                         <  4:685685/149‑1 (MQ=255)
        aaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                                                                                                        <  2:365938/148‑1 (MQ=255)
         aGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                                                                                                      <  4:782895/149‑1 (MQ=255)
          gATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCg                                                                                                                                     <  3:508674/149‑1 (MQ=255)
           atatttttGCACGCTAAAATTGCGGTTATCACCAAAAGCGTTTATCAACATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCg                                                                                                                                     <  1:168445/142‑1 (MQ=255)
             ataATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                                                                                                   >  1:628865/1‑148 (MQ=255)
                 ttGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATcc                                                                                                                               <  1:175791/148‑1 (MQ=255)
                    cACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCa                                                                                                                            <  2:443014/148‑1 (MQ=255)
                    cACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCa                                                                                                                            <  2:296340/148‑1 (MQ=255)
                    cACGCGAAAATTGCGGTTATCACCAAAAGCGTAGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCa                                                                                                                            <  2:534724/148‑1 (MQ=255)
                     aCGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                                                                                          >  4:898047/1‑149 (MQ=255)
                              ttGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAatat                                                                                                                 >  3:734905/1‑149 (MQ=255)
                                 cGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                              <  3:579362/149‑1 (MQ=255)
                                   gTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGa                                                                                                            >  3:881261/1‑149 (MQ=255)
                                   gTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtat                                                                                                                                      >  1:453823/1‑119 (MQ=255)
                                   gTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtat                                                                                                                                      <  2:453823/123‑5 (MQ=255)
                                    ttATCACCAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                                                                                                          <  3:403207/149‑1 (MQ=255)
                                      aTCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                                                                                                   >  3:599536/1‑123 (MQ=255)
                                      aTCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                                                                                                   <  4:599536/123‑1 (MQ=255)
                                        cACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatga                                                                                                        >  1:501610/1‑148 (MQ=255)
                                        cACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatga                                                                                                        <  2:433446/148‑1 (MQ=255)
                                           cAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                                                                                                     >  1:599697/1‑148 (MQ=255)
                                            aaaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCg                                                                                                    >  1:501644/1‑148 (MQ=255)
                                             aaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGa                                                                                                   >  1:164446/1‑148 (MQ=255)
                                             aaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGa                                                                                                   >  1:48674/1‑148 (MQ=255)
                                             aaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                  <  4:734905/149‑1 (MQ=255)
                                                  gTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAATCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATAtt                                                                                             <  3:522978/149‑1 (MQ=255)
                                                   ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGCAAAATAAATAATCAGCAGATCCAACCCATAGACAATGGtttgtatc                                                                                                                                     >  3:690442/1‑103 (MQ=255)
                                                   ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGt                                                                                                                                                   >  1:451059/1‑94 (MQ=255)
                                                   ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGt                                                                                                                                                   <  2:451059/94‑1 (MQ=255)
                                                   ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATAtt                                                                                             <  1:654171/148‑1 (MQ=255)
                                                   ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtatc                                                                                                                                     <  4:690442/108‑6 (MQ=255)
                                                     gATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCAt                                                                                          >  4:67780/1‑149 (MQ=255)
                                                        catcatTTTTATCACTTTCCCCAAGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaa                                                                                        <  2:137611/148‑1 (MQ=255)
                                                         atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAg                                                                                      >  4:304149/1‑149 (MQ=255)
                                                           catTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt                                                                                     >  1:446400/1‑148 (MQ=255)
                                                             tttttATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa                                                                                   <  2:501644/148‑1 (MQ=255)
                                                                ttATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCa                                                                               >  3:802450/1‑149 (MQ=255)
                                                                  aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                             >  4:356616/1‑149 (MQ=255)
                                                                  aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaa                                                                              <  2:599697/148‑1 (MQ=255)
                                                                  aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaa                                                                              >  2:763412/1‑148 (MQ=255)
                                                                    cACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATCATCAGCAGATACAAACCATTg                                                                                                                                                    >  2:782798/1‑76 (MQ=255)
                                                                    cACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                            >  1:730935/1‑148 (MQ=255)
                                                                    cACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATc                                                                           >  3:533377/1‑149 (MQ=255)
                                                                    cACTTTCCCCATGTACCAAAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTg                                                                                                                                                    <  1:782798/76‑1 (MQ=255)
                                                                     aCTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCt                                                                          >  3:94681/1‑149 (MQ=255)
                                                                      cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTg                                                                         >  4:115186/1‑149 (MQ=255)
                                                                       tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTg                                                                         <  1:501991/148‑1 (MQ=255)
                                                                       tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTg                                                                         <  2:48674/148‑1 (MQ=255)
                                                                       tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                        <  4:94681/149‑1 (MQ=255)
                                                                       tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                        <  3:356616/149‑1 (MQ=255)
                                                                           cccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTac                                                                    <  4:533377/149‑1 (MQ=255)
                                                                            ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa                                                                                   <  4:380533/133‑1 (MQ=255)
                                                                            ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa                                                                                   >  3:380533/1‑133 (MQ=255)
                                                                             cATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTaca                                                                   <  2:628865/148‑1 (MQ=255)
                                                                                    cacaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCaa                                                            >  1:208228/1‑148 (MQ=255)
                                                                                     acaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAat                                                           >  2:893464/1‑148 (MQ=255)
                                                                                       aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAatat                                                         >  1:221912/1‑148 (MQ=255)
                                                                                        gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                        >  2:700513/1‑148 (MQ=255)
                                                                                        gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                        <  1:893464/148‑1 (MQ=255)
                                                                                        gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                        <  2:731368/148‑1 (MQ=255)
                                                                                          gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGt                                                     <  3:288167/149‑1 (MQ=255)
                                                                                           aTAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTc                                                    >  4:565753/1‑149 (MQ=255)
                                                                                            tAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                       <  4:605949/145‑1 (MQ=255)
                                                                                            tAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                       >  3:605949/1‑145 (MQ=255)
                                                                                              ataaCACCGCGAACCAGGAAAAAGAAATAATCAGCGGATACAAACCATAGTCTATGGTTTTTCCGACATGCGCACCTTTAAATATTCCGAAAAATGCGACATATTCATAAAGTGAGCCAAATGTGACTACTCTGTCAATATGGGTCaa                                                  <  2:255535/146‑1 (MQ=255)
                                                                                                   acCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                            >  4:179095/1‑149 (MQ=255)
                                                                                                    ccGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                            <  2:221912/148‑1 (MQ=255)
                                                                                                      gcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAAcc                                         >  4:559906/1‑149 (MQ=255)
                                                                                                      gcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTGAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAAc                                          >  2:380771/1‑148 (MQ=255)
                                                                                                         aaCCAGGAAAAATAAATAATCTACGGATACAAACCATTGTATAGGGTTTTTCCGACATCCGCACCGTTAAATATGCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                       <  1:25486/148‑1 (MQ=255)
                                                                                                          aCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGAcc                                     <  4:812091/149‑1 (MQ=255)
                                                                                                              ggAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCAtaa                                 >  4:875827/1‑147 (MQ=255)
                                                                                                               gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAc                                >  4:435374/1‑149 (MQ=255)
                                                                                                                      aaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCgg                          >  2:146106/1‑148 (MQ=255)
                                                                                                                             cagcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACa                   >  2:699164/1‑148 (MQ=255)
                                                                                                                               gcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAc                        >  3:646641/1‑141 (MQ=255)
                                                                                                                               gcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAc                        <  4:646641/141‑1 (MQ=255)
                                                                                                                                    tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGat            <  1:763412/148‑1 (MQ=255)
                                                                                                                                      cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataat         <  3:898047/149‑1 (MQ=255)
                                                                                                                                       aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataata        <  3:304149/149‑1 (MQ=255)
                                                                                                                                       aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataat         >  1:298810/1‑148 (MQ=255)
                                                                                                                                           cATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACACTACGATAATatgat    <  3:67780/149‑1 (MQ=255)
                                                                                                                                           cATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATatgat    >  3:852823/1‑149 (MQ=255)
                                                                                                                                             ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATAtgatga  >  4:775249/1‑149 (MQ=255)
                                                                                                                                             ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATAtgatga  >  3:364557/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGA  >  NC_000913/2129761‑2130050

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