Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,886 C→A 6.2% intergenic (+33/‑116) garD → / → prlF D‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,276,8860CA 6.2% 736.3 / 27.0 211intergenic (+33/‑116)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (100/98);  new base A (7/6);  total (107/104)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 217 total aligned reads displayed.

CTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTG  >  NC_000913/3276739‑3277034
                                                                                                                                                    |                                                                                                                                                    
  ggAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTa                                                                                                                                                   >  3:905651/1‑149 (MQ=255)
  ggAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTa                                                                                                                                                   <  2:483554/149‑1 (MQ=255)
   gAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAg                                                                                                                                                  >  4:151300/1‑149 (MQ=255)
       ttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGcc                                                                                                                                              >  2:276528/1‑149 (MQ=255)
           tCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGt                                                                                                                                          >  3:1070584/1‑149 (MQ=255)
            ccACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATg                                                                                                                                           <  1:1268737/147‑1 (MQ=255)
            ccACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGtt                                                                                                                                         <  3:1257769/149‑1 (MQ=255)
              aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGCACTAAG‑GCCATGtttt                                                                                                                                       <  2:192451/149‑1 (MQ=255)
              aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGtttt                                                                                                                                       >  4:400726/1‑149 (MQ=255)
              aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtttt                                                                                                                                       <  3:427743/149‑1 (MQ=255)
              aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtttt                                                                                                                                       <  4:1057208/149‑1 (MQ=255)
                tttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTa                                                                                                                                     <  3:1225389/149‑1 (MQ=255)
                   aTTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTg                                                                                                                                  <  1:932085/149‑1 (MQ=255)
                    ttCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGt                                                                                                                                 >  4:1004694/1‑149 (MQ=255)
                      ctcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGttt                                                                                                                               >  4:507572/1‑149 (MQ=255)
                       tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtt                                                                                                                                         <  3:1298337/138‑1 (MQ=255)
                       tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAACTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTa                                                                                                                              >  1:1473782/1‑149 (MQ=255)
                         gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                            <  2:1220936/149‑1 (MQ=255)
                            gTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcc                                                                                                                         >  1:341725/1‑149 (MQ=255)
                             tCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccc                                                                                                                        >  3:636736/1‑149 (MQ=255)
                              cGCCAGCGGGAGGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTTATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGcccc                                                                                                                       >  3:1110607/1‑149 (MQ=255)
                                  aGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACg                                                                                                                   >  3:551401/1‑149 (MQ=255)
                                      ggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCa                                                                                                               <  4:1084264/149‑1 (MQ=255)
                                      ggAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCa                                                                                                               <  1:1723360/149‑1 (MQ=255)
                                       gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                              <  1:879500/149‑1 (MQ=255)
                                        aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTTGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                             <  3:785883/149‑1 (MQ=255)
                                          gaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                           <  1:119323/149‑1 (MQ=255)
                                          gaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                           >  2:262387/1‑149 (MQ=255)
                                          gaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaa                                                                                                           >  4:456057/1‑149 (MQ=255)
                                            agaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaata                                                                                                         >  2:1319372/1‑149 (MQ=255)
                                            agaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaata                                                                                                         <  1:276528/149‑1 (MQ=255)
                                             gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGTCCTGCTTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATTTTTTTATTGTTTAAAGccccc                                                                                                                      >  2:335208/1‑135 (MQ=255)
                                             gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                        >  3:191361/1‑149 (MQ=255)
                                             gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                        >  3:1661224/1‑149 (MQ=255)
                                             gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtttt                                                                                                                                       >  4:391805/1‑118 (MQ=255)
                                             gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtttt                                                                                                                                       <  3:391805/118‑1 (MQ=255)
                                              aaaaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                       <  2:1283993/149‑1 (MQ=255)
                                                  ccTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcat                                                                                                   <  2:1473782/149‑1 (MQ=255)
                                                   cTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                                  >  2:494566/1‑149 (MQ=255)
                                                          gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                           >  4:192411/1‑149 (MQ=255)
                                                           atcgacgggggCTGCATAACCAGCTGGAGGCGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                          <  4:1457872/143‑1 (MQ=255)
                                                           aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTCGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                          >  2:79029/1‑149 (MQ=255)
                                                            tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAcc                                                                                         <  3:1412160/149‑1 (MQ=255)
                                                            tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCAAATGCATTatt                                                                                         <  4:57608/149‑3 (MQ=255)
                                                            tCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttt                                                                                                                                        >  2:1160500/1‑102 (MQ=255)
                                                               aTGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                      <  4:1275021/149‑1 (MQ=255)
                                                                             aCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATaa                                                                        <  1:1172119/149‑1 (MQ=255)
                                                                                aGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                     <  1:262387/149‑1 (MQ=255)
                                                                                  cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                   <  2:946972/149‑1 (MQ=255)
                                                                                  cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                   <  4:848637/149‑1 (MQ=255)
                                                                                        gtgtTTCCCCCGGCACCGGTGTCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTATAGCCACCACATCCATTAATAATGCATTTGCATTACCTTAGTTCGAGCTTATAATTTGAGCAGaa                                                             >  4:585338/1‑149 (MQ=255)
                                                                                          gtTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGaaaa                                                           <  1:578178/149‑1 (MQ=255)
                                                                                           tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                          >  1:164495/1‑149 (MQ=255)
                                                                                             tAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                        <  2:919826/149‑1 (MQ=255)
                                                                                             tAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                        <  1:149120/149‑1 (MQ=255)
                                                                                             tAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGA‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAg                                                        <  4:1032406/149‑1 (MQ=255)
                                                                                              aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCtt                                                                                       <  1:685781/117‑1 (MQ=255)
                                                                                                cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                     <  2:1252413/149‑1 (MQ=255)
                                                                                                cccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTAAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCt                                                     >  2:1352781/1‑149 (MQ=255)
                                                                                                      aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGc‑c                                                                                                                                                      >  2:1077304/1‑46 (MQ=255)
                                                                                                      aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGc‑c                                                                                                                                                      <  1:1077304/46‑1 (MQ=255)
                                                                                                      aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                               >  1:986925/1‑149 (MQ=255)
                                                                                                      aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                               <  2:986925/149‑1 (MQ=255)
                                                                                                           tGACCTGATTTCTTTTGGAATGTCCCTCGCAAAACCTGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCACACGTCCATTAATAAGGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                          <  4:1670596/149‑1 (MQ=255)
                                                                                                           tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                          <  1:630628/149‑1 (MQ=255)
                                                                                                           tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                          <  3:151300/149‑1 (MQ=255)
                                                                                                             aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                   <  1:1607724/122‑1 (MQ=255)
                                                                                                             aCCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                                   >  2:1607724/1‑122 (MQ=255)
                                                                                                               cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTta                                                                           >  2:1484908/1‑112 (MQ=255)
                                                                                                                   tttCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                      <  4:1416636/97‑1 (MQ=255)
                                                                                                                     tCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATg                                <  3:40031/149‑1 (MQ=255)
                                                                                                                       ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                  >  2:1602002/1‑113 (MQ=255)
                                                                                                                       ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAg                                                                  <  1:1602002/113‑1 (MQ=255)
                                                                                                                       ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                     >  3:495931/1‑142 (MQ=255)
                                                                                                                       ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGcc                              >  2:835398/1‑149 (MQ=255)
                                                                                                                        tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGccc                             >  3:1580558/1‑149 (MQ=255)
                                                                                                                        tttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGccc                             >  3:259801/1‑149 (MQ=255)
                                                                                                                            gAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGtt                                                                                   <  3:1506872/91‑1 (MQ=255)
                                                                                                                            gAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTa                         >  2:859241/1‑149 (MQ=255)
                                                                                                                             aaCGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTaa                        >  3:1008843/1‑149 (MQ=255)
                                                                                                                              aCGTCCCTCGCAAAACATGGA‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATGTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                       <  1:79029/149‑1 (MQ=255)
                                                                                                                              aCGTCCCTCGCAAAACATGAC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                       <  4:826615/149‑1 (MQ=255)
                                                                                                                                gTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                       <  3:731484/147‑1 (MQ=255)
                                                                                                                                gTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                       >  4:731484/1‑147 (MQ=255)
                                                                                                                                 tCCCTCGGAAAACATGGC‑CGTACTGCCATGTTTTTATTGTTTAAAACCCCCACGTCCATTAATAATTCAATTGCATTACCTTAGTTGAAGCTTATAATTTGAGCGGATAACAGGCTGTAAAAGGACAGGGAATCTTGCCCGCTAAAgct                    >  2:1025388/1‑149 (MQ=255)
                                                                                                                                  cccTCGCATAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                     >  3:1556150/1‑131 (MQ=255)
                                                                                                                                  cccTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGaa                                     <  4:1556150/131‑1 (MQ=255)
                                                                                                                                    cTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAg                                                                                                                           <  2:679891/43‑1 (MQ=255)
                                                                                                                                    cTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc                 <  1:1012780/149‑1 (MQ=255)
                                                                                                                                    cTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc                 >  2:635535/1‑149 (MQ=255)
                                                                                                                                      cGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgctc               <  4:1324217/149‑1 (MQ=255)
                                                                                                                                       gCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct              <  1:494566/149‑1 (MQ=255)
                                                                                                                                        aaaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCtctc             >  4:1172194/2‑149 (MQ=255)
                                                                                                                                         aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                          >  2:1435684/1‑119 (MQ=255)
                                                                                                                                         aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAAt                                    >  1:1151330/1‑125 (MQ=255)
                                                                                                                                             cATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCt        >  4:25853/1‑149 (MQ=255)
                                                                                                                                              aTGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTg       >  1:8364/1‑149 (MQ=255)
                                                                                                                                                  c‑cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACt   <  1:1522633/149‑1 (MQ=255)
                                                                                                                                                    cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTg  >  3:296108/1‑149 (MQ=255)
                                                                                                                                                    cTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTg  >  2:589298/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTGGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTGTACTG  >  NC_000913/3276739‑3277034

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