Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,955,071 C→A 7.7% E311D (GAG→GAT recB ← exonuclease V (RecBCD complex), beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,955,0710CA 7.7% 297.4 / 11.7 91E311D (GAG→GATrecBexonuclease V (RecBCD complex), beta subunit
Reads supporting (aligned to +/- strand):  ref base C (40/44);  new base A (3/4);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.89e-01

TCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGC  >  NC_000913/2954928‑2955217
                                                                                                                                               |                                                                                                                                                  
tCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGGTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCa                                                                                                                                                 >  6:73920/1‑147 (MQ=255)
 cGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCaaa                                                                                                                                               >  5:131641/1‑148 (MQ=255)
  gAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCaaa                                                                                                                                               >  8:205242/1‑147 (MQ=255)
   aGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCaaa                                                                                                                                               >  4:210337/1‑146 (MQ=255)
   aGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACa                                                                                                                                             <  1:33254/148‑1 (MQ=255)
     ccGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGt                                                                                                                                           <  3:234855/148‑1 (MQ=255)
     ccGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGt                                                                                                                                           <  6:85800/148‑1 (MQ=255)
       gACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTgg                                                                                                                                         <  1:176905/148‑1 (MQ=255)
       gACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTgg                                                                                                                                         <  1:78387/148‑1 (MQ=255)
       gACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTgg                                                                                                                                         <  6:266614/148‑1 (MQ=255)
                tGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCgg                                                                                                                                <  8:23765/148‑1 (MQ=255)
                tGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCgg                                                                                                                                >  6:254002/1‑148 (MQ=255)
                 gTCATCAAAACCCAAATCGCCACGCCGGCGTATTTCACGCGCTACTGTATCGCGGATCACAGCAAAAGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCggg                                                                                                                               >  6:162256/1‑148 (MQ=255)
                    atcaAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGtt                                                                                                                            <  5:254002/148‑1 (MQ=255)
                        aaaCCCAATTCGCCACGGCGGCGTTTTTCACGCTCTACTGTTTCGCGGATCTCAGCCAATGAGCGGGTGATCACCAGATAGCGGATCGACAATGGTTCTGAAAGAAGTTTATCGATCGCCTAAAACAGTGGATGACGCGGGGTTTccc                                                                                                                        >  1:302876/1‑148 (MQ=255)
                          aCCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccccc                                                                                                                      >  7:158726/1‑148 (MQ=255)
                          aCCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTccccc                                                                                                                      >  3:276196/1‑148 (MQ=255)
                             cAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCg                                                                                                                   >  2:211140/1‑148 (MQ=255)
                               aTTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGc                                                                                                                 >  8:40184/1‑148 (MQ=255)
                                        ggcggcGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCg                                                                                                        <  3:2457/148‑1 (MQ=255)
                                               tttttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttc                                                                                                  >  3:220148/1‑147 (MQ=255)
                                                ttttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttct                                                                                                 >  2:79546/1‑147 (MQ=255)
                                                  ttCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtcttct                                                                                                 >  1:45411/1‑145 (MQ=255)
                                                    cACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGa                                                                                             <  8:158726/147‑1 (MQ=255)
                                                     aCGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGaaa                                                                                           >  1:319560/1‑148 (MQ=255)
                                                     aCGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGaa                                                                                            <  8:106974/147‑1 (MQ=255)
                                                     aCGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGaa                                                                                            <  2:319560/147‑1 (MQ=255)
                                                       gcgcTACTGTTTCGCGGATCTCAGCCAATTCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACg                                                                                         <  1:79546/148‑1 (MQ=255)
                                                       gcgcTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACg                                                                                         <  6:269799/148‑1 (MQ=255)
                                                         tcTAGTGTTTCGCTGATCTCAGCCAATGCGCGGGTGTTCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTTGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCt                                                                                       <  5:145340/147‑1 (MQ=255)
                                                         gcTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCt                                                                                       <  4:220148/148‑1 (MQ=255)
                                                           tACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCt                                                                                       <  1:211140/146‑1 (MQ=255)
                                                            aCTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAt                                                                                                      >  1:153753/1‑130 (MQ=255)
                                                            aCTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAt                                                                                                      <  2:153753/130‑1 (MQ=255)
                                                                 ttggcgGATCTCAGCCAATGCGCGGTTGATCACGAGATCGCGGATCTACAATGGCGCTTCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGGTTCCACCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGaa                                                                                <  4:42419/144‑1 (MQ=255)
                                                                    gcgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttt                                                                             <  2:45411/147‑1 (MQ=255)
                                                                    gcgGATCTCAGCCAATGCGCGGGTGATCACCAGATAGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttt                                                                             <  8:127067/147‑1 (MQ=255)
                                                                     cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt                                                                           >  3:297075/1‑148 (MQ=255)
                                                                     cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt                                                                           <  6:131641/148‑1 (MQ=255)
                                                                     cgGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAAttttt                                                                           >  2:201499/1‑148 (MQ=255)
                                                                             aGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGAc                                                                   >  4:177306/1‑148 (MQ=255)
                                                                             aGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTGTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCg                                                                     <  5:73920/146‑1 (MQ=255)
                                                                                 aaTGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCg                                                               >  5:161232/1‑148 (MQ=255)
                                                                                  aTGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCgg                                                              <  3:127380/148‑1 (MQ=255)
                                                                                   tGCGCGGGTGATCACCAGATCGCGGATCGCCAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGGGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCATCGACTCAgg                                                              >  1:300922/1‑147 (MQ=255)
                                                                                    gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                             >  7:117526/1‑147 (MQ=255)
                                                                                    gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                             >  7:271562/1‑147 (MQ=255)
                                                                                    gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                             >  4:298122/1‑147 (MQ=255)
                                                                                    gcgcGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                             >  3:127448/1‑147 (MQ=255)
                                                                                    gcgcGGGTGATCACCAGATAGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGGGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAACCGCTGGGAGAATTTTTCCAGCGACTCCGGc                                                             >  6:100279/1‑147 (MQ=255)
                                                                                     cgcgGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCaa                                                           >  6:154578/1‑148 (MQ=255)
                                                                                            gATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGAc                                                                   <  1:192895/133‑1 (MQ=255)
                                                                                             aTCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTATCCCCCCCGCCGTTCGTCGGATCTTCTAACAAACGCTGGGAGAATTTTTCAAGCGACt                                                                  >  2:192895/1‑133 (MQ=255)
                                                                                               cACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTg                                                 >  8:77406/1‑148 (MQ=255)
                                                                                               cACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTg                                                 >  7:275824/1‑148 (MQ=255)
                                                                                                aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTg                                                 <  7:40184/147‑1 (MQ=255)
                                                                                                aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGt                                                <  4:276196/148‑1 (MQ=255)
                                                                                                aCCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGt                                                <  4:127448/148‑1 (MQ=255)
                                                                                                 ccAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGtt                                               <  7:205242/148‑1 (MQ=255)
                                                                                                  cAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGtt                                               <  2:171464/147‑1 (MQ=255)
                                                                                                     aTCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtc                                           <  3:210337/148‑1 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATa                                                     <  5:217864/137‑1 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATa                                                     >  6:217864/1‑137 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAAc                                                   >  8:190809/1‑139 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAAc                                                   <  7:190809/139‑1 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                         >  2:297812/1‑148 (MQ=255)
                                                                                                      tCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCGGCCTTCGTGCCATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                         >  2:302351/1‑148 (MQ=255)
                                                                                                       cgctGATCGACAATGTGTCTGCAAGCAGTAGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCGCCCCGGCCTGCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGAATCCGGCAACTGATAACTGTTTGtct                                          <  8:93868/147‑1 (MQ=255)
                                                                                                       cgcgGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGtctc                                         >  7:302809/1‑148 (MQ=255)
                                                                                                               gACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCtcttcttct                                  <  3:298122/147‑1 (MQ=255)
                                                                                                                aCAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGc                                <  5:154578/148‑1 (MQ=255)
                                                                                                                    tGGTTCTGCAAGCAGTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAg                            <  5:302127/148‑1 (MQ=255)
                                                                                                                         cTGCAAGCAGGTGATCGCTCGCCGCAAACAGTGGCTGTCGCGGGGGTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGgcgc                        <  6:199338/147‑1 (MQ=255)
                                                                                                                           gCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGtt                                               <  2:238654/122‑1 (MQ=255)
                                                                                                                           gCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTCCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGtt                                               >  1:238654/1‑122 (MQ=255)
                                                                                                                            cAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGt                                                <  8:7608/120‑1 (MQ=255)
                                                                                                                            cAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGt                                                >  7:7608/1‑120 (MQ=255)
                                                                                                                             aaGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATc                   <  6:161232/148‑1 (MQ=255)
                                                                                                                                 aGTTGATCGATCGCCTCAAACAGTTGATGTCGCGGGGGTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGt               <  4:155046/148‑1 (MQ=255)
                                                                                                                                  gTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCGGCCTTCGTGCTATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCg             <  1:297812/148‑1 (MQ=255)
                                                                                                                                  gTTGATCGATCGCATCAAACAGTGGATGTCGCGGGGTTTCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCg             <  1:302351/148‑1 (MQ=255)
                                                                                                                                   ttGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTc              >  4:102149/1‑147 (MQ=255)
                                                                                                                                    tGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa            >  7:1668/1‑148 (MQ=255)
                                                                                                                                     gatcgatcGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGa            <  7:77406/147‑1 (MQ=255)
                                                                                                                                      atcgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtct                                                                                                    <  3:127356/58‑1 (MQ=255)
                                                                                                                                      atcgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGAtct                                                                                                    >  4:127356/1‑58 (MQ=255)
                                                                                                                                      atcgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTg                      <  5:73768/136‑1 (MQ=255)
                                                                                                                                      atcgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTg                      >  6:73768/1‑136 (MQ=255)
                                                                                                                                        cgatcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCa        >  1:320408/1‑148 (MQ=255)
                                                                                                                                          atcgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTTCCCAGGCGCTGATCTTGTCGATCCAtt      <  3:177306/148‑1 (MQ=255)
                                                                                                                                            cgCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCAtt      >  3:295657/1‑146 (MQ=255)
                                                                                                                                              ccTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGc  <  8:302809/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCGAGCCGACTTAACATGTCATCAAAACCCAATTCGCCACGGCGGCGTTTTTCACGCGCTACTGTTTCGCGGATCTCAGCCAATGCGCGGGTGATCACCAGATCGCGGATCGACAATGGTTCTGCAAGCAGTTGATCGATCGCCTCAAACAGTGGATGTCGCGGGGTTTCCCCCCCGGCCTTCGTGCGATCTTCTAAGAAACGCTGGGAGAATTTTTCCAGCGACTCCGGCAACTGATAACTGTTTGTCTCTTCTTCTGCCCAGGCGCTGATCTTGTCGATCCATTTAGC  >  NC_000913/2954928‑2955217

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