Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,276,886 C→A 7.3% 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 7.3% 476.1 / 21.1 137intergenic (+33/‑116)garD/prlFD‑galactarate dehydrogenase/antitoxin of the SohA(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (66/61);  new base A (7/3);  total (73/64)
Fisher's exact test for biased strand distribution p-value = 3.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 137 total aligned reads displayed.

GGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTG  >  NC_000913/3276741‑3277029
                                                                                                                                                  |                                                                                                                                               
 gAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAg                                                                                                                                             >  1:220003/1‑149 (MQ=255)
    gttgttCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGc                                                                                                                                          >  3:254007/1‑149 (MQ=255)
          ccACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGtt                                                                                                                                    <  3:847541/149‑1 (MQ=255)
           cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGc                                                                                                                                          >  1:1031289/1‑142 (MQ=255)
           cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGc                                                                                                                                          <  2:1031289/142‑1 (MQ=255)
           cACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGttt                                                                                                                                   >  2:628123/1‑149 (MQ=255)
            aCTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAATCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACGTGGC‑CTTAGTGCCATGtttt                                                                                                                                  <  1:376750/149‑1 (MQ=255)
             cTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGttttt                                                                                                                                 >  4:177021/1‑149 (MQ=255)
             cTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttt                                                                                                                                   >  3:802837/1‑147 (MQ=255)
               ttATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTAt                                                                                                                               <  1:849423/149‑1 (MQ=255)
                     tcGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTa                                                                                                                         <  1:926331/149‑1 (MQ=255)
                       gACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                       >  4:622187/1‑149 (MQ=255)
                        aCGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAg                                                                                                                      >  4:592709/1‑149 (MQ=255)
                         cGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGc                                                                                                                     <  1:628123/149‑1 (MQ=255)
                         cGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGc                                                                                                                     <  1:930926/149‑1 (MQ=255)
                             gCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccccc                                                                                                                 <  4:511205/149‑1 (MQ=255)
                             gCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGccccc                                                                                                                 >  2:243004/1‑149 (MQ=255)
                                     gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                       <  2:671080/135‑1 (MQ=255)
                                     gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTaaa                                                                                                                       >  1:671080/1‑135 (MQ=255)
                                     gaagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAt                                                                                                         <  2:894537/149‑1 (MQ=255)
                                      aagaagaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCAtt                                                                                                        >  4:272486/1‑149 (MQ=255)
                                          agaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTCCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAATCCCCCACGTCCATtaata                                                                                                    >  3:730404/1‑149 (MQ=255)
                                           gcaacccTTCTCGCATCAATGTGTGCTGCATGACCAGCGGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTTGAACGTCCCTCGCAAAACATGGA‑CTTAGTCCCATGTTTTTATTGTTTAAAGCCCCGACGTCCATtaataa                                                                                                   <  3:261976/144‑1 (MQ=255)
                                           gaaAACCTTCTCGGTTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                   >  3:81019/1‑149 (MQ=255)
                                           gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                   >  3:494823/1‑149 (MQ=255)
                                           gaaAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATtaataa                                                                                                   >  3:503003/1‑149 (MQ=255)
                                            aaaaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTaataat                                                                                                  >  3:621345/1‑149 (MQ=255)
                                              aaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCAc                                                                                                               <  3:1024579/134‑1 (MQ=255)
                                              aaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCAc                                                                                                               >  4:1024579/1‑134 (MQ=255)
                                              aaCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgc                                                                                                <  4:52990/149‑1 (MQ=255)
                                               aCCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgca                                                                                               >  3:541372/1‑149 (MQ=255)
                                                 cTTCTCGGATCAATGGGGGCTGCATGACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                             >  4:696688/1‑149 (MQ=255)
                                                 cTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatt                                                                                             <  4:689593/149‑1 (MQ=255)
                                                    ctcGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcatttgc                                                                                          <  1:1127185/149‑1 (MQ=255)
                                                        gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGttt                                                                                                                                   >  3:21173/1‑104 (MQ=255)
                                                        gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                      >  2:429391/1‑149 (MQ=255)
                                                        gATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTa                                                                                      >  4:1101604/1‑149 (MQ=255)
                                                         aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                     >  3:898680/1‑149 (MQ=255)
                                                         aTCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTAc                                                                                     <  3:177021/149‑1 (MQ=255)
                                                                ggggCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGtt                                                                              >  2:1040546/1‑149 (MQ=255)
                                                                        aTAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTta                                                                      >  4:453017/1‑149 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                            <  1:964945/125‑1 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                            >  2:964945/1‑125 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgca                                                                                               >  1:1140225/1‑122 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGt                                                                               >  4:178333/1‑138 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTtata                                                                    <  4:621345/149‑1 (MQ=255)
                                                                          aaCCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACACGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgca                                                                                               <  2:1140225/122‑1 (MQ=255)
                                                                            ccAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                  <  4:730404/149‑1 (MQ=255)
                                                                            ccAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                  <  3:619440/149‑1 (MQ=255)
                                                                                cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                              <  1:351808/149‑1 (MQ=255)
                                                                                cTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGCCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGa                                                              <  2:373811/149‑1 (MQ=255)
                                                                                   gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCAc                                                                                                               >  1:754981/1‑97 (MQ=255)
                                                                                   gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCAc                                                                                                               <  2:754981/97‑1 (MQ=255)
                                                                                   gCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                           <  4:254007/149‑1 (MQ=255)
                                                                                     ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttttgcga                                                                                                                              >  4:427812/1‑76 (MQ=25)
                                                                                     ggTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑ACTAAGGCCATGttttgcga                                                                                                                              <  3:427812/80‑5 (MQ=25)
                                                                                         tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTa                                                                                 <  3:592139/121‑1 (MQ=255)
                                                                                         tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                     >  4:562636/1‑149 (MQ=255)
                                                                                         tttAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAAc                                                     >  4:927627/1‑149 (MQ=255)
                                                                                            aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTtat                                                                     >  2:1049888/1‑130 (MQ=255)
                                                                                            aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                  <  3:453017/149‑1 (MQ=255)
                                                                                            aaCCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAgg                                                  <  4:584089/149‑1 (MQ=255)
                                                                                               ccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCa                                                           >  3:1127276/1‑137 (MQ=255)
                                                                                               ccGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                               <  2:1148570/149‑1 (MQ=255)
                                                                                                  gCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                  >  2:54044/1‑127 (MQ=255)
                                                                                                  gCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAt                                                                  <  1:54044/127‑1 (MQ=255)
                                                                                                    aCCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTaaaa                                          <  3:927627/149‑1 (MQ=255)
                                                                                                     ccGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                         <  1:1040546/149‑1 (MQ=255)
                                                                                                         tGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACa                                     <  3:1101604/149‑1 (MQ=255)
                                                                                                             cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                 >  2:132978/1‑149 (MQ=255)
                                                                                                             cTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGa                                 >  2:371713/1‑149 (MQ=255)
                                                                                                              tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                            <  1:186714/89‑1 (MQ=255)
                                                                                                              tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAtgcattt                                                                                            >  2:186714/1‑89 (MQ=255)
                                                                                                              tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                               >  2:208522/1‑134 (MQ=255)
                                                                                                              tGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTg                                               <  1:208522/134‑1 (MQ=255)
                                                                                                                aTTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATc                              >  4:398971/1‑149 (MQ=255)
                                                                                                                    cTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCt                                                                                   >  3:226281/1‑92 (MQ=255)
                                                                                                                    cTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCt                                                                                   <  4:226281/92‑1 (MQ=255)
                                                                                                                     ttttGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAAttt                                                                <  3:358997/110‑1 (MQ=255)
                                                                                                                           aaCGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTaa                   >  4:299439/1‑149 (MQ=255)
                                                                                                                            aCGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAAGGCATTTGCATTACCTTAGTTCAAGCTTATACTTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAAt                  <  3:634024/149‑1 (MQ=255)
                                                                                                                                cccTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctc              <  4:1093736/149‑1 (MQ=255)
                                                                                                                                cccTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctc              <  3:696688/149‑1 (MQ=255)
                                                                                                                                 ccTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcg             >  4:515084/1‑149 (MQ=255)
                                                                                                                                  cTCGCAAAATATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgc            >  3:1073413/1‑149 (MQ=255)
                                                                                                                                   tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct           >  1:900462/1‑149 (MQ=255)
                                                                                                                                   tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgct           >  4:1136031/1‑149 (MQ=255)
                                                                                                                                   tCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcg             >  2:880067/1‑147 (MQ=255)
                                                                                                                                    cGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATgctcgctc          <  4:898680/149‑1 (MQ=255)
                                                                                                                                     taaaaaCATGGCACTAAG‑GCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct         >  3:643001/3‑149 (MQ=255)
                                                                                                                                     taaaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct         >  1:832420/3‑149 (MQ=255)
                                                                                                                                     taaaaaCATGGC‑ACTAAGGACATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct         >  1:1013563/3‑149 (MQ=255)
                                                                                                                                     gCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGctct         <  2:905485/149‑1 (MQ=255)
                                                                                                                                       aaaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAg                                         <  1:88716/115‑1 (MQ=255)
                                                                                                                                       aaaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCa       >  2:802036/1‑149 (MQ=255)
                                                                                                                                        aaaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGACCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc      <  3:986753/149‑1 (MQ=255)
                                                                                                                                        aaaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCAc      >  2:649728/1‑149 (MQ=255)
                                                                                                                                         aaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg     >  1:1078140/1‑149 (MQ=255)
                                                                                                                                         aaCATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg     >  3:188043/1‑149 (MQ=255)
                                                                                                                                         aaCATGGC‑ACTAAGGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACg     >  4:223765/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGAAGTTGTTCCACTTTATTCTCGACGTCGCCAGCGGGAAGAAGAAAACCTTCTCGGATCAATGGGGGCTGCATAACCAGCTGGCGGTGTTTAACCCGGCACCGGTGACCTGATTTCTTTTGGAACGTCCCTCGCAAAACATGGC‑CTTAGTGCCATGTTTTTATTGTTTAAAGCCCCCACGTCCATTAATAATGCATTTGCATTACCTTAGTTCAAGCTTATAATTTGAGCAGAAAACAGGCTGTAAAAGGACAGTGAATCATGCCCGCTAATGCTCGCTCTCACGCTG  >  NC_000913/3276741‑3277029

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: ‑