Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,467,904 C→T 6.3% P265L (CCG→CTG)  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,9040CT 6.3% 282.2 / 4.6 95P265L (CCG→CTG) queFNADPH‑dependent 7‑cyano‑7‑deazaguanine reductase
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base T (3/3);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

TCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGC  >  NC_002947/2467757‑2468049
                                                                                                                                                   |                                                                                                                                                 
tCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGcc                                                                                                                                                   >  4:95064/1‑148 (MQ=255)
   gCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGAt                                                                                                                                                >  4:367592/1‑148 (MQ=255)
    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                                                                                                                                       >  1:93458/1‑142 (MQ=255)
                 ttCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTGATCGGCCCggt                                                                                                                                       <  2:93458/143‑2 (MQ=255)
                 ttCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACa                                                                                                                                  >  7:321594/1‑148 (MQ=255)
                   ccATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAAc                                                                                                                               >  4:314446/1‑149 (MQ=255)
                       gAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATc                                                                                                                                               >  4:327968/1‑129 (MQ=255)
                       gAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATc                                                                                                                                               <  3:327968/129‑1 (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:112962/1‑147 (MQ=255)
                                        gATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATa                                                                                                            >  3:101355/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                                                                                                     >  1:27629/1‑138 (MQ=255)
                                                        aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAggggg                                                                                                     <  2:27629/138‑1 (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                                                                                                                 >  2:216642/1‑100 (MQ=255)
                                                                                  gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAg                                                                                                                 <  1:216642/100‑1 (MQ=255)
                                                                                  gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                >  6:290504/1‑149 (MQ=255)
                                                                                  gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                >  4:34660/1‑149 (MQ=255)
                                                                                  gACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCt                                                                >  3:340494/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                                                    <  8:158767/139‑1 (MQ=255)
                                                                                                        cgccgcGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGtt                                                    >  7:158767/1‑139 (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                                                                            <  6:348932/92‑1 (MQ=255)
                                                                                                                               cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCa                                                                            >  5:348932/1‑92 (MQ=255)
                                                                                                                               cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCggg                   >  1:333464/1‑149 (MQ=255)
                                                                                                                               cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                    <  8:147349/148‑1 (MQ=255)
                                                                                                                               cccGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCgg                    >  7:147349/1‑148 (MQ=255)
                                                                                                                                  gTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGggtg                >  3:376006/1‑149 (MQ=255)
                                                                                                                                  gTATCGCAGCACCGGGCCGATCAGCCCGGATAAAAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGACCCCGCTATTTCAGATCACCATGCTGTGCATCTAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGggtg                >  2:75542/1‑149 (MQ=255)
                                                                                                                                     tCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTc             <  2:333464/149‑1 (MQ=255)
                                                                                                                                        cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc          >  7:308604/1‑149 (MQ=255)
                                                                                                                                        cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc          >  6:153452/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      >  8:271128/1‑147 (MQ=255)
                                                                                                                                              cGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCg      >  2:348998/1‑147 (MQ=255)
                                                                                                                                              cGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGaa    <  1:391555/149‑1 (MQ=255)
                                                                                                                                                ggCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAgc  >  5:433902/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCCGCCAGCATGCCGACTTCCATGAGCAGTGTGTGGAGCGGATTTACCTGGACTTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGC  >  NC_002947/2467757‑2468049

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: