Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,307,391 G→A 5.0% intergenic (‑92/+94) recC ← / ← PP_5SG multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,307,3910GA 5.0% 361.8 / 6.0 119intergenic (‑92/+94)recC/PP_5SGmultifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease recombinase subunit/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (60/52);  new base A (2/4);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 119 total aligned reads displayed.

GCGGTTGCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGA  >  NC_002947/5307246‑5307538
                                                                                                                                                 |                                                                                                                                                   
gCGGTTGCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGcc                                                                                                                                                 >  1:1172800/1‑150 (MQ=255)
 cGGTTGCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGccc                                                                                                                                                >  2:423568/1‑150 (MQ=255)
       ccGCGGCCGATCATGAAGCCAGGGTGGGGATTGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGGCACGGCCGCTccc                                                                                                                                          <  1:1411629/150‑1 (MQ=255)
          tGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGGCACGGCCGCTCCCGCa                                                                                                                                       <  1:128948/150‑1 (MQ=255)
               tctCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGTTCTGCCGTGCCCGCTCCCGCAGCGAt                                                                                                                                  <  1:2718610/148‑1 (MQ=255)
                aTCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTc                                                                                                                                            >  2:37199/1‑139 (MQ=255)
                aTCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTc                                                                                                                                            <  1:37199/139‑1 (MQ=255)
                 tCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTc                                                                                                                                >  1:814604/1‑150 (MQ=255)
                   aTGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTcgc                                                                                                                              >  1:952521/1‑150 (MQ=255)
                   aTGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACACGACAACTCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGAttaga                                                                                                                              >  2:2718610/1‑147 (MQ=255)
                        gCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGGCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAgg                                                                                                                         >  2:83832/1‑150 (MQ=255)
                         ccAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCAGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGt                                                                                                                        >  1:829736/1‑150 (MQ=255)
                                 gagaTGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGtttt                                                                                                                >  2:1117523/1‑150 (MQ=255)
                                 gagaTGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGtttt                                                                                                                >  1:983984/1‑150 (MQ=255)
                                      gggCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGcc                                                                                                           <  2:1647028/150‑1 (MQ=255)
                                        gCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGc                                                                                                                                              >  1:2406665/1‑113 (MQ=255)
                                        gCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCca                                                                                                         <  2:1141432/150‑1 (MQ=255)
                                             gtTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGGCACGGCCGCTCCCGCa                                                                                                                                       <  1:1499600/115‑1 (MQ=39)
                                             gtTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGGCACGGCCGCTCCCGCa                                                                                                                                       >  2:1499600/1‑115 (MQ=39)
                                             gtTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAg                                                                                                    >  1:2581630/1‑150 (MQ=255)
                                              ttGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGa                                                                                                   <  2:829736/150‑1 (MQ=255)
                                               tGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGAGATCGGACTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGCTTCGCGCAGGCCGAGTTGTTTGCCCACAAAGAc                                                                                                  >  1:2065230/1‑150 (MQ=255)
                                                   cATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCa                                                                                                                                       <  1:453781/109‑1 (MQ=255)
                                                   cATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCa                                                                                                                                       >  2:453781/1‑109 (MQ=255)
                                                    aTCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGc                                                                                                            >  2:928522/1‑135 (MQ=255)
                                                    aTCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGc                                                                                                            <  1:928522/135‑1 (MQ=255)
                                                     tCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTATACTAGACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAAcc                                                                                            <  1:809906/150‑1 (MQ=255)
                                                     tCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAAcc                                                                                            >  1:629409/1‑150 (MQ=255)
                                                        gAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGtttttt                                                                                                              <  1:2784955/129‑1 (MQ=255)
                                                            gTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCt                                                                                     >  1:945012/1‑150 (MQ=255)
                                                               tctGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcg                                                                                  >  1:1210000/1‑150 (MQ=255)
                                                               tctGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcg                                                                                  >  2:2426183/1‑150 (MQ=255)
                                                                 tGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCTGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAAcccc                                                                                          <  1:2045752/140‑1 (MQ=255)
                                                                 tGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCTGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAAcccc                                                                                          >  2:2045752/1‑140 (MQ=255)
                                                                       cgcTGATGGCTGAGAGGGGACCGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCGGGCTCCCGCAGAGATTTCAATACTCCTAGCGGCTGTCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGa                                                                          <  1:1302850/150‑1 (MQ=255)
                                                                        gcTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTCCCTGCGTGCGCAGAt                                                                         >  2:2949657/1‑150 (MQ=255)
                                                                        gcTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCAGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTGTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCCGAt                                                                         <  1:1012091/150‑1 (MQ=255)
                                                                          tGATGGCTGGGAGTGTAACTATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAg                                                                       <  2:952521/150‑1 (MQ=255)
                                                                          tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGc                                                                                                                                                  <  1:2312900/75‑1 (MQ=255)
                                                                          tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGc                                                                                                                                                  >  2:2312900/1‑75 (MQ=255)
                                                                          tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAg                                                                       <  2:1210000/150‑1 (MQ=255)
                                                                          tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAg                                                                       <  1:83832/150‑1 (MQ=255)
                                                                          tGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAg                                                                       >  1:1049569/1‑150 (MQ=255)
                                                                            aTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                     >  1:65923/1‑150 (MQ=255)
                                                                            aTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAggg                                                                     <  2:2168811/150‑1 (MQ=255)
                                                                             tGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAgggg                                                                    <  1:1117523/150‑1 (MQ=255)
                                                                                 tGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAGGGGtttt                                                                <  2:15948/150‑1 (MQ=39)
                                                                                 tGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGtttt                                                                >  1:2146222/1‑150 (MQ=255)
                                                                                   agagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGc                                                              <  1:1682755/150‑1 (MQ=255)
                                                                                      ggggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAGGGGtttt                                                                <  2:1968796/145‑1 (MQ=39)
                                                                                      ggggAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTATCTGCGCACGCAGGTAGGGGtttt                                                                >  1:1968796/1‑145 (MQ=255)
                                                                                           aCGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGc                                                                                                            <  2:1508200/96‑1 (MQ=255)
                                                                                             gATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATc                                                    <  1:302060/150‑1 (MQ=255)
                                                                                              atatCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                          <  2:249110/143‑1 (MQ=255)
                                                                                               tatCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATATGGGTTTTGCGAAATGAATCtt                                                  <  2:673084/150‑1 (MQ=255)
                                                                                               tatCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGGATCtt                                                  <  2:2146222/150‑1 (MQ=255)
                                                                                                 tCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                >  2:1606541/1‑150 (MQ=255)
                                                                                                 tCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGa                                                <  2:1049569/150‑1 (MQ=255)
                                                                                                  cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGAc                                               <  1:2426183/150‑1 (MQ=255)
                                                                                                  cACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGAc                                               <  2:629409/150‑1 (MQ=255)
                                                                                                   aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGttt                                                                 >  1:21196/1‑131 (MQ=255)
                                                                                                   aCGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGttt                                                                 <  2:21196/131‑1 (MQ=255)
                                                                                                     gcgTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCa                                                                                                                           <  2:2823493/71‑1 (MQ=255)
                                                                                                     gcgTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCa                                                                                                                           >  1:2823493/1‑71 (MQ=255)
                                                                                                        tACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAg                                                                                                                                      <  2:39064/57‑1 (MQ=255)
                                                                                                        tACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAg                                                                                                                                      >  1:39064/1‑57 (MQ=255)
                                                                                                        tACACGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGAc                                         >  2:2397379/1‑150 (MQ=255)
                                                                                                           aTGAAACGACCTAATTCCATCATGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTa                                      >  2:2348654/1‑150 (MQ=255)
                                                                                                            tGAAACGGCCGAATGCCATCAGGCTGCGGGAACGGCCTTGCCCATTCCCGCAGCGATTCGCGCAGGGCCATGGTTTTGCACCCCAATACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAACCGAAACAGTACGATTACCTAc                                     >  1:218682/1‑150 (MQ=16)
                                                                                                             gAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTAct                                    >  2:2639950/1‑150 (MQ=255)
                                                                                                               aaCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                          <  2:2923722/126‑1 (MQ=255)
                                                                                                               aaCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGaaa                                                          >  1:2923722/1‑126 (MQ=255)
                                                                                                                aCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGtttttt                                                                                                              >  2:2939471/1‑73 (MQ=255)
                                                                                                                aCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGtttttt                                                                                                              <  1:2939471/73‑1 (MQ=255)
                                                                                                                aCGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACtctc                                 >  1:897779/1‑150 (MQ=255)
                                                                                                                  ggCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGttttt                                                                                                               <  2:973635/70‑1 (MQ=255)
                                                                                                                  ggCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAACAACAACACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTcac                               >  1:53953/1‑150 (MQ=255)
                                                                                                                   gCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGATAATAATCTTGACGATGACCTACTCTcaca                              <  1:1606541/150‑1 (MQ=255)
                                                                                                                     cGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTc                                                                                                                                            <  1:2995985/38‑1 (MQ=255)
                                                                                                                     cGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTc                                                                                                                                            >  2:2995985/1‑38 (MQ=255)
                                                                                                                      gAATGCCATTAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGTGTTCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                           <  2:753506/150‑1 (MQ=255)
                                                                                                                      gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGtt                                                                  >  2:2313970/1‑111 (MQ=255)
                                                                                                                      gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGtt                                                                  <  1:2313970/111‑1 (MQ=255)
                                                                                                                      gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                           >  1:1983660/1‑150 (MQ=255)
                                                                                                                      gAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                           >  1:1884961/1‑150 (MQ=255)
                                                                                                                       aaTGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATggg                          <  2:2581630/150‑1 (MQ=255)
                                                                                                                         tGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGa                        >  1:123120/1‑150 (MQ=255)
                                                                                                                          gCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATTTTGGCAATTACCTACTCTCACATGGGGaa                       >  1:778405/1‑150 (MQ=255)
                                                                                                                           ccATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCtt                                                  <  2:182155/122‑1 (MQ=255)
                                                                                                                           ccATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCtt                                                  >  1:182155/1‑122 (MQ=255)
                                                                                                                           ccATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGaaa                      >  2:879499/1‑150 (MQ=255)
                                                                                                                               cAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAcccc                  <  1:1800318/150‑1 (MQ=255)
                                                                                                                                  gCTGCGGGAGCGGGCACGGCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAgg                                                                      >  2:446304/1‑95 (MQ=21)
                                                                                                                                  gCTGCGGGAGCGGGCACGGCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAgg                                                                      <  1:446304/95‑1 (MQ=21)
                                                                                                                                       gggAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTAcc          >  1:651908/1‑150 (MQ=255)
                                                                                                                                             ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcgtgcg                                                                              >  2:1102835/1‑76 (MQ=255)
                                                                                                                                             ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCtgcgtgcg                                                                              <  1:1102835/76‑1 (MQ=255)
                                                                                                                                              gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCGCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGa                                                       >  1:938699/1‑98 (MQ=255)
                                                                                                                                              gCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCGCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGa                                                       <  2:938699/98‑1 (MQ=255)
                                                                                                                                               ccGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGa  <  2:2933523/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCGGTTGCCGTGGACGATCATGAAGCCAGGGTGGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGGCGA  >  NC_002947/5307246‑5307538

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: