Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,950,165 A→T 7.1% L367Q (CTG→CAG)  oguA ← 8‑oxoguanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,950,1650AT 7.1% 311.1 / 8.9 100L367Q (CTG→CAG) oguA8‑oxoguanine deaminase
Reads supporting (aligned to +/- strand):  ref base A (48/44);  new base T (4/3);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GGGGTGCAGAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCC  >  NC_002947/2950021‑2950313
                                                                                                                                                  |                                                                                                                                                    
ggggTGCAGAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTg                                                                                                                                                 <  2:1871496/150‑1 (MQ=255)
   gTGCAGAACTCCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCa                                                                                                                                              <  1:1217083/150‑1 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTTAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtctc                                                                                                                                         <  1:1774709/150‑1 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCTGGGTGGcgcgaaa                                                                                                                                         >  2:2232738/1‑144 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCTGGGTGGcgcgaaa                                                                                                                                         <  2:2132582/150‑7 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCTGGGTGGcgcgaaa                                                                                                                                         <  1:2232738/150‑7 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGt                                                                                                                                                  >  1:2075941/1‑141 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGt                                                                                                                                                  <  2:2075941/141‑1 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtctc                                                                                                                                         >  2:1662850/1‑150 (MQ=255)
        gAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGCCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCTGGGTGGcgcgaaa                                                                                                                                         >  1:2132582/1‑144 (MQ=255)
         aaCACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGCAGCCGGCCGGGGCGCGGTGTCCAAGATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTcc                                                                                                                                        >  1:281444/1‑150 (MQ=255)
           cacCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      >  2:2235965/1‑150 (MQ=255)
           cacCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      >  2:2235167/1‑150 (MQ=255)
           cacCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      <  1:487938/150‑1 (MQ=255)
           cacCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      <  1:1662850/150‑1 (MQ=255)
                ggggtcgatggctctgcctccaacgacggcgccagcatgatcggtgaagtgcgtcaggccctgctgctgcaacgggtgggctttggccctgacgcgatgaccgcacgcgAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCtt                                                                                                                                 >  1:103544‑M1/110‑150 (MQ=255)
                   gtcgatggctctgcctccaacgacggcgccagcatgatcggtgaagtgcgtcaggccctgctgctgcaacgggtgggctttggccctgacgcgatgaccgcacgcgAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTcgc                                                                                                                              >  1:732092‑M1/107‑150 (MQ=255)
                     gccAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGt                                                                                                                            <  1:657255/150‑1 (MQ=255)
                     gccAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGAGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGt                                                                                                                            >  2:1252696/1‑150 (MQ=255)
                     gccAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCAACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGt                                                                                                                            <  1:1633947/150‑1 (MQ=255)
                        tggctctgcctccaacgacggcgccagcatgatcggtgaagtgcgtcaggccctgctgctgcaacgggtgggctttggccctgacgcgatgaccgcacgcgAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCg                                                                                                                         >  2:1838000‑M1/102‑150 (MQ=255)
                        aaTGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtctc                                                                                                                                         <  2:714971/134‑1 (MQ=255)
                        aaTGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtctc                                                                                                                                         >  1:714971/1‑134 (MQ=255)
                        aaTGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCg                                                                                                                         <  1:2415457/150‑1 (MQ=255)
                            ggTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCa                                                                                                                     >  2:1938920/1‑150 (MQ=255)
                                gtggtgGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGAGTGCGGTCATcgc                                                                                                                 <  1:1733618/150‑1 (MQ=255)
                                      gccgccGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAgg                                                                                                           >  1:827431/1‑150 (MQ=255)
                                      gccgccGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAgg                                                                                                           <  1:2317346/150‑1 (MQ=255)
                                       ccgccgGCATACGCCACATGGCCCAGGTCGAAGGCCATAAATTCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAggg                                                                                                          <  1:723203/150‑1 (MQ=255)
                                        cgccgGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGc                                                                                                         <  1:1499379/150‑1 (MQ=255)
                                        cgccgGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGc                                                                                                         >  1:487981/1‑150 (MQ=255)
                                        cgccgGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGc                                                                                                         <  1:2047231/150‑1 (MQ=255)
                                                   gCCACATTGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAc                                                                                               <  1:636233/149‑1 (MQ=255)
                                                   gCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtct                                                                                                                                          >  1:981804/1‑106 (MQ=255)
                                                   gCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAAtct                                                                                                                                          <  2:981804/106‑1 (MQ=255)
                                                   gCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAcc                                                                                              >  1:1329076/1‑150 (MQ=255)
                                                   gCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAcc                                                                                              >  1:185441/1‑150 (MQ=255)
                                                   gCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAc                                                                                               >  2:636233/1‑149 (MQ=255)
                                                    ccACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAccc                                                                                             >  2:1198524/1‑150 (MQ=255)
                                                      acaTGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGt                                                                                           >  1:1218628/1‑150 (MQ=255)
                                                        aTGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTg                                                                                         >  2:200026/1‑150 (MQ=255)
                                                          ggCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTgca                                                                                       <  1:2083923/150‑1 (MQ=255)
                                                          ggCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTgca                                                                                       <  1:848508/150‑1 (MQ=255)
                                                             ccAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATc                                                                                                                   >  2:935548/1‑119 (MQ=255)
                                                             ccAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATc                                                                                                                   <  1:935548/119‑1 (MQ=255)
                                                              cAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTgcagcag                                                                                   >  2:740884/1‑150 (MQ=255)
                                                              cAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTgcagcag                                                                                   >  1:1186908/1‑150 (MQ=255)
                                                                    gAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGc                                                                             <  1:1936294/150‑1 (MQ=255)
                                                                    gAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAACCCCACCCGTTGCAGCAGCAGGGc                                                                             <  2:1186908/150‑1 (MQ=255)
                                                                            cgaAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCa                                                                     >  2:828320/1‑150 (MQ=255)
                                                                             gaAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCg                                                                                                                         <  2:1357546/97‑1 (MQ=255)
                                                                             gaAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCg                                                                                                                         >  1:1357546/1‑97 (MQ=255)
                                                                                  cGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCAc                                                               <  1:828320/150‑1 (MQ=255)
                                                                                      gACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      >  2:1417359/1‑75 (MQ=255)
                                                                                      gACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAg                                                                                                                                      <  1:1417359/75‑1 (MQ=255)
                                                                                       aCCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAcc                                                                                              <  2:2248031/114‑1 (MQ=255)
                                                                                       aCCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCAcc                                                                                              >  1:2248031/1‑114 (MQ=255)
                                                                                       aCCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATc                                                          <  1:740884/150‑1 (MQ=255)
                                                                                         cATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGt                                                                                                               <  2:532031/95‑1 (MQ=255)
                                                                                         cATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGt                                                                                                               >  1:532031/1‑95 (MQ=255)
                                                                                          aTGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTGCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATg                                                       <  2:456452/150‑1 (MQ=255)
                                                                                            gCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCt                                                     <  1:1539210/150‑1 (MQ=255)
                                                                                              ccGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTgg                                                   >  1:1612843/1‑150 (MQ=255)
                                                                                              ccGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTgg                                                   >  1:1978886/1‑150 (MQ=255)
                                                                                                 gggCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGgcgc                                                <  1:1683210/150‑1 (MQ=255)
                                                                                                  ggCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGcc                                               <  1:1252696/150‑1 (MQ=255)
                                                                                                    ccAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCATGACGCACTTCACCGATCATGCTGGCGCcgt                                             >  2:1246441/1‑150 (MQ=255)
                                                                                                       cctgacgcgatgaccgcacgcgAAGCACTGGAGATTGCCACCCTGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACGCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCcgtc                                            <  1:1119731‑M1/128‑1 (MQ=255)
                                                                                                        ggCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGtt                                         >  2:399263/1‑150 (MQ=255)
                                                                                                         gcgtCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCGGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTg                                        <  2:1407352/150‑1 (MQ=255)
                                                                                                         gcgcCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTg                                        >  1:1112386/1‑150 (MQ=255)
                                                                                                         gcgcCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTg                                        >  1:1781261/1‑150 (MQ=255)
                                                                                                           gcCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGa                                      <  2:185441/150‑1 (MQ=255)
                                                                                                             cAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAgg                                    >  1:2240566/1‑150 (MQ=255)
                                                                                                              aaTATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGc                                   >  2:1869629/1‑150 (MQ=255)
                                                                                                                 atCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCaga                                <  1:2400704/150‑1 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCAt                                                        >  2:901486/1‑117 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCAt                                                        <  1:901486/117‑1 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAc                       >  1:56719/1‑150 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAc                       >  2:1050630/1‑150 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAc                       >  2:2268700/1‑150 (MQ=255)
                                                                                                                          gTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGGCATCGCGTCAGGGCGTAATCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAc                       >  1:1355424/1‑150 (MQ=255)
                                                                                                                             AAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAcccc                    <  1:1494161‑M1/150‑1 (MQ=255)
                                                                                                                              AGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAccccc                   <  2:1864001‑M1/150‑1 (MQ=255)
                                                                                                                              AGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGAccccc                   >  2:1494161‑M1/1‑150 (MQ=255)
                                                                                                                               GCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAgg                                                                               >  1:1814604‑M1/1‑89 (MQ=255)
                                                                                                                               GCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAgg                                                                               <  2:1814604‑M1/89‑1 (MQ=255)
                                                                                                                               GCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCa                  <  1:1629868‑M1/150‑1 (MQ=255)
                                                                                                                               GCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCa                  >  1:1869629‑M1/1‑150 (MQ=255)
                                                                                                                                CACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAg                 >  2:699220‑M1/1‑150 (MQ=255)
                                                                                                                                   TTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGcc              <  1:1915826‑M1/150‑1 (MQ=255)
                                                                                                                                   TTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGcc              <  1:1326088‑M1/150‑1 (MQ=255)
                                                                                                                                     gagattgcCACCCTGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGcc              >  1:2010423/7‑150 (MQ=255)
                                                                                                                                        c‑‑gcCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCATCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGTGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCg         <  1:1453811/150‑1 (MQ=255)
                                                                                                                                           gcCACCCTGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCgg        >  1:1703688/1‑150 (MQ=255)
                                                                                                                                            ccACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGt       <  1:86277/150‑1 (MQ=255)
                                                                                                                                               cccAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACg    >  2:2043954/1‑150 (MQ=255)
                                                                                                                                                ccAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGc   >  1:1645865/1‑150 (MQ=255)
                                                                                                                                                ccAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATAATGGAGAAATCAATTGATTCATAGCCATCGACTCCCAGGCCCACCGGTACga   >  2:1841565/1‑149 (MQ=255)
                                                                                                                                                 cAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGcc  <  2:1781261/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGGGTGCAGAACACCAGTGCCGCCAATGGGTCGTGGTGGCCGCCGGCATACGCCACATGGCCCAGGTCGAAGGCGACGAAATCGGCGACCATGCCCGGGGCCAGGGCGCCAATATCATTGCGGTTAAGCACTTTCGC‑‑GCCACCCAGGGTGGCAATCTCCAGTGCTTCGCGTGCGGTCATCGCGTCAGGGCCAAAGCCCACCCGTTGCAGCAGCAGGGCCTGACGCACTTCACCGATCATGCTGGCGCCGTCGTTGGAGGCAGAGCCATCGACCCCCAGGCCCACCGGTACGCC  >  NC_002947/2950021‑2950313

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑