Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,961,248 A→T 6.1% L394Q (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,2480AT 6.1% 310.7 / 5.3 115L394Q (CTG→CAG) recCexonuclease V (RecBCD complex), gamma chain
Reads supporting (aligned to +/- strand):  ref base A (57/49);  new base T (3/5);  total (61/54)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

ACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAG  >  NC_000913/2961100‑2961396
                                                                                                                                                    |                                                                                                                                                    
 cgacgGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCcagc                                                                                                                                                    >  2:1368452/1‑150 (MQ=255)
  gacgGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCctgct                                                                                                                                                   <  1:520455/150‑2 (MQ=255)
     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)
               ggCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGATTATACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGCTCTCGCGCGGAGTAAGTGTCTCGTCTTCCTCCAGCATCGCCAGCACGCGAACGTATaa                                                                                                                                      >  1:338044/1‑150 (MQ=255)
               ggCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACAATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCTGCTGGCGAtgctggag                                                                                                                                      <  1:974754/150‑8 (MQ=255)
                gCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGGTATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTaaa                                                                                                                                     >  1:1438819/1‑150 (MQ=255)
                     aggtaggtAACGATCCGCATGTGCGCTACCAAACCCAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAaacttc                                                                                                                                <  2:1356094/150‑1 (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                                                                                                                      >  2:813991/1‑150 (MQ=255)
                               cGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGc                                                                                                                      >  1:421535/1‑150 (MQ=255)
                                          tGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCTCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATg                                                                                                           >  2:112546/1‑150 (MQ=255)
                                           gcgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                          <  2:622905/150‑1 (MQ=255)
                                           gcgcTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATgg                                                                                                          <  2:932657/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)
                                                    aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTgg                                                                                                 >  1:806055/1‑150 (MQ=255)
                                                    aaCACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTgg                                                                                                 <  2:1795633/150‑1 (MQ=255)
                                                     acacAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGa                                                                                                <  1:70494/150‑1 (MQ=255)
                                                         cgCCTGAACAAACGGACTGTAGCTGTCTATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCTCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGt                                                                                            <  2:658951/149‑1 (MQ=255)
                                                          gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                           >  2:940239/1‑150 (MQ=255)
                                                          gCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTg                                                                                           >  2:4857/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:1484933/150‑1 (MQ=255)
                                                                               cTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCa                                                                      <  2:1653324/150‑1 (MQ=255)
                                                                                 gTCGATATCAGCCACCATCACGATGATGTCGCTCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCCAACGTGGAAGGTGATACTGCTATCCAGTTGATCAAg                                                                    >  2:1389418/1‑150 (MQ=255)
                                                                                  tCGATATCAGCCACCATCACGTTGATGTCGCGCGTAGTAAGTGTCGTGTCTTCCTCCAGCATCGCCTGCAGGCGATCGTGTAAAACTTCTACTTCACTCTGCGGTCTATGGCAAACGTGTTATGTGATACTGCTATCCATTTGATCAAGc                                                                   <  2:1317648/150‑1 (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)
                                                                                                           tgTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGaaa                                          >  1:246328/1‑150 (MQ=255)
                                                                                                            gTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGaaaa                                         >  1:798656/1‑150 (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)
                                                                                                                            gtCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTcaca                         >  2:928494/1‑150 (MQ=255)
                                                                                                                              cGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCACTGGATCAAGCGGGCGATAGTTATCGCTACGGGAAAACTCTTCGATGTTCACAcc                       >  2:526695/1‑150 (MQ=255)
                                                                                                                               gggTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAAATTCAACTTCACGCTGCGGGCTATGGCAAACGTCGAAGGTGATACTGCTATCCAGTGAATCAAGCGGGCGTTTGTTATCGCTACGCGAAAACTCTACGAAGTTAGCAccc                      >  1:1424145/1‑149 (MQ=255)
                                                                                                                                ggTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACATCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAg                     <  1:125091/150‑1 (MQ=255)
                                                                                                                                      cctccAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACgg               >  2:1184727/1‑150 (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)
                                                                                                                                          cacgaTCGCCTGCTGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGcgcg           >  1:1347698/5‑150 (MQ=255)
                                                                                                                                          cAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCTGTCTATGGCAAACGTGGCAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGcgcg           <  2:1124173/150‑1 (MQ=255)
                                                                                                                                              aTCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCCCGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGAGCAAGCGGGGGTTTGTTAGCGCTACGGGAAACCTCTTCGATGTTCCCACCAGCAACGGCGCGGttt       >  1:1650234/1‑150 (MQ=255)
                                                                                                                                                cGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttc     >  1:141226/1‑150 (MQ=255)
                                                                                                                                                 gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc    >  2:562195/1‑150 (MQ=255)
                                                                                                                                                 gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc    >  2:664049/1‑150 (MQ=255)
                                                                                                                                                 gCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttcc    <  2:289591/150‑1 (MQ=255)
                                                                                                                                                   cagcagGCGATCGTGTAAAACTTCAACTTCCCGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag  <  1:1235132/150‑1 (MQ=255)
                                                                                                                                                   cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCTAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag  <  2:538693/150‑1 (MQ=255)
                                                                                                                                                   cagcagGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTttccag  <  2:1545315/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGACGGTCGGAAATGGCGTAAGGTAGGTAACGATCCGCAGGTGCGCTACCAAACACAGCCTGAATAAACGGACTGTAGCTGTCGATATCAGCCACCATCACGATGATGTCGCGCGGAGTAAGTGTCGGGTCTTCCTCCAGCATCGCCAGCAGGCGATCGTGTAAAACTTCAACTTCACGCTGCGGGCTATGGCAAACGTGGAAGGTGATACTGCTATCCAGTGGATCAAGCGGGCGTTTGTTATCGCTACGGGAAAACTCTTCGATGTTCACACCAGCAACGGCGCGGTTTTCCAG  >  NC_000913/2961100‑2961396

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: