Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,998,273 G→A 5.5% L97L (TTG→TTA uvrD → DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,2730GA 5.5% 238.5 / 2.5 72L97L (TTG→TTAuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base A (2/2);  total (33/39)
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.59e-01

GCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACATC  >  NC_000913/3998125‑3998420
                                                                                                                                                    |                                                                                                                                                   
gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCCCGCGCCAGGGCGGTATGTGGGTCGGCCCCTTCCCCGGGCTGGCGCCGCGTTtg                                                                                                                                                     >  6:120457/1‑149 (MQ=255)
gCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                                                                                                                                                     >  1:290397/1‑149 (MQ=255)
  gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgct                                                                                                                                                   <  1:99151/149‑1 (MQ=255)
                gCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                                                                                                                     <  7:285473/149‑1 (MQ=255)
                gCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                                                                                                                     >  5:3337/1‑149 (MQ=255)
                gCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                                                                                                                     >  5:113862/1‑149 (MQ=255)
                gCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCata                                                                                                                                     <  4:307618/149‑1 (MQ=255)
                gCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCa                                                                                                                                       <  4:359985/147‑1 (MQ=255)
                  cATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATg                                                                                                                                   <  7:89868/149‑1 (MQ=255)
                    tACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGa                                                                                                                                 >  1:351898/1‑149 (MQ=255)
                      cTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCTTCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                                                                                                                               <  4:200425/149‑1 (MQ=255)
                      cTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCTTCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                                                                                                                               >  3:200425/1‑149 (MQ=255)
                      cTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACg                                                                                                                               >  4:272617/1‑149 (MQ=255)
                       tCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc                                                                                                                              >  4:246663/1‑149 (MQ=255)
                         gATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGc                                                                                                                              <  4:201669/147‑1 (MQ=255)
                              tGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCatat                                                                                                                                    <  7:204869/136‑1 (MQ=255)
                              tGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCatat                                                                                                                                    >  8:204869/1‑136 (MQ=255)
                                   gTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCa                                                                                                                  <  2:124403/149‑1 (MQ=255)
                                     gACGTTTACCAACAAAGCGGCGGCGGAGATGCTTCATCGTTTGGGGCATCTGATGGGAACTAGCCAGTGCGGTATGTGGGCCGGAACCTTCCACGGGCTTGCTCACCGTTTGCTGCGGGCGCACCATATGGACGCCTATCTGCCGCAgg                                                                                                                <  5:153365/149‑1 (MQ=255)
                                         tttACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAttt                                                                                                            >  2:99686/1‑149 (MQ=255)
                                           tACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc                                                                                                           <  8:280100/148‑1 (MQ=255)
                                             ccAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTc                                                                                                           >  5:352709/1‑146 (MQ=255)
                                              caacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                                                                                                                                               >  7:153468/1‑109 (MQ=255)
                                              caacaaAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                                                                                                                                               <  8:153468/109‑1 (MQ=255)
                                                  aaaGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCt                                                                                                   <  3:160289/149‑1 (MQ=255)
                                                    aGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGCCGCCAATCTGCCGCAGGATTTCCAGATCCTc                                                                                                  >  7:307869/1‑148 (MQ=255)
                                                           gcggAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTc                                                                                                  <  7:112406/141‑1 (MQ=255)
                                                           gcggAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTc                                                                                                  >  8:112406/1‑141 (MQ=255)
                                                           gcggAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGa                                                                                          <  1:326825/149‑1 (MQ=255)
                                                            cggAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGaa                                                                                         <  1:99686/149‑1 (MQ=255)
                                                            cggAGAGGCGTCATCGTATCGGGCACCTGAGGGGCACGCGCCAGGGCGGTATTTGGGGAGGCACCTTCCACGTGCTGGCGCACCGTTTGCTGCGGGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGaa                                                                                         <  4:26616/149‑1 (MQ=255)
                                                               agaTGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGAc                                                                                      <  6:3337/149‑1 (MQ=255)
                                                                 aTGCGTCATCGTATCGGGCAACTGATGGGCACGATCCAGGGCGGTATGTGGGTCGGCACCTTGCACGGGCTGTCGGACCGTTTGGTGCTTGCGCACCATATGGACGCCATGCTGCCGCAGGATTTCCAGATCCTCGTCAGTGAAGACCa                                                                                    <  6:248280/149‑1 (MQ=255)
                                                                   gCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgct                                                                                                                                                   >  5:130986/1‑84 (MQ=255)
                                                                   gCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgct                                                                                                                                                   <  6:130986/84‑1 (MQ=255)
                                                                    cGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGc                                                                                  >  5:45345/1‑148 (MQ=255)
                                                                     gTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTa                                                                                <  3:246663/149‑1 (MQ=255)
                                                                     gTCATCGTAGCGGGCAACTGATGGGCCCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTa                                                                                <  2:402300/149‑1 (MQ=255)
                                                                        aTCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGCTCCTCGACAGTGAAGACCAGCTACGc                                                                             <  3:26027/149‑1 (MQ=255)
                                                                        aTCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGc                                                                             >  2:218859/1‑149 (MQ=255)
                                                                           gTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTACTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc                                                                                                    <  6:278059/123‑1 (MQ=255)
                                                                           gTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTCGCGCACCGTTTACTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATcc                                                                                                    >  5:278059/1‑123 (MQ=255)
                                                                            tATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGc                                                                         >  6:338069/1‑149 (MQ=255)
                                                                                gggCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTaa                                                                     <  4:288650/149‑1 (MQ=255)
                                                                                 ggCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAg                                                                    >  8:49106/1‑149 (MQ=255)
                                                                                  gCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGc                                                                   <  1:218859/149‑1 (MQ=255)
                                                                                    aaCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGt                                                                 >  8:389634/1‑149 (MQ=255)
                                                                                       tGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa                                                                                              >  6:257532/1‑117 (MQ=255)
                                                                                       tGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACa                                                                                              <  5:257532/117‑1 (MQ=255)
                                                                                               aCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGc                                                      <  7:87103/149‑1 (MQ=255)
                                                                                               aCGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGc                                                      >  1:78623/1‑149 (MQ=255)
                                                                                                cGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGcc                                                     <  4:163763/149‑1 (MQ=255)
                                                                                                 gAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCa                                                    <  1:338416/149‑1 (MQ=255)
                                                                                                      aGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAAc                                               <  2:290397/149‑1 (MQ=255)
                                                                                                      aGGGCGGTATGTGGGTCGGCACCGTCCACGGGCTGGCGCCCCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGAGTTCCAGATCCTCGACAGTGAAGCCCAGCTACGCCTGCTGACGCGTCGGATCAAAGCCATGAAc                                               >  1:134790/1‑149 (MQ=255)
                                                                                                       gggCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAAc                                               >  6:51197/1‑148 (MQ=255)
                                                                                                            gTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTcgac                                         <  3:272617/149‑1 (MQ=255)
                                                                                                                  gggTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAg                                   >  6:110210/1‑149 (MQ=255)
                                                                                                                   ggTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGACCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAg                                   <  4:92014/148‑1 (MQ=255)
                                                                                                                    ggcGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGGGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCa                                 <  4:260131/147‑1 (MQ=255)
                                                                                                                    gTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAg                                   >  1:271443/1‑147 (MQ=255)
                                                                                                                     tCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAg                                <  2:223808/149‑1 (MQ=255)
                                                                                                                          aCCTTCCACGGGCTGGCGCACCGTTTACTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGAt                                                            >  6:222919/1‑116 (MQ=255)
                                                                                                                          aCCTTCCACGGGCTGGCGCACCGTTTACTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGAt                                                            <  5:222919/116‑1 (MQ=255)
                                                                                                                           ccTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGgccg                          >  5:266863/1‑149 (MQ=255)
                                                                                                                           ccTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGgccg                          >  3:190486/1‑149 (MQ=255)
                                                                                                                            cTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGgccgc                         >  3:378435/1‑149 (MQ=255)
                                                                                                                            cTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGgccgc                         >  7:388862/1‑149 (MQ=255)
                                                                                                                             ttCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGcc                                                     >  1:222931/1‑120 (MQ=255)
                                                                                                                             ttCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGcc                                                     <  2:222931/120‑1 (MQ=255)
                                                                                                                                     gCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCggcagg                <  4:355208/149‑1 (MQ=255)
                                                                                                                                           gcACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGCCGAGAACCAGTGGCCGCCGCGGCAGGCAAagt          >  8:75220/1‑149 (MQ=255)
                                                                                                                                             aCCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTgg        <  6:352709/149‑1 (MQ=255)
                                                                                                                                             aCCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTgg        <  2:351898/149‑1 (MQ=255)
                                                                                                                                               cGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTa      >  3:372773/1‑149 (MQ=255)
                                                                                                                                               cGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTa      >  1:272359/1‑149 (MQ=255)
                                                                                                                                                 tttGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACa    <  6:266863/149‑1 (MQ=255)
                                                                                                                                                    gctgcGGGCGCACCATAGGGACGCCAAGCTGCCGCAGGATTTCCAGATCCTCGACAGGGACGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACATc  <  5:120457/148‑1 (MQ=255)
                                                                                                                                                    gctgcCTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGCTCCTCGACTGTGAAGACCAGCTACGCCTGCTTACGCGTCTGATCAAAGCCATGAACCTCGACGACAAGCAGTGGCCGCCGCGGCAGGCAATGCGGTACATc  <  4:88261/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGACGCCAATCTGCCGCAGGATTTCCAGATCCTCGACAGTGAAGACCAGCTACGCCTGCTTAAGCGTCTGATCAAAGCCATGAACCTCGACGAGAAGCAGTGGCCGCCGCGGCAGGCAATGTGGTACATC  >  NC_000913/3998125‑3998420

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