Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,036,655 C→T 5.2% L484L (CTG→CTA recJ ← ssDNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,036,6550CT 5.2% 244.7 / 2.0 76L484L (CTG→CTArecJssDNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base C (38/34);  new base T (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

GCAATACCATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACA  >  NC_000913/3036506‑3036803
                                                                                                                                                     |                                                                                                                                                    
gCAATACCATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCgcagc                                                                                                                                                      >  2:1154359/1‑150 (MQ=255)
   aTACCATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCagcagc                                                                                                                                                   >  2:1133042/1‑150 (MQ=255)
    tACCATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCt                                                                                                                                                  >  2:1143513/1‑150 (MQ=255)
      ccATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTgc                                                                                                                                                >  2:1097097/1‑150 (MQ=255)
        aTCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGTAGCTgcgc                                                                                                                                              >  2:712682/1‑150 (MQ=255)
        aTCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGTAGCTgcgc                                                                                                                                              >  2:707582/1‑150 (MQ=255)
        aTCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTgcgc                                                                                                                                              >  1:254652/1‑150 (MQ=255)
         tCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGcc                                                                                                                                             <  2:1601198/150‑1 (MQ=255)
         tCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGcc                                                                                                                                             <  2:142988/150‑1 (MQ=255)
          ccAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCa                                                                                                                                            <  1:1097097/150‑1 (MQ=255)
           cagcagTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCTGCTgcgcgat                                                                                                                                           <  2:782391/150‑4 (MQ=255)
           cagcagTGGACCGCCGCCGACCGGTTCGACCATAACCTTCAAATGACGTTCGCCCACCAGCCGCTGATGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGcccc                                                                                                                                           >  1:744809/1‑148 (MQ=255)
                gTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCa                                                                                                                                      >  1:1326368/1‑150 (MQ=255)
                     ccgccgccgACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGTGGCTCCGGGAACATCTGCCCCCCCGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTc                                                                                                                                 >  2:571817/1‑150 (MQ=255)
                      cgccgccgACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCa                                                                                                                                <  2:1187256/150‑1 (MQ=255)
                         cgccgACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTgcg                                                                                                                                               <  2:1139720/132‑1 (MQ=255)
                         cgccgACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTgcg                                                                                                                                               >  1:1139720/1‑132 (MQ=255)
                          gccgACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTc                                                                                                                            >  2:791804/1‑150 (MQ=255)
                             gACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGc                                                                                                                         >  1:910002/1‑150 (MQ=255)
                             gACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGc                                                                                                                         <  1:801383/150‑1 (MQ=255)
                               ccGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCg                                                                                                                       >  2:270415/1‑150 (MQ=255)
                                    tCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCtt                                                                                                                  >  1:1555443/1‑150 (MQ=255)
                                      gACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTaa                                                                                                                >  2:1237661/1‑150 (MQ=255)
                                          aTCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGa                                                                                                            >  1:601183/1‑150 (MQ=255)
                                             aCCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCGCGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCg                                                                                                         <  2:1286582/150‑1 (MQ=255)
                                               cTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTc                                                                                                       <  1:791804/150‑1 (MQ=255)
                                               cTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTc                                                                                                       >  1:1562731/1‑150 (MQ=255)
                                                  cAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGa                                                                                                    <  1:1157765/150‑1 (MQ=255)
                                                  cAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGa                                                                                                    >  1:1173222/1‑150 (MQ=255)
                                                   aaaTGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGAt                                                                                                   <  1:1133042/150‑1 (MQ=255)
                                                     aTGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATAc                                                                                                 >  1:77723/1‑150 (MQ=255)
                                                      tGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATAcc                                                                                                <  1:676256/150‑1 (MQ=255)
                                                       gACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCa                                                                                               >  1:1067484/1‑150 (MQ=255)
                                                              gCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGcc                                                                                        >  2:335788/1‑150 (MQ=255)
                                                               cccACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCt                                                                                       <  1:629519/150‑1 (MQ=255)
                                                                ccaccaGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCtt                                                                                      >  2:1336389/1‑150 (MQ=255)
                                                                  accaGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCGCTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGc                                                                                    >  2:281793/1‑150 (MQ=255)
                                                                     aGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAAt                                                                                 <  1:467156/150‑1 (MQ=255)
                                                                     aGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAAt                                                                                 <  1:26757/150‑1 (MQ=255)
                                                                     aGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACAGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAAt                                                                                 <  1:897164/150‑1 (MQ=255)
                                                                       ccGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAg                                                                               >  2:1288987/1‑150 (MQ=255)
                                                                          cTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGGCACTTCCATGGTGATTTCGGCCGGGCTTAACGGACCGTCTGATACAACTTCGCCTTGCATTAGCGa                                                                            <  2:876500/150‑1 (MQ=255)
                                                                            gTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAg                                                                          >  2:1514466/1‑150 (MQ=255)
                                                                            gTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAg                                                                          <  2:1457326/150‑1 (MQ=255)
                                                                               gcagcaGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCGCTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGt                                                                       <  1:281793/150‑1 (MQ=255)
                                                                               gcagcaGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGt                                                                       >  1:1245898/1‑150 (MQ=255)
                                                                                 agcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTcc                                                                     >  2:1077300/1‑150 (MQ=255)
                                                                                  gcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCa                                                                    >  1:804599/1‑150 (MQ=255)
                                                                                  gcagACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCa                                                                    <  2:254652/150‑1 (MQ=255)
                                                                                    agACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGc                                                                  <  1:87562/150‑1 (MQ=255)
                                                                                        ggAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCAcg                                                              >  1:507853/1‑149 (MQ=255)
                                                                                        ggAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACt                                                              <  2:1562731/150‑1 (MQ=255)
                                                                                         gAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGACGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTc                                                             <  2:128569/150‑1 (MQ=255)
                                                                                                ccttCAAACAGCGGCTCAGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAAcc                                                      >  2:1302289/4‑150 (MQ=255)
                                                                                                ccGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAAcc                                                      <  1:270415/150‑1 (MQ=255)
                                                                                                ccGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAAcc                                                      <  1:1091850/150‑1 (MQ=255)
                                                                                                 cGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCa                                                     >  2:1061794/1‑150 (MQ=255)
                                                                                                  gTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAg                                                    >  1:569415/1‑150 (MQ=255)
                                                                                                     aaaCAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTc                                                 >  2:1169731/1‑150 (MQ=255)
                                                                                                               tCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGGGATGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCgt                                       <  1:1293798/150‑1 (MQ=255)
                                                                                                                ccGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCTACTTGCAATAGCGAAGGGTTCAGCAACTCAGTAACAATTTCGTTAAATCgtt                                      <  2:1413830/150‑1 (MQ=255)
                                                                                                                 cGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGTAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCgttg                                     <  1:707582/150‑1 (MQ=255)
                                                                                                                 cGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGTAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCgttg                                     <  1:712682/150‑1 (MQ=255)
                                                                                                                        aTCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAag                              <  1:931219/150‑1 (MQ=255)
                                                                                                                        aTCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAag                              >  1:295187/1‑150 (MQ=255)
                                                                                                                         tCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAaga                             >  1:1289196/1‑150 (MQ=255)
                                                                                                                            gCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCACACCGTTGTTGAAAGAGtt                          >  1:1445197/1‑150 (MQ=255)
                                                                                                                            gCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGtt                          >  2:1260647/1‑150 (MQ=255)
                                                                                                                            gCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGtt                          >  1:374783/1‑150 (MQ=255)
                                                                                                                                 cACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAAt                     <  2:1580793/150‑1 (MQ=255)
                                                                                                                                 cACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAAt                     <  1:675728/150‑1 (MQ=255)
                                                                                                                                    gggCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTa                  <  1:1154359/150‑1 (MQ=255)
                                                                                                                                     ggCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTAt                 >  2:51461/1‑150 (MQ=255)
                                                                                                                                     ggCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTAt                 <  1:571817/150‑1 (MQ=255)
                                                                                                                                      gCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATc                <  2:910002/150‑1 (MQ=255)
                                                                                                                                       ccAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATcc               >  1:1432799/1‑150 (MQ=255)
                                                                                                                                         gaagtggcgcAGCTGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCctc             >  2:680874/7‑150 (MQ=255)
                                                                                                                                                cgcAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGc      <  2:1245898/150‑1 (MQ=255)
                                                                                                                                                 gcagcaGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCg     <  1:1143513/150‑1 (MQ=255)
                                                                                                                                                   agcagcTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGAc   <  1:1237661/150‑1 (MQ=255)
                                                                                                                                                   agcagcTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGAc   <  1:1077300/150‑1 (MQ=255)
                                                                                                                                                    gcagcTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACa  >  2:1022471/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAATACCATCCAGCAGTGGACCGCCGCCGACCGGTTCGACCATCACCTTCAAATGACGTTCGCCCACCAGCCGCTGTTGCAGCAGACGGAAATGACCGTCAAACAGCGGCTCCGGGAACATCTGCCCCCACGGGCCAGCATCGCGCAGCAGCTGCGCCACTTCCATGGTCATTTCGGCCGGGCTTAACGGACCGTCTGATACCACTTCGCCTTGCAATAGCGAAGGGTCCAGCCACTCAGTAACCAGTTCGCCAAACCGTTGTTGAAAGAGTTTGAATTTATCCTCTTCCAGCGACA  >  NC_000913/3036506‑3036803

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