Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,628 A→G 7.4% intergenic (‑176/‑135) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,6280AG 7.4% 238.7 / 6.6 81intergenic (‑176/‑135)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/29);  new base G (2/4);  total (48/33)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCA  >  NC_002947/1639480‑1639775
                                                                                                                                                    |                                                                                                                                                   
gATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                                    >  2:1610955/1‑150 (MQ=255)
  ttCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACgg                                                                                                                                                  >  2:1326030/1‑150 (MQ=255)
       GACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                                 >  2:2856729‑M1/1‑146 (MQ=255)
       GACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                                 <  1:2856729‑M1/146‑1 (MQ=255)
       GACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTa                                                                                                                                             >  1:1810032‑M1/1‑150 (MQ=255)
          CCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAccc                                                                                                                                          >  1:1406698‑M1/1‑150 (MQ=255)
           CTGTAGTAACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGTTCCTACAGGCGCGGGTttactcg                                                                                                                                         >  1:1472644‑M1/1‑147 (MQ=255)
               aGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGaa                                                                                                                                     <  1:1367163/150‑1 (MQ=25)
               aGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                     <  1:1326030/150‑1 (MQ=255)
                 tGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAaga                                                                                                                                   >  1:2046004/1‑150 (MQ=255)
                   aCATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAgg                                                                                                                                 <  2:1088657/150‑1 (MQ=25)
                      tCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGccc                                                                                                                              <  2:44544/150‑1 (MQ=255)
                            tcaagCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACg                                                                                                                                                   <  1:1231935‑M2/118‑1 (MQ=255)
                            tcaagCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACg                                                                                                                                                   >  2:1231935‑M2/6‑123 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAgg                                                                                                                                 >  1:1961149/1‑133 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAgg                                                                                                                                 <  2:1961149/133‑1 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                                >  1:1444945/1‑134 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                                <  2:1444945/134‑1 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                      >  1:222375/1‑144 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                                      <  2:222375/144‑1 (MQ=255)
                                    cccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacaca                                                                                                                >  1:2246583/1‑150 (MQ=255)
                                                gCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGc                                                                                                    >  2:2172304/1‑150 (MQ=255)
                                                  atcaaaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGTTGTCTGTTCGGGCCCCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAg                                                                                                  <  1:780837/147‑1 (MQ=255)
                                                  aTGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTAGCCGGGAAGAGGCCCGAACAGTCGACACAAGGGCACATTGCAg                                                                                                  >  2:780837/1‑150 (MQ=255)
                                                     aaaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGAc                                                                                               <  2:2046004/150‑1 (MQ=255)
                                                       aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcg                                                                                                                                         >  1:688427/1‑106 (MQ=255)
                                                       aaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcg                                                                                                                                         <  2:688427/106‑1 (MQ=255)
                                                        aaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGt                                                                                            <  1:1610955/150‑1 (MQ=255)
                                                              cTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCGTACAGGCGCGGGTTTACTCGCGAagag                                                                                                                                  <  1:2093276/106‑1 (MQ=17)
                                                              cTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAagag                                                                                                                                  >  2:2093276/1‑106 (MQ=14)
                                                                        gATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCt                                                                            >  2:927896/1‑150 (MQ=255)
                                                                            aCGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                                        <  2:2246583/150‑1 (MQ=255)
                                                                                aTTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                     >  1:892486/1‑85 (MQ=255)
                                                                                aTTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGaa                                                                                                                                     <  2:892486/85‑1 (MQ=255)
                                                                                              tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                                      >  1:2192416/1‑150 (MQ=255)
                                                                                              tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTa                                                      >  1:1360557/1‑150 (MQ=255)
                                                                                                   tCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATg                                                 >  1:7254/1‑150 (MQ=255)
                                                                                                   tCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATg                                                 <  2:7254/150‑1 (MQ=255)
                                                                                                    cTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                                <  2:41728/70‑1 (MQ=255)
                                                                                                    cTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGc                                                                                                                                >  1:41728/1‑70 (MQ=255)
                                                                                                          gtctgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                                          >  1:750388/1‑150 (MQ=255)
                                                                                                           tctgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTc                                         >  1:1657453/1‑150 (MQ=255)
                                                                                                             tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgagatcg                                       >  1:2969283/1‑144 (MQ=255)
                                                                                                              gtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCATCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                                      >  1:2507282/1‑150 (MQ=255)
                                                                                                               ttCGGGCCTCTTGCCGAGTAAACGCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                                     <  2:2478724/150‑1 (MQ=255)
                                                                                                               ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAc                                                                              >  1:934599/1‑109 (MQ=255)
                                                                                                               ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAc                                                                              <  2:934599/109‑1 (MQ=255)
                                                                                                               ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                                     >  1:829016/1‑150 (MQ=255)
                                                                                                               ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAt                                     >  1:2587103/1‑150 (MQ=255)
                                                                                                                tCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATc                                    <  1:927896/150‑1 (MQ=255)
                                                                                                                     ccTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGtta                               <  2:2278843/150‑1 (MQ=255)
                                                                                                                      ctctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttat                              >  2:981086/1‑150 (MQ=255)
                                                                                                                       tctTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatt                             >  2:1785740/1‑150 (MQ=255)
                                                                                                                         ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgt                                            >  2:156387/1‑133 (MQ=255)
                                                                                                                         ttCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgt                                            <  1:156387/133‑1 (MQ=255)
                                                                                                                          tCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACg                          <  1:2815050/150‑1 (MQ=255)
                                                                                                                           cgcgAGTAAACACGCGCCAGCAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa                         <  2:2192416/150‑1 (MQ=255)
                                                                                                                             cgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAgg                       <  2:1360557/150‑1 (MQ=255)
                                                                                                                              gAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc                      >  2:1957461/1‑150 (MQ=255)
                                                                                                                                gTAAACCCGTTCCTACAGGCGCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCgg                                  <  2:2952042/136‑1 (MQ=37)
                                                                                                                                gTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCAGTGTCTGGTTGGTTATTACGAGGCaa                    >  2:1476028/1‑150 (MQ=255)
                                                                                                                                 tAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCaaa                   >  1:2964347/1‑150 (MQ=255)
                                                                                                                                   aaCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGc                                               <  1:1547751/120‑1 (MQ=255)
                                                                                                                                   aaCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGc                                               >  2:1547751/1‑120 (MQ=255)
                                                                                                                                    aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac                >  2:2244375/1‑150 (MQ=255)
                                                                                                                                    aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac                >  1:277427/1‑150 (MQ=255)
                                                                                                                                    aCCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac                >  2:717359/1‑150 (MQ=255)
                                                                                                                                       cgcgCCTGTAGGAACGGGTTTACCCGCGGAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGAt             <  1:2172304/150‑1 (MQ=255)
                                                                                                                                        gcgcCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                              >  2:2772057/1‑52 (MQ=255)
                                                                                                                                        gcgcCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAg                                                                                                              <  1:2772057/52‑1 (MQ=255)
                                                                                                                                        gcgcCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATccc                                                                                      >  1:2810572/1‑76 (MQ=255)
                                                                                                                                        gcgcCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATccc                                                                                      <  2:2810572/76‑1 (MQ=255)
                                                                                                                                            cTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgt                                            >  2:1275080/1‑114 (MQ=255)
                                                                                                                                            cTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgt                                            <  1:1275080/114‑1 (MQ=255)
                                                                                                                                              gTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTg      >  2:775217/1‑150 (MQ=255)
                                                                                                                                              gTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTg      >  1:892306/1‑150 (MQ=255)
                                                                                                                                              gTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTg      >  1:1282207/1‑150 (MQ=255)
                                                                                                                                               tAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGa     <  1:1785740/150‑1 (MQ=255)
                                                                                                                                                  gAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCCAACAAGGACACATTGCAGGAGTGTGATGCCCGGATCGGCTTGACGTGAGTGAACCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAAGACGATGAAGCTGACCa  >  1:1238897/1‑150 (MQ=255)
                                                                                                                                                  gAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCa  >  1:1387748/1‑150 (MQ=255)
                                                                                                                                                  gAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCCGATCGGTTATTACGAGGCAAACACGATGAACCTGACCa  >  1:2941211/1‑150 (MQ=255)
                                                                                                                                                   aaCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                                        >  2:270351/1‑79 (MQ=255)
                                                                                                                                                   aaCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                                        <  1:270351/79‑1 (MQ=255)
                                                                                                                                                   aaCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc                           >  2:1948100/1‑124 (MQ=255)
                                                                                                                                                   aaCGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc                           <  1:1948100/124‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGATGAACCTGACCA  >  NC_002947/1639480‑1639775

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: