Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,994,367 G→A 5.8% R114C (CGT→TGT)  uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,994,3670GA 5.8% 229.5 / 3.6 69R114C (CGT→TGT) uvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (29/36);  new base A (2/2);  total (31/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.92e-01

GCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTAC  >  NC_000913/1994221‑1994505
                                                                                                                                                  |                                                                                                                                          
gcgaTTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGAcgcg                                                                                                                                          >  4:33146/1‑149 (MQ=255)
    ttGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgggtg                                                                                                                                      <  2:181953/149‑1 (MQ=255)
     tGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgggtgg                                                                                                                                     <  3:579388/149‑1 (MQ=255)
       cGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgggtgggt                                                                                                                                   <  2:484200/149‑1 (MQ=255)
       cGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCgggtgggt                                                                                                                                   <  3:44572/149‑1 (MQ=255)
        gATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTa                                                                                                                                  >  3:72943/1‑149 (MQ=255)
           aaaCACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATca                                                                                                                               <  2:689822/149‑1 (MQ=255)
            aaCACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATcac                                                                                                                              >  1:224336/1‑149 (MQ=255)
             acacTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGACAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATcacc                                                                                                                             <  1:628988/149‑1 (MQ=255)
              cacTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATcacca                                                                                                                            >  2:491535/1‑149 (MQ=255)
                  aTTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCa                                                                                                                        <  2:211576/149‑1 (MQ=255)
                     ttCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGa                                                                                                                     <  4:381693/149‑1 (MQ=255)
                          aCTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGTTATCACCACTCAGGAAg                                                                                                                   >  3:568679/1‑146 (MQ=255)
                          aCTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAg                                                                                                                   <  4:568679/146‑1 (MQ=255)
                          aCTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATa                                                                                                                >  4:305527/1‑149 (MQ=255)
                            tGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATaaa                                                                                                              >  4:247779/1‑149 (MQ=255)
                             ggCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATaaaa                                                                                                             <  1:315145/149‑1 (MQ=255)
                              gCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAg                                                                                                            >  4:190167/1‑149 (MQ=255)
                                gAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGa                                                                                                          <  2:293758/149‑1 (MQ=255)
                                       gAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGTGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGAtt                                                                                                   <  2:418666/149‑1 (MQ=255)
                                       gAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGat                                                                                                         <  1:607887/143‑1 (MQ=255)
                                       gAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGat                                                                                                         >  2:607887/1‑143 (MQ=255)
                                        aaaaTCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACACGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGAttt                                                                                                  >  3:73492/1‑149 (MQ=255)
                                        aaaaTCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACACGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGAttt                                                                                                  >  1:674610/1‑149 (MQ=255)
                                             cTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTatcat                                                                                             <  3:647093/149‑1 (MQ=255)
                                              ttttGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTatcatc                                                                                            <  4:511846/149‑1 (MQ=255)
                                                   cAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTa                                                                                       <  3:327784/149‑1 (MQ=255)
                                                   cAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTa                                                                                       <  3:690509/149‑1 (MQ=255)
                                                     gTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCa                                                                                                                        <  3:540193/114‑1 (MQ=255)
                                                     gTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCa                                                                                                                        >  4:540193/1‑114 (MQ=255)
                                                     gTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGc                                                                                     >  2:630893/1‑149 (MQ=255)
                                                       aGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCaa                                                                                   >  2:29648/1‑149 (MQ=255)
                                                          gcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAAc                                                                                >  2:140566/1‑149 (MQ=255)
                                                          gcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACACGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAAc                                                                                <  4:73492/149‑1 (MQ=255)
                                                          gcCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACACGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAAc                                                                                <  2:674610/149‑1 (MQ=255)
                                                              gtgtTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTg                                                                            >  3:529572/1‑149 (MQ=255)
                                                                 tttCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGtt                                                                             >  3:599363/1‑145 (MQ=255)
                                                                 tttCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTaa                                                                         <  4:382794/149‑1 (MQ=255)
                                                                 tttCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTaa                                                                         >  3:148213/1‑149 (MQ=255)
                                                                 tttCACGTACGGCATAGCCATTCCGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGtt                                                                             <  4:599363/145‑1 (MQ=255)
                                                                      cGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCgg                                                                    <  3:247779/149‑1 (MQ=255)
                                                                          cGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGa                                                                <  4:72943/149‑1 (MQ=255)
                                                                           ggCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGAt                                                               <  3:455924/149‑1 (MQ=255)
                                                                              aTAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga                                                            <  4:529572/149‑1 (MQ=255)
                                                                                 gCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGt                                                                           >  4:623394/1‑131 (MQ=255)
                                                                                 gCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGt                                                                           <  3:623394/131‑1 (MQ=255)
                                                                                 gCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga                                                            >  2:75285/1‑146 (MQ=255)
                                                                                 gCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATaga                                                            <  1:75285/146‑1 (MQ=255)
                                                                                   cATTCGGGAACGGGCCGAAATATTCACCTTTGCCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTg                                                       <  1:487518/149‑1 (MQ=255)
                                                                                    aTTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGa                                                      >  4:137319/1‑149 (MQ=255)
                                                                                    aTTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGa                                                      <  3:305527/149‑1 (MQ=255)
                                                                                     ttCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGAt                                                     >  3:101520/1‑149 (MQ=255)
                                                                                        gggAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTa                                                  <  4:101520/149‑1 (MQ=255)
                                                                                        gggAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTa                                                  <  3:368153/149‑1 (MQ=255)
                                                                                           aaCGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGtt                                               >  4:283292/1‑149 (MQ=255)
                                                                                             cGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTtgt                                             <  3:282215/149‑1 (MQ=255)
                                                                                                  cGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTcc                                        >  1:286165/1‑149 (MQ=255)
                                                                                                   gAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCa                                       >  4:222528/1‑149 (MQ=255)
                                                                                                         ttCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTGCAGCAACa                                 >  2:596487/1‑149 (MQ=255)
                                                                                                           cACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAgc                               >  1:624578/1‑149 (MQ=255)
                                                                                                              cTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCt                            >  1:157016/1‑149 (MQ=255)
                                                                                                              cTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCt                            >  2:240466/1‑149 (MQ=255)
                                                                                                                       gCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCACTTCGGTctct                   <  1:596487/149‑1 (MQ=255)
                                                                                                                         ttCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCtgt                 <  1:491535/149‑1 (MQ=255)
                                                                                                                           cGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCtgtgt               <  3:283292/149‑1 (MQ=255)
                                                                                                                             caccacGATGCATCGCCAGACGCGGCTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGa             <  3:451511/149‑1 (MQ=255)
                                                                                                                                cacGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTa          >  4:642629/1‑149 (MQ=255)
                                                                                                                                        aTCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTAc  >  2:403937/1‑149 (MQ=255)
                                                                                                                                           gCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGttag                                                       >  3:390037/1‑91 (MQ=255)
                                                                                                                                           gCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTg                                                       <  4:390037/93‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
GCGATTGCGATAAACACTATTTTCGCACTGGCGAATGGGGAAAATCTTTTGCAGTAGCGCCAGTGTTTCACGTACGGCATAGCCATTCGGGAACGGGCCGAAATATTCACCTTTGGCATGCTTCGCACCACGATGCATCGCCAGACGCGGGTGGGTATCACCACTCAGGAAGATAAAAGGATATGATTTATCATCGCGTAGCAAAACGTTGTAACGCGGCTGATAGAGTTTGATGTAGTTGTGTTCCAGCAACAGCGCTTCGGTTTCTGTGTGAGTAACCGTTAC  >  NC_000913/1994221‑1994505

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