Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,467,909 A→G 6.6% S267G (AGC→GGC)  queF → NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,467,9090AG 6.6% 262.7 / 6.2 91S267G (AGC→GGC) queFNADPH‑dependent 7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base G (3/3);  total (47/44)
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.82e-01

CCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATT  >  NC_002947/2467761‑2468057
                                                                                                                                                    |                                                                                                                                                    
ccAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCa                                                                                                                                                      <  1:324737/149‑1 (MQ=255)
 cAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAg                                                                                                                                                     >  6:161074/1‑149 (MQ=255)
  aGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGc                                                                                                                                                    >  4:255534/1‑149 (MQ=255)
    cATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGccc                                                                                                                                                  <  5:297092/149‑1 (MQ=255)
     aTGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCg                                                                                                                                                 >  7:73916/1‑149 (MQ=255)
        ccGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggtg                                                                                                                                              <  1:209837/149‑3 (MQ=255)
        ccGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggtg                                                                                                                                              <  5:340691/149‑3 (MQ=255)
         cGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATa                                                                                                                                             <  6:158720/149‑1 (MQ=255)
         cGACGTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGAt                                                                                                                                              >  8:177450/1‑148 (MQ=255)
             ttCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggt                                                                                                                                               <  2:93458/143‑2 (MQ=255)
             ttCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggt                                                                                                                                               >  1:93458/1‑142 (MQ=255)
             ttCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACa                                                                                                                                          >  7:321594/1‑148 (MQ=255)
               ccATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAAc                                                                                                                                       >  4:314446/1‑149 (MQ=255)
                      cAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATATGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGt                                                                                                                                >  2:338771/1‑149 (MQ=255)
                        gtgtgtGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGGTAACAAACGGCTGGTcc                                                                                                                              >  2:115251/1‑149 (MQ=255)
                           tgtgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCgg                                                                                                                            <  8:156061/148‑1 (MQ=255)
                           tgtgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCgg                                                                                                                            <  3:314446/148‑1 (MQ=255)
                             tgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCa                                                                                                                          <  1:389006/148‑1 (MQ=255)
                             tgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         <  5:161074/149‑1 (MQ=255)
                             tgGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         >  5:373275/1‑149 (MQ=255)
                               gAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGt                                                                                                                        <  3:15578/148‑1 (MQ=255)
                               gAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGATGGTCCGGCAGt                                                                                                                        >  6:111917/1‑148 (MQ=255)
                                    gATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATa                                                                                                                    >  3:101355/1‑147 (MQ=255)
                                    gATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATa                                                                                                                    >  3:112962/1‑147 (MQ=255)
                                        tACCTGCACTTGAAGAATCTGGTGCAGCCGTAGCACCTGAGCGTCTATGCGCGGTATGTGCGCCGGGGTGGGCTGGTCATCAACCGGTATCGCAGCACCGGGCCGTTCAGCGCGGATAACTAACGGCTGGTCCGGCAGTAATACAggg                                                                                                               <  8:390420/148‑1 (MQ=255)
                                         aCCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAggggg                                                                                                             >  7:14873/1‑149 (MQ=255)
                                               aCTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCtt                                                                                                       >  8:216162/1‑149 (MQ=255)
                                                   gaagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGc                                                                                                    <  4:112962/148‑1 (MQ=255)
                                                    aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAggggg                                                                                                             <  2:27629/138‑1 (MQ=255)
                                                    aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAggggg                                                                                                             >  1:27629/1‑138 (MQ=255)
                                                       aaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATAGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGCTCCCGCAGTAATACAGGGGGCCGCTTTGCGGccc                                                                                               >  5:54531/1‑149 (MQ=255)
                                                        aTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGcccc                                                                                              <  8:14873/149‑1 (MQ=255)
                                                        aTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCCTATCGCAGCACCGCGCCGATCAGCCCGGATAACAAACGGCTCGTGCGGCAGTAATACAGGGGGCCGCATTGCGGcccc                                                                                              >  2:426223/1‑149 (MQ=255)
                                                             ctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCTATCAGCCCGTATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCg                                                                                            <  1:190831/146‑1 (MQ=255)
                                                             ctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCg                                                                                            <  3:270319/146‑1 (MQ=255)
                                                             ctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGc                                                                                           <  1:385819/147‑1 (MQ=255)
                                                              tgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACTGGGGGCCGCTTTGCGGCCCCCGCTAt                                                                                        <  8:73916/149‑1 (MQ=255)
                                                                cAGCGGGAGCACGTGAGGTTTTTTGCGCGCTTTGTGCGCGGGGGTGGGCTGGACATATAGCCTTATCGCAGCTCGGGGCCTATCAGCCCGTATAACAATGGGCTGGTCGGGCAGTATTACTGGGGGCCGTTTTGCGGCCCCCGCTAttt                                                                                      <  6:54531/149‑1 (MQ=255)
                                                                cAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCGGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAt                                                                                        <  6:373275/147‑1 (MQ=255)
                                                                  gCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    <  3:397141/149‑1 (MQ=255)
                                                                   ccGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAttt                                                                                      >  1:99969/1‑146 (MQ=255)
                                                                   ccGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAg                                                                                   <  7:177450/149‑1 (MQ=255)
                                                                    cGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGa                                                                                  <  3:255534/149‑1 (MQ=255)
                                                                          aCCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCa                                                                            >  4:132910/1‑149 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         <  1:216642/100‑1 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                         >  2:216642/1‑100 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  3:340494/1‑149 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  6:290504/1‑149 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                        >  4:34660/1‑149 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCCGCGGGGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGTATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    <  6:225829/137‑1 (MQ=255)
                                                                              gACGGTCTATGCGCGCTATGTGCGCAGCGGTGGGCTCGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    >  5:225829/1‑137 (MQ=255)
                                                                               aCGGTCTATGCGCGCTATGTGCGCCGAGGTGGGCTGGACATCAACACGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCACGCTATTTCAGATCCCCATGCtg                                                                       >  2:222413/1‑149 (MQ=255)
                                                                                 ggTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGGTAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                     <  1:115251/149‑1 (MQ=255)
                                                                                       tGCGCGCTATGTGCGCCGCGGTGGGCTGTACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGTATAACAAACGGCTGGTCCGTCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGc                                                                 <  1:338771/147‑1 (MQ=255)
                                                                                           cgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGGTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTg                                                           <  4:101355/149‑1 (MQ=255)
                                                                                                tgtgCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAAt                                                      <  7:216162/149‑1 (MQ=255)
                                                                                                    cgccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGtt                                                            >  7:158767/1‑139 (MQ=255)
                                                                                                    cgccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGtt                                                            <  8:158767/139‑1 (MQ=255)
                                                                                                     gccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGccc                                                                                               >  3:218259/1‑103 (MQ=255)
                                                                                                     gccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGccc                                                                                               <  4:218259/103‑1 (MQ=255)
                                                                                                       cgcgGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTgcgc                                               <  3:134399/149‑1 (MQ=255)
                                                                                                           gTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAggg                                           <  8:321594/149‑1 (MQ=255)
                                                                                                            tGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTTAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGt                                          >  2:359893/1‑149 (MQ=255)
                                                                                                            tGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGt                                          <  4:340494/149‑1 (MQ=255)
                                                                                                            tGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGACCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGt                                          <  3:367592/149‑1 (MQ=255)
                                                                                                                   taCATCAACCCGTTTCGCAGCACCGGGCCGTTCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACTGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATg                                                                          <  4:159749/109‑1 (MQ=255)
                                                                                                                   gACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATg                                                                          >  3:159749/1‑110 (MQ=255)
                                                                                                                   gACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGt                                    >  7:73345/1‑148 (MQ=255)
                                                                                                                    aCATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGa                                  <  1:169547/149‑1 (MQ=255)
                                                                                                                    aCATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGa                                  <  8:73345/149‑1 (MQ=255)
                                                                                                                      aTCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAgg                                >  8:328946/1‑149 (MQ=255)
                                                                                                                         aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGTCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                             >  6:194611/1‑149 (MQ=255)
                                                                                                                         aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                             >  5:174872/1‑149 (MQ=255)
                                                                                                                         aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCGATTTCAGATCCCCAGGCTGTGCAGCGAGTGGGCAATGCTGCGCAGGCGGGCAGTGAGGGCg                             >  5:439644/1‑149 (MQ=255)
                                                                                                                           cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    >  5:348932/1‑92 (MQ=255)
                                                                                                                           cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                                    <  6:348932/92‑1 (MQ=255)
                                                                                                                           cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCggg                           >  1:333464/1‑149 (MQ=255)
                                                                                                                           cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                            >  7:147349/1‑148 (MQ=255)
                                                                                                                           cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                            <  8:147349/148‑1 (MQ=255)
                                                                                                                              gTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGggtg                        >  3:376006/1‑149 (MQ=255)
                                                                                                                              gTATCGCAGCACCGGGCCGATCAGCCCGGATAAAAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGACCCCGCTATTTCAGATCACCATGCTGTGCATCTAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGggtg                        >  2:75542/1‑149 (MQ=255)
                                                                                                                                 tCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTc                     <  2:333464/149‑1 (MQ=255)
                                                                                                                                    cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc                  >  6:153452/1‑149 (MQ=255)
                                                                                                                                    cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc                  >  7:308604/1‑149 (MQ=255)
                                                                                                                                     ttatccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGcc                 >  5:248662/5‑149 (MQ=255)
                                                                                                                                        tccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTc               >  7:117177/2‑148 (MQ=255)
                                                                                                                                        aCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCGGGCAGTAATACAGGGGGCCTCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTTTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              <  5:111917/149‑1 (MQ=255)
                                                                                                                                        aCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              <  4:376006/149‑1 (MQ=255)
                                                                                                                                          cGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              >  2:348998/1‑147 (MQ=255)
                                                                                                                                          cGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg              >  8:271128/1‑147 (MQ=255)
                                                                                                                                          cGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGaa            <  1:391555/149‑1 (MQ=255)
                                                                                                                                            ggCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAgc          >  5:433902/1‑149 (MQ=255)
                                                                                                                                                    aGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGAtt  <  5:290504/149‑1 (MQ=255)
                                                                                                                                                    aGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGAtt  >  2:222619/1‑149 (MQ=255)
                                                                                                                                                    aGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGAtt  <  1:222619/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATT  >  NC_002947/2467761‑2468057

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