Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,961,251 A→T 6.6% L393Q (CTG→CAG)  recC ← exonuclease V (RecBCD complex), gamma chain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,961,2510AT 6.6% 324.7 / 7.6 115L393Q (CTG→CAG) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base A (57/50);  new base T (4/4);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAGTTCC  >  NC_000913/2961105‑2961400
                                                                                                                                                  |                                                                                                                                                     
ggTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGACAGCTGGc                                                                                                                                                    >  1:302605/1‑150 (MQ=255)
ggTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGAGATCAGCCACCATCACGATGATGTCGCGCGGAGTAATTGTCGGGTGTTCCTCCAGCATCGCACGCAggg                                                                                                                                                    >  2:248677/1‑149 (MQ=255)
 gTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGAGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCg                                                                                                                                                   >  1:1317648/1‑150 (MQ=255)
   cGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGAt                                                                                                                                                 <  2:29162/150‑1 (MQ=255)
     gAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAACCGGACTGTAGCTGCCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCg                                                                                                                                               <  1:1116901/150‑1 (MQ=255)
     gAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCg                                                                                                                                               <  2:1700848/150‑1 (MQ=255)
     gAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCg                                                                                                                                               <  1:652879/150‑1 (MQ=255)
     gAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCg                                                                                                                                               <  1:1399876/150‑1 (MQ=255)
       actGGCGTAAGGTAGGTAACCATCCGCAGGTGCGCTACCAAACACCGCCTGACTAACCGGACTGTAGCTGGCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGGCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAtgctg                                                                                                                                             <  1:313530/148‑5 (MQ=255)
       aagggCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGGCTTCCTCCAGCATCGCCTGCTGGCGAtgctg                                                                                                                                             <  1:1131957/147‑5 (MQ=255)
       aagggCGTAAGCTAGGTAACGATCCGCAGGCGCGCTACCAACCACAGCCTGAATAAACGGACTGTACCTGGCGATAGCAGCCACCATCACGATGATGGCGCGCGGAGTAAGTGTCGGGGCTTCCTCCAGCATCGCCAGCAGGCGATCgtg                                                                                                                                             <  1:1283809/147‑1 (MQ=255)
         tGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTa                                                                                                                                           >  2:779309/1‑150 (MQ=255)
           gCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGGTATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTaaa                                                                                                                                         >  1:1438819/1‑150 (MQ=255)
                aggtaggtAACGATCCGCATGTGCGCTACCAAACCCAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttc                                                                                                                                    <  2:1356094/150‑1 (MQ=255)
                aggtaggtAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttc                                                                                                                                    >  1:803548/1‑150 (MQ=255)
                   taggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaac                                                                                                                                 <  1:1387337/150‑1 (MQ=255)
                    aggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaact                                                                                                                                <  2:1081850/150‑1 (MQ=255)
                     ggtaACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCGTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttcaactt                                                                                                                               >  1:289591/1‑150 (MQ=255)
                       taACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAacttcaacttca                                                                                                                             <  2:1019252/150‑1 (MQ=255)
                       taACGATCCGCAGGTGAGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAacttcaacttca                                                                                                                             <  1:1095821/150‑1 (MQ=255)
                        aaCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAACCGGACTGTAGCTGTCGAAATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCAc                                                                                                                            <  1:697615/150‑1 (MQ=255)
                        aaCGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCAc                                                                                                                            <  1:790454/150‑1 (MQ=255)
                          cGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGc                                                                                                                          >  1:421535/1‑150 (MQ=255)
                          cGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGc                                                                                                                          >  2:813991/1‑150 (MQ=255)
                                     tGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCTCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATg                                                                                                               >  2:112546/1‑150 (MQ=255)
                                      gcgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                              <  2:932657/150‑1 (MQ=255)
                                      gcgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                              <  2:622905/150‑1 (MQ=255)
                                       cgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTACAGCATAGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGc                                                                                                             >  2:1269780/1‑150 (MQ=255)
                                       cgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATCATGACGCGCGGAGAAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCTGACGATCGTGTAAAACATCAACTTCACGCTGCGGCCtcagga                                                                                                             >  1:229610/1‑145 (MQ=255)
                                        gcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCGTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCAATGGCa                                                                                                            >  2:1610358/1‑150 (MQ=255)
                                         cTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGGAAAACTTCAACTTCACGCTGCGGGCTATGGCaa                                                                                                           >  2:1299861/1‑150 (MQ=255)
                                           aCCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAAc                                                                                                         >  1:1484933/1‑150 (MQ=255)
                                            ccAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACg                                                                                                        <  2:1236825/150‑1 (MQ=255)
                                               aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTgg                                                                                                     <  2:1795633/150‑1 (MQ=255)
                                               aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTgg                                                                                                     >  1:806055/1‑150 (MQ=255)
                                                acacAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGa                                                                                                    <  1:70494/150‑1 (MQ=255)
                                                    cgCCTGAACAAACGGACTGTAGCTGTCTATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCTCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGt                                                                                                <  2:658951/149‑1 (MQ=255)
                                                     gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                               >  2:4857/1‑150 (MQ=255)
                                                     gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                               >  2:940239/1‑150 (MQ=255)
                                                     gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACATCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                               >  2:558501/1‑150 (MQ=255)
                                                       cTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGAt                                                                                             <  2:806055/150‑1 (MQ=255)
                                                         gAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATAc                                                                                           >  1:538693/1‑150 (MQ=255)
                                                          aaTCAACGGACTGTAGCTGTCGATATCCGCCACCATCACGATGATGTCGCGCGGAGTAAGTGGCGGGGCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACt                                                                                          <  1:914614/150‑1 (MQ=255)
                                                            tAAACGGACTGTAGCTGTCGATATCAGCCACCATCCCGATGATGTCGCGCGGAGTATGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAATACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGc                                                                                        <  1:298294/150‑1 (MQ=255)
                                                                      gTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTgg                                                                              >  1:1068105/1‑150 (MQ=255)
                                                                        aGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGAt                                                                            <  2:1562548/150‑1 (MQ=255)
                                                                          cTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCa                                                                          >  2:1817532/1‑150 (MQ=255)
                                                                          cTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCa                                                                          <  2:1653324/150‑1 (MQ=255)
                                                                          cTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCa                                                                          <  2:1484933/150‑1 (MQ=255)
                                                                            gTCGATATCAGCCACCATCACGATGATGTCGCTCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCCAACGTGGAAGGTGATACTGCTATCCAGTTGATCAAg                                                                        >  2:1389418/1‑150 (MQ=255)
                                                                             tCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCGTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGc                                                                       <  1:925349/150‑1 (MQ=255)
                                                                              cggtatCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGCTCGTGTAAAACTTGAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCg                                                                      >  2:84371/4‑150 (MQ=255)
                                                                               gATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCgg                                                                     <  1:1368452/150‑1 (MQ=255)
                                                                               gATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCgg                                                                     <  1:1136581/150‑1 (MQ=255)
                                                                               gATATCAGACACAATCTCGTTGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCgg                                                                     <  2:1669407/150‑1 (MQ=255)
                                                                                 tatCAGCCACCACCACGATGCTGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCACGCTATCGTGTAAAACTTCCACTTACCGCTGAGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGAGCAAGCGGGc                                                                   >  2:72417/1‑150 (MQ=255)
                                                                                   tCAGCCACCATCACGATGATGTCGCGCGGTGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATTCTGCTATCCAGAGGATCAAGCGGGCGt                                                                 >  2:362528/1‑150 (MQ=255)
                                                                                    cAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGtt                                                                <  2:1677089/150‑1 (MQ=255)
                                                                                      gCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTg                                                              <  2:185747/150‑1 (MQ=255)
                                                                                        caccaTCACGATGATGTCGCGCGGAGTAAGGGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGtt                                                            <  2:569379/150‑1 (MQ=255)
                                                                                           caTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATc                                                         >  2:276433/1‑150 (MQ=255)
                                                                                            aTCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCg                                                        >  1:414762/1‑150 (MQ=255)
                                                                                                cGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGGAGGCGATCGTGTTAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATGAAGCGGGCGTTTGTTATCGCTAc                                                    >  1:123536/1‑150 (MQ=255)
                                                                                                   tgatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACggg                                                 <  1:1630625/150‑1 (MQ=255)
                                                                                                   tgatgTCGCGCGGAGTAAGTGTCGGGGCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGATGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACggg                                                 <  1:1456607/150‑1 (MQ=255)
                                                                                                    gatgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGa                                                >  2:592304/1‑150 (MQ=255)
                                                                                                     atgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTAAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGaa                                               <  2:1288756/150‑1 (MQ=255)
                                                                                                        tCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAAc                                            <  2:471640/150‑1 (MQ=255)
                                                                                                           ctcgGAGTAAGTGGCGGGGCTGCCTTCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGTGCTATGGCAAACGTGTAAGTTGATACTTCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAActct                                         <  2:1600345/148‑1 (MQ=255)
                                                                                                           cgcgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAActct                                         >  1:536088/1‑150 (MQ=255)
                                                                                                            gcgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCtt                                        <  2:24984/150‑1 (MQ=255)
                                                                                                             cgGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTc                                       >  1:1433952/1‑150 (MQ=255)
                                                                                                                aGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGAt                                    <  2:553965/150‑1 (MQ=255)
                                                                                                                    aGTGTCGGGTCTTCCTCCAGCCTCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGACGGGGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTc                                <  2:83945/150‑1 (MQ=255)
                                                                                                                      tgtCGGGTCGTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGGGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTcac                              <  2:302605/150‑1 (MQ=255)
                                                                                                                       gtCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAACCTTCAACTTAACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCATTCCAGTGGATCAAGCGGGAGTTTGTTATCGCTACGGGAAATCTCATCGATCTTcacc                             >  1:350080/1‑149 (MQ=255)
                                                                                                                       gtCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTcaca                             >  2:928494/1‑150 (MQ=255)
                                                                                                                         cGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCACTGGATCAAGCGGGCGATAGTTATCGCTACGGGAAAACTCTTCGATGTTCACAcc                           >  2:526695/1‑150 (MQ=255)
                                                                                                                          gggTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAAATTCAACTTCACGCTGCGGGCTATGGCAAACGTCGAAGGTGATACTGCTATCCAGTGAATCAAGCGGGCGTTTGTTATCGCTACGCGAAAACTCTACGAAGTTAGCAccc                          >  1:1424145/1‑149 (MQ=255)
                                                                                                                           ggTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAg                         >  1:982658/1‑150 (MQ=255)
                                                                                                                           ggTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACATCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAg                         <  1:125091/150‑1 (MQ=255)
                                                                                                                                 cctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACgg                   >  2:1671900/1‑150 (MQ=255)
                                                                                                                                  ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATTCTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgc                  >  2:57534/1‑150 (MQ=255)
                                                                                                                                  ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACCCCACCAACGgc                  >  2:1268908/1‑150 (MQ=255)
                                                                                                                                  ctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgc                  >  2:960341/1‑150 (MQ=255)
                                                                                                                                   tacacgaTCGCCTGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAATGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgcg                 >  2:1433437/7‑150 (MQ=255)
                                                                                                                                    ccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACATCACGCTGCGGTCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGgcgc                <  2:1319641/150‑1 (MQ=255)
                                                                                                                                     cacgaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGcgcg               >  1:1347698/5‑150 (MQ=255)
                                                                                                                                     cAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCTGTCTATGGCAAACGTGGCAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGcgcg               <  2:1124173/150‑1 (MQ=255)
                                                                                                                                       gCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGt             >  1:275628/1‑150 (MQ=255)
                                                                                                                                         aTCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCCCGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGAGCAAGCGGGGGTTTGTTAGCGCTACGGGAAACCTCTTCGATGTTCCCACCAGCAACGGCGCGGttt           >  1:1650234/1‑150 (MQ=255)
                                                                                                                                           cGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc         >  1:141226/1‑150 (MQ=255)
                                                                                                                                            gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc        >  2:664049/1‑150 (MQ=255)
                                                                                                                                            gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc        <  2:289591/150‑1 (MQ=255)
                                                                                                                                            gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc        >  2:562195/1‑150 (MQ=255)
                                                                                                                                              cagcagGCGATCGTGTAAAACTTCAACTTCCCGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag      <  1:1235132/150‑1 (MQ=255)
                                                                                                                                              cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCTAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag      <  2:538693/150‑1 (MQ=255)
                                                                                                                                              cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag      <  2:1545315/150‑1 (MQ=255)
                                                                                                                                                gcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGGTttccagtt    <  1:1578943/150‑1 (MQ=255)
                                                                                                                                                 cagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCAATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGGTATCGCTACGGGATAACTCTTCGATGTTCACACCAGCACCGGCGCGGTTttccagttc   <  1:722694/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAGTTCC  >  NC_000913/2961105‑2961400

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