Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,594,291 T→C 7.6% T392A (ACC→GCC)  kgtP ← alpha‑ketoglutarate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,594,2910TC 7.6% 264.3 / 9.3 92T392A (ACC→GCC) kgtPalpha‑ketoglutarate permease
Reads supporting (aligned to +/- strand):  ref base T (43/42);  new base C (3/4);  total (46/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 = 2.71e-01

TTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCA  >  NC_002947/1594145‑1594439
                                                                                                                                                  |                                                                                                                                                    
ttACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTgcc                                                                                                                                                   <  2:1466173/150‑1 (MQ=255)
ttACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTgcc                                                                                                                                                   <  1:1410745/150‑1 (MQ=255)
  aCTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTgccgc                                                                                                                                                 >  1:1039271/1‑150 (MQ=255)
    tCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGccgccg                                                                                                                                               >  2:703085/1‑150 (MQ=255)
    tCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGccgccg                                                                                                                                               >  1:860541/1‑150 (MQ=255)
    tCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCAccgccg                                                                                                                                               >  2:225162/1‑150 (MQ=255)
      aGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGaa                                                                                                                                             <  1:703085/150‑1 (MQ=255)
       gTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaat                                                                                                                                            <  1:225162/150‑2 (MQ=255)
          gTAATCCGCGAGTGCTGCTTGGTGTATTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATc                                                                                                                                         >  2:356980/1‑150 (MQ=255)
           tAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCCGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGCACCGCCGaatacgt                                                                                                                                        <  2:1184684/150‑6 (MQ=255)
              tCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGaca                                                                                                                                     >  1:481603/1‑150 (MQ=255)
                     gTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCa                                                                                                                               <  1:1609674/149‑1 (MQ=255)
                     gTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCa                                                                                                                               >  2:1609674/1‑149 (MQ=255)
                     gTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAg                                                                                                                              <  1:85941/150‑1 (MQ=255)
                     gTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAg                                                                                                                              <  2:916209/150‑1 (MQ=255)
                      tgctgcTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                             <  2:1109213/150‑1 (MQ=255)
                      ggctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAgg                                                                                                                             <  2:153547/149‑1 (MQ=255)
                       gctgctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAggg                                                                                                                            <  2:860541/150‑1 (MQ=255)
                           ctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                                        >  2:431022/1‑150 (MQ=255)
                           ctTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATa                                                                                                                        >  1:740839/1‑150 (MQ=255)
                            ttGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATAc                                                                                                                       >  2:1497907/1‑150 (MQ=255)
                                tgtCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCa                                                                                                                   <  1:1497907/150‑1 (MQ=255)
                                tgtCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCa                                                                                                                   <  1:92607/150‑1 (MQ=255)
                                    gttCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGcc                                                                                                               >  1:558888/2‑150 (MQ=255)
                                              gcgTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                     >  1:99121/1‑150 (MQ=255)
                                                      cagcagCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAAt                                                                                             <  2:740839/150‑1 (MQ=255)
                                                             gcgCAGGCGATACAGGCGGTGACGTAACAGTAGAAGCCGCTATCGCTCCCGGCGCGCTTGCACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGTCACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTa                                                                                      <  2:496847/148‑1 (MQ=255)
                                                             gAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTa                                                                                      <  2:1126228/150‑1 (MQ=255)
                                                                     gATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAg                                                                              >  2:892852/1‑150 (MQ=255)
                                                                     gATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATGCGGCGGTGCCGAGTAACACTTACAGCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAg                                                                              >  1:779312/1‑150 (MQ=255)
                                                                          aGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACAcc                                                                                                                                   >  1:460081/1‑92 (MQ=255)
                                                                          aGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACAcc                                                                                                                                   <  2:460081/92‑1 (MQ=255)
                                                                              ggTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGa                                                                                                                                       <  2:1416489/84‑1 (MQ=255)
                                                                              ggTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGa                                                                                                                                       >  1:1416489/1‑84 (MQ=255)
                                                                                tGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAGCACCCAGTGCTCGAATTTCGGTGGGGAACAGCTCGGCTTTGa                                                                   >  2:676900/1‑150 (MQ=255)
                                                                                  aCGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     <  2:103545/146‑1 (MQ=255)
                                                                                  aCGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     >  1:103545/1‑146 (MQ=255)
                                                                                    gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                               <  1:431022/150‑1 (MQ=255)
                                                                                    gTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCAc                                                               >  1:1525771/1‑150 (MQ=255)
                                                                                     tACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACg                                                              >  1:1587236/1‑150 (MQ=255)
                                                                                      aCCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         <  1:865055/122‑1 (MQ=255)
                                                                                      aCCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         >  2:865055/1‑122 (MQ=255)
                                                                                       ccAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         >  1:1018/1‑121 (MQ=255)
                                                                                       ccAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         <  2:1018/121‑1 (MQ=255)
                                                                                         aGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGc                                                                                                                           >  2:1520687/1‑85 (MQ=255)
                                                                                         aGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGc                                                                                                                           <  1:1520687/85‑1 (MQ=255)
                                                                                          gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                     >  2:1052885/1‑106 (MQ=255)
                                                                                          gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                     <  1:1052885/106‑1 (MQ=255)
                                                                                          gTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTAGGCACTACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                         >  1:611498/1‑150 (MQ=255)
                                                                                          gTAGAAGCCGCTCTCCACCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGtt                                                         >  2:1198003/1‑150 (MQ=255)
                                                                                            aGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGa                                                       <  2:99121/150‑1 (MQ=255)
                                                                                             gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGaa                                                                                              <  2:344519/110‑1 (MQ=255)
                                                                                             gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGaa                                                                                              >  1:344519/1‑110 (MQ=255)
                                                                                             gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                          >  2:1300848/1‑146 (MQ=255)
                                                                                             gAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGt                                                          <  1:1300848/146‑1 (MQ=255)
                                                                                               aGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCg                                                    >  1:892599/1‑150 (MQ=255)
                                                                                                gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     >  2:1614916/1‑132 (MQ=255)
                                                                                                gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     <  1:1614916/132‑1 (MQ=255)
                                                                                                gCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTCCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGa                                                   >  2:1329926/1‑150 (MQ=255)
                                                                                                 ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTGGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                                  <  1:676900/150‑1 (MQ=255)
                                                                                                 ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                                  <  2:558888/150‑1 (MQ=255)
                                                                                                 ccGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                                  <  1:892852/150‑1 (MQ=255)
                                                                                                 ccGCCCTCCATACCGGAGCTTTTGAATCACAGGGGCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAg                                                  <  2:1587236/150‑1 (MQ=255)
                                                                                                  cGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAgg                                                 <  2:611498/150‑1 (MQ=255)
                                                                                                   gCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                                >  1:1058134/1‑150 (MQ=255)
                                                                                                    ctctCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGta                                               <  2:892599/150‑1 (MQ=255)
                                                                                                       tcCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAg                                            >  2:58786/1‑150 (MQ=255)
                                                                                                       tcCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAg                                            >  1:772062/1‑150 (MQ=255)
                                                                                                        ccATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGc                                           <  2:1039271/150‑1 (MQ=255)
                                                                                                         cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     <  1:98044/123‑1 (MQ=255)
                                                                                                         cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     >  2:98044/1‑123 (MQ=255)
                                                                                                         cATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGcc                                          >  1:705329/1‑150 (MQ=255)
                                                                                                          aTCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCg                                         >  1:1201272/1‑150 (MQ=255)
                                                                                                               ggCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         >  1:1010969/1‑97 (MQ=255)
                                                                                                               ggCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCg                                                                                         <  2:1010969/97‑1 (MQ=255)
                                                                                                                gcgcTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACg                                   <  2:481603/150‑1 (MQ=255)
                                                                                                                     ctTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACgatgat                              >  1:1538849/1‑150 (MQ=255)
                                                                                                                      ttGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATc                             >  2:322775/1‑150 (MQ=255)
                                                                                                                        gAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                     >  2:974593/1‑76 (MQ=255)
                                                                                                                        gAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAgg                                                                                                     <  1:974593/76‑1 (MQ=255)
                                                                                                                             acaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGtat                                              >  2:1017921/1‑126 (MQ=255)
                                                                                                                             acaGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGtat                                              <  1:1017921/126‑1 (MQ=255)
                                                                                                                                   ccACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGa                >  2:1348882/1‑150 (MQ=255)
                                                                                                                                      gatcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACAGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     <  2:1397975/90‑1 (MQ=255)
                                                                                                                                      gatcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACAGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCttt                                                                     >  1:1397975/5‑94 (MQ=255)
                                                                                                                                      cGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCa             <  1:356980/150‑1 (MQ=255)
                                                                                                                                       gTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAg            <  2:772062/150‑1 (MQ=255)
                                                                                                                                       atcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCAt                  >  2:718702/4‑144 (MQ=255)
                                                                                                                                       atcttCGGCGGCACCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCAt                  <  1:718702/141‑1 (MQ=255)
                                                                                                                                        tATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                                <  2:1168867/113‑1 (MQ=255)
                                                                                                                                        tATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGt                                                >  1:1168867/1‑113 (MQ=255)
                                                                                                                                                 gTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCa  <  1:322775/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TTACTCAGTGGTAATCCGCGAGTGCTGCTTGGTGTCTTTCATGGTCGCGTAAACCAGCAGCGAGCAGGCGATACAGGCGGTGACGTACCAGTAGAAGCCGCTCTCCATCCCGGCGCTCTTGAACCACAGGGCCACGTATTCGGCGGTGCCGCCGAAGATCGACACCGTCAGGGCATACGGCAGGCCAACGCCCAGGGCGCGAATTTCGGTAGGGAACAGCTCGGCTTTGACCACGGCGTTGATCGAGGTATAGCCGCTTACGATGATCAGTGCGGCCATGATCAGGAAGAACGCA  >  NC_002947/1594145‑1594439

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