Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,459 A→T 10.4% intergenic (+113/+29) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,4590AT10.4% 622.3 / 43.4 242intergenic (+113/+29)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (118/94);  new base T (18/5);  total (143/99)
Fisher's exact test for biased strand distribution p-value = 4.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 243 total aligned reads displayed.

GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4340313‑4340607
                                                                                                                                                  |                                                                                                                                                    
ggCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                   <  2:5941522/150‑1 (MQ=255)
    ggCAACGAACTTAACGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  2:9153276/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCcgc                                                                                                                                               <  2:4183763/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  1:2153100/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  2:2317699/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  1:8850692/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCcgc                                                                                                                                               >  2:1917293/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTACACCCGCGaa                                                                                                                                            >  2:2990215/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaa                                                                                                                                            >  2:5315509/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTCAACCCCCGaa                                                                                                                                            >  2:4099653/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCAGTCGCGGGGAAAACcgcgca                                                                                                                                            >  1:6877365/1‑148 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                           <  2:4982765/150‑1 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                           <  2:8507079/150‑1 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                           <  1:1880172/150‑1 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                           <  2:3770401/150‑1 (MQ=255)
         cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                          <  1:1204901/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                         >  1:6690729/1‑150 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                         >  2:5223083/1‑150 (MQ=255)
              ttAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCg                                                                                                                                     >  2:8355730/1‑150 (MQ=255)
                aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               >  2:2226116/1‑138 (MQ=255)
                aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTGGCGGGGAAAACCCCGACAGGGCCCGc                                                                                                                                   >  1:124398/1‑150 (MQ=255)
                  aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCaa                                                                                                                                 >  1:8981696/1‑150 (MQ=255)
                   gTGCATCTGCTGTGTGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGGGGGTAAACCCGCGAAAGGGCCGGCaaa                                                                                                                                <  2:81651/150‑1 (MQ=255)
                   gTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCcaa                                                                                                                                >  1:5185203/1‑150 (MQ=255)
                     gCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCCATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGGTGGCGGGTAAAACCGCGAAAGGGGCgggccaaa                                                                                                                              >  1:5340460/1‑144 (MQ=255)
                       aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTAAAACCCCGAAAGGGCCCGCaaaaaat                                                                                                                            >  1:6450321/1‑146 (MQ=255)
                         cTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGc                                                                                                                          >  2:2187448/1‑150 (MQ=255)
                            cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTAGCGGGTAAACCCGCGCAAGGGCCGGCAAATACTGCTCa                                                                                                                       >  2:768795/1‑150 (MQ=255)
                               tGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAggg                                                                                                                    >  2:3409819/1‑150 (MQ=255)
                                        aaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAccc                                                                                                           >  2:4443514/2‑150 (MQ=255)
                                           cAATAGCCCCGCCTGGGGGGGACTTTTGCTTCTATTTGCGGAAACCGCCCGCAGGCTACAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        <  1:6040696/150‑1 (MQ=255)
                                           cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATACAGCTCAGGGCTCAAACCCTTc                                                                                                        >  1:3642686/1‑150 (MQ=255)
                                           cAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGCAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        >  2:4440430/1‑150 (MQ=255)
                                            aaaaGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCTCAAACCCATcc                                                                                                       >  2:1069915/4‑150 (MQ=255)
                                            aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTAACCCGCGGAAGGGCCGGCACATCCTGCTCAGGGCGCAAACCCTTcc                                                                                                       >  1:8361961/1‑150 (MQ=255)
                                            aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAACCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       >  1:8191892/1‑150 (MQ=255)
                                            aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGAAACCCCCGAAAGGGCCGGCAAATCCGGCGCAGGGCGCAAACCCTTcc                                                                                                       >  2:6195042/1‑150 (MQ=255)
                                             aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                      >  2:6043039/1‑150 (MQ=255)
                                                 aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                  <  2:8066506/150‑1 (MQ=255)
                                                    cgcgTGGGGGGTTCTTTTGCTTCTTTTTGGGTAAACCGCCCGCTGGCTTCAATGAATCTGCCTGGGGTGCGGCTGGTCGGGCCGTTTCGGGGGTAACCCGGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt                                                                                               <  1:2659685/150‑1 (MQ=255)
                                                     gCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATc                                                                                              >  2:4309772/1‑150 (MQ=255)
                                                           ggggtttttttGCTTCTATTTGCGGAACCAGCCCGCCGGCTACAATTATTCTGCCAGGGGTGCGGATGGTCGGGCCCTTTCGGGTGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                        <  1:6268024/143‑1 (MQ=255)
                                                           ggggttattttGATTCTTTTTGCGTAACCCGCCCGCTGGCTTCAATGAATCTGCCAGCGGTGGGGCCGGTCCGGCCCTTTCGGGGGTAAACCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                        <  1:522421/143‑1 (MQ=255)
                                                           gggggTTTTTTGCTTCTTTTTGCGCAAACCGCCCGCAGGCTTATTTTAATCTGCCTGGGGGGCGGCTGGTCCGGCCCTTTGGCGGTTAACCCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                        <  1:8500251/148‑1 (MQ=255)
                                                            cGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTAGCGGGTAAACCCGCGAAAGGGCCGGCCAAACCAGCGCAGGGCTCAAACCCCTCCACAAAGGTCAACACAg                                                                                       >  2:415598/1‑150 (MQ=255)
                                                                  ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCACAGCCTTCg                                                                                 >  2:5162114/1‑150 (MQ=255)
                                                                   tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGGAAACCCCCGAAAGGGCCGGCCAAAACAGCTCAGGGCGCAAAACCATCCACAAAGATAACCGCAGccatagc                                                                                >  1:4470686/1‑145 (MQ=255)
                                                                      gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAc                                                                             >  2:5710333/1‑150 (MQ=255)
                                                                       cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCCGGTAAACCCGCGGAAGGGCCGGCAAAAACAGCGCAGGGCTCAAACCCATCCACAATGATACCCACAGCCTTAGGCAcc                                                                            >  2:3812188/1‑150 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTAGCGGGGGAACCCGCGAAAGGGGCCGCCAAAACAGCTCAGGGCGCAAAAACATACACAAAGATCAACGCAGCCCTGGCCAccc                                                                           >  1:3477499/1‑150 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCCGGGAAAAACCCCAAAGGGCCCGCAAAAAAAGCTAAGGGCACAAAACCCAACAAAAATATCACCACAGCCATCCCCAccc                                                                           >  1:7319101/1‑150 (MQ=255)
                                                                         tCTATTTGCGGAAACCGCCGGCTGGCTTCAATGAATCTGCCCGGGGGGCGCCAGGGCGGGCCCTTTCGCGGGTAACCCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAcccc                                                                          <  1:4208010/150‑1 (MQ=255)
                                                                                  gtaacccgCCCGCTGGCTTCAATGAAGCTGCCCGCGGTGCGGCTGGTCGGGCCCTTTCGCGGGAAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATg                                                                 <  1:1762858/143‑1 (MQ=255)
                                                                                  ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCGTCGCCACCCCTTGGCGATg                                                                 >  2:6232366/1‑150 (MQ=255)
                                                                                   gAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACCATGATAACCTCAGCCGTCGCCACCCCATGGCGATgg                                                                >  1:7617332/1‑150 (MQ=255)
                                                                                       cAGCCCGCCGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGc                                                            <  2:5090214/150‑1 (MQ=255)
                                                                                        aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCCGGGAAACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCGCAAAACCATACACAAAGATCACCACAGCCATCGCCACCCCCAGGCGCAGGCGGCa                                                           >  1:5926421/1‑150 (MQ=255)
                                                                                             gCAGGCTTCAAGGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCct                                                      <  1:5210255/150‑1 (MQ=255)
                                                                                                 gattCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGAAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  1:5500165/148‑1 (MQ=255)
                                                                                                 gCTTCAATTAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  1:2839451/150‑1 (MQ=255)
                                                                                                 gCTTCAATGACTCTGCCTGCGGTGCGGCTGGTCCGGCCATTTCGGGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  1:8751986/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  1:2050037/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  2:6946065/150‑1 (MQ=255)
                                                                                                  cTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                                 <  2:2822183/150‑1 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATTACCTCAGCCTTCGCCACCCCTTGGTGATGGCTGCATGCCGCTTGGt                                                >  2:6581964/1‑150 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGt                                                >  2:6744054/1‑150 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTACCGGGTAAACCCGCGTAAGGGCCGGCAAATCCTGCTCAGGGCGCAAAACCTTCCACAAAGATCACCTCAGCCGTCGCCACCCCTTGGCGATGGCAGCATGCCGCTTGGt                                                >  2:2102821/1‑150 (MQ=255)
                                                                                                      aaTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACa                                             <  2:1474004/150‑1 (MQ=255)
                                                                                                       aTGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCTTCCACAATGATCACCACACCCATCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAg                                            >  2:3507381/1‑150 (MQ=255)
                                                                                                            tCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGCAAGGGCCGGCCAATCCTGCTCAGGGCGCAAACCCGTCCACCATGATCACCTCAGCCGTCGCCCCCCCGTGGCGGTGGCTGCCGGCCGCGTGGTACAGCGCgg                                       >  2:4589143/1‑150 (MQ=255)
                                                                                                            tCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTTAACCCGCGGAAGGGCCGGCCAATACTGCGCAGGGGGCAAACCCGTACACAATGGTCACCTCAGCCGTCGCCACCCCGTGGCGGTGGCTGCATGCCGCGTGGTACAGCGCgg                                       >  2:1634241/1‑150 (MQ=255)
                                                                                                               gCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGc                                    <  2:1703344/150‑1 (MQ=255)
                                                                                                                 cAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                  <  2:3212301/150‑1 (MQ=255)
                                                                                                                  aGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                                 <  2:8302710/150‑1 (MQ=255)
                                                                                                                      gTGCGGCAGGTCCGGCCCTTTCGGGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                             <  1:3158040/150‑1 (MQ=255)
                                                                                                                         cGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGc                          <  2:7834947/150‑1 (MQ=255)
                                                                                                                           gCAGGTCGGGCCCTTTCGCGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGa                        <  1:2518630/150‑1 (MQ=255)
                                                                                                                           gCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGa                        >  2:6852386/1‑150 (MQ=255)
                                                                                                                                   ggCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAACCCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgc                <  2:8981696/150‑1 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAAACCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcg               <  1:4822402/150‑1 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcg               >  2:2037695/1‑150 (MQ=255)
                                                                                                                                     cccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc              >  2:2495040/1‑150 (MQ=255)
                                                                                                                                     cccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcgc              >  2:8388681/1‑150 (MQ=255)
                                                                                                                                      ccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg             >  1:1657499/1‑150 (MQ=255)
                                                                                                                                      ccTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGcgcg             >  1:7756916/1‑150 (MQ=255)
                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt            >  1:67907/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGtc           >  1:7808323/1‑149 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  1:8118221/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  2:2946195/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  2:5516483/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           <  2:1828242/150‑1 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  1:1106270/1‑150 (MQ=255)
                                                                                                                                         ttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCTGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc          >  1:2720390/1‑150 (MQ=255)
                                                                                                                                            gcgGTTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGaa       <  1:6169417/150‑1 (MQ=255)
                                                                                                                                             cgGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt      >  1:5870464/1‑150 (MQ=255)
                                                                                                                                              gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgcgc                                         <  2:8868848/114‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTAGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc     <  1:1917757/150‑1 (MQ=255)
                                                                                                                                              gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCATCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc     >  1:4220822/1‑150 (MQ=255)
                                                                                                                                                gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGa   >  2:7079373/1‑150 (MQ=255)
                                                                                                                                                 gaaaCCCGCGAAAGGGCCGGCAAATCCTGCTCCGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGaa  <  1:2990215/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4340313‑4340607

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