Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,094,547 G→C 7.2% intergenic (‑273/+64) astC ← / ← PP_4482 acetylornithine aminotransferase/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,094,5470GC 7.2% 451.1 / 19.5 138intergenic (‑273/+64)astC/PP_4482acetylornithine aminotransferase/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (70/58);  new base C (7/3);  total (77/61)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 155 total aligned reads displayed.

CAAGGCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAA  >  NC_002947/5094398‑5094696
                                                                                                                                                     |                                                                                                                                                     
cAAGGCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTgg                                                                                                                                                       >  1:1825094/1‑150 (MQ=255)
 aaGGCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTggg                                                                                                                                                      >  1:939478/1‑150 (MQ=255)
   ggCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcg                                                                                                                                                    <  1:1270850/150‑1 (MQ=255)
    gCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTgg                                                                                                                                                       <  1:1913729/146‑1 (MQ=255)
    gCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcgg                                                                                                                                                   <  1:3547797/150‑1 (MQ=255)
     cTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcggc                                                                                                                                                  <  1:464902/150‑1 (MQ=255)
            cAATGCGACAAACCCTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGACCGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGAt                                                                                                                                           <  1:517254/150‑1 (MQ=255)
             aaTGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGAtt                                                                                                                                          <  2:438412/150‑1 (MQ=255)
              aTGCGACAAACACTGTCGCAATCGCTCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGAttt                                                                                                                                         <  2:2358410/150‑1 (MQ=255)
                gCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGCGGGGGCGGGGTGAACGGGTGCTCACACAGGATTGGGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTgc                                                                                                                                       <  1:227995/150‑1 (MQ=255)
                  gACAAACCCTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTgcgc                                                                                                                                     <  2:2972047/150‑1 (MQ=255)
                  gACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTgcgc                                                                                                                                     <  2:2493014/150‑1 (MQ=255)
                  gACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTgcgc                                                                                                                                     <  2:2972053/150‑1 (MQ=255)
                     aaaCACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGcccc                                                                                                                                  >  2:4592386/1‑150 (MQ=255)
                      aaCACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCa                                                                                                                                 >  1:2334027/1‑150 (MQ=255)
                       acacTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGAttt                                                                                                                                         <  2:708143/141‑1 (MQ=255)
                           tGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTggg                                                                                                                                                      >  1:3471005/1‑124 (MQ=255)
                           tGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcgg                                                                                                                                                   <  2:2304120/127‑1 (MQ=255)
                           tGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcgg                                                                                                                                                   >  1:2304120/1‑127 (MQ=255)
                           tGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGggcgg                                                                                                                                                   <  2:2303785/127‑1 (MQ=255)
                             tCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCaa                                                                                                                          >  1:4337924/1‑150 (MQ=255)
                             tCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCaa                                                                                                                          >  1:397533/1‑150 (MQ=255)
                                   tCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCaa                                                                                                                    >  2:1958869/1‑150 (MQ=255)
                                   tCGCGCAGTGTACGGGAGAGAGGGGGCGGGGTGACCGGGTGCTCCCACAGGAGTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCaa                                                                                                                    <  1:814499/150‑1 (MQ=255)
                                     gcgcAGTGTACGGGAGAGCGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAgg                                                                                                                  <  1:4745035/150‑1 (MQ=255)
                                     gcgcAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAaca                                                                                                                              <  2:3575543/138‑1 (MQ=255)
                                     gcgcAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAaca                                                                                                                              >  1:3575543/1‑138 (MQ=255)
                                      cgcAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGc                                                                                                                 <  2:5429926/150‑1 (MQ=255)
                                         aGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATa                                                                                                              <  2:2374964/150‑1 (MQ=255)
                                         aGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATa                                                                                                              <  2:3520631/150‑1 (MQ=255)
                                         aGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATa                                                                                                              <  2:3409303/150‑1 (MQ=255)
                                           tgtACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATaca                                                                                                            <  1:3624465/150‑1 (MQ=255)
                                             tACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACAcc                                                                                                          >  2:1837376/1‑150 (MQ=255)
                                               cGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACAcccc                                                                                                        <  1:3058586/150‑1 (MQ=255)
                                                  gagagaGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGaa                                                                                                     >  1:4158489/1‑150 (MQ=255)
                                                        gggggCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAacccacc                                                                                               >  2:4942881/1‑150 (MQ=255)
                                                            gCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGaa                                                                                           >  2:1722007/1‑150 (MQ=255)
                                                                ggTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtct                                                                                       <  1:1417556/150‑1 (MQ=255)
                                                                ggTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtct                                                                                       >  1:5364375/1‑150 (MQ=255)
                                                                ggTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtct                                                                                       >  1:5396759/1‑150 (MQ=255)
                                                                 gTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAAtctc                                                                                      >  2:2263159/1‑150 (MQ=255)
                                                                  tGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGCCGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     <  1:554104/150‑1 (MQ=255)
                                                                  tGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     <  2:559386/150‑1 (MQ=255)
                                                                  tGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     <  1:554146/150‑1 (MQ=255)
                                                                           gCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTc                                                                            >  2:2786867/1‑150 (MQ=255)
                                                                             tCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAg                                                                          >  2:1645780/1‑150 (MQ=255)
                                                                             tCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAg                                                                          >  1:3739129/1‑150 (MQ=255)
                                                                               acacaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGcc                                                                        >  1:2422072/1‑150 (MQ=255)
                                                                                 acaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCgg                                                                      >  2:2705689/1‑150 (MQ=255)
                                                                                 acaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCgg                                                                      >  2:1425814/1‑150 (MQ=255)
                                                                                 acaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCgg                                                                      >  2:1425807/1‑150 (MQ=255)
                                                                                 acaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCgg                                                                      >  1:5240358/1‑150 (MQ=255)
                                                                                 acaGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCgg                                                                      >  2:3238142/1‑150 (MQ=255)
                                                                                   aGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGgcg                                                                    >  1:4235485/1‑150 (MQ=255)
                                                                                      aTTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCg                                                                 <  1:4306515/150‑1 (MQ=255)
                                                                                       ttGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGa                                                                <  2:3773056/150‑1 (MQ=255)
                                                                                       ttGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGa                                                                <  2:3482401/150‑1 (MQ=255)
                                                                                          tgAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCt                                                             <  2:5307696/150‑1 (MQ=255)
                                                                                          tgAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCt                                                             >  1:2187787/1‑150 (MQ=255)
                                                                                          tgAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCt                                                             >  2:4632640/1‑150 (MQ=255)
                                                                                           gAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTc                                                            >  1:2518366/1‑150 (MQ=255)
                                                                                                ggCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTc                                                       >  1:2265847/1‑150 (MQ=255)
                                                                                                   aGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     <  2:3200516/117‑1 (MQ=255)
                                                                                                   aGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     >  1:3200516/1‑117 (MQ=255)
                                                                                                   aGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCaaa                                                    >  1:457356/1‑150 (MQ=255)
                                                                                                   aGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTGTTCAGCCGGCGCCGAGCTCAATTCaaa                                                    >  2:5086190/1‑150 (MQ=255)
                                                                                                    ggATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAAc                                                   <  1:1070329/150‑1 (MQ=255)
                                                                                                            ggCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgct                                           >  2:4866071/1‑150 (MQ=255)
                                                                                                            ggCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGgctgct                                           >  1:3280093/1‑150 (MQ=255)
                                                                                                               gATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     >  1:4716131/1‑105 (MQ=255)
                                                                                                               gATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCa                                                                                     <  2:4716131/105‑1 (MQ=255)
                                                                                                                aTGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGctgctgctg                                       >  1:1861373/1‑150 (MQ=255)
                                                                                                                  gTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCg                                     >  2:2206681/1‑150 (MQ=255)
                                                                                                                  gTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCg                                     >  2:4285144/1‑150 (MQ=255)
                                                                                                                      aTGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGc                                                              >  2:5411278/1‑121 (MQ=255)
                                                                                                                         tgtTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGa                              <  2:939478/150‑1 (MQ=255)
                                                                                                                         tgtTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGa                              >  2:1120216/1‑150 (MQ=255)
                                                                                                                               ggcgcaaaTCGAGCGCCGCCCACGCGGCGCTCGATCTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCgg                        >  2:156912/8‑150 (MQ=38)
                                                                                                                               cGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCgg                        <  1:4632640/150‑1 (MQ=255)
                                                                                                                                ggggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGa                                    <  2:2599966/137‑1 (MQ=255)
                                                                                                                                ggggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGa                                    >  1:2599966/1‑137 (MQ=255)
                                                                                                                                  ggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAAtt                                                        >  1:3906900/1‑115 (MQ=255)
                                                                                                                                  ggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAAtt                                                        <  2:3906900/115‑1 (MQ=255)
                                                                                                                                  ggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCt                     >  1:3025434/1‑150 (MQ=255)
                                                                                                                                  ggAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCt                     >  1:1573964/1‑150 (MQ=255)
                                                                                                                                    agaTCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTc                   <  1:4592386/150‑1 (MQ=255)
                                                                                                                                     aaTCGAGCGCCGCCCACGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCt                                                             <  1:2464586/106‑1 (MQ=39)
                                                                                                                                      aTCGCGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTccc                 <  2:4420093/150‑1 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCACGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg                >  1:3440220/1‑150 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCACGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg                >  2:4055735/1‑150 (MQ=255)
                                                                                                                                       tCGAGCGCCGCCCACGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcg                >  1:1019624/1‑150 (MQ=255)
                                                                                                                                        cGAGCGCCGCCCACGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcgc               >  1:4646823/1‑150 (MQ=255)
                                                                                                                                         gAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCcgcg              >  2:2469016/1‑150 (MQ=255)
                                                                                                                                             gccgcGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGt          >  1:1833729/1‑150 (MQ=255)
                                                                                                                                               cgcgTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTg        <  2:459668/150‑1 (MQ=255)
                                                                                                                                               cgcgTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTg        <  1:2746276/150‑1 (MQ=255)
                                                                                                                                                gcgTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGc       <  1:3597365/150‑1 (MQ=255)
                                                                                                                                                 cgTGGGCGGCGCTCGATGTGCGCCCCAACACCAACCCCAAGGCAGACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCa      <  2:3325448/150‑1 (MQ=255)
                                                                                                                                                   tGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCAcc    <  1:1157463/150‑1 (MQ=255)
                                                                                                                                                     ggcggcGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCaa  >  2:683094/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAAGGCTTTCGGCAATGCGACAAACACTGTCGCAATCGCGCAGTGTACGGGAGAGAGGGGGCTGGGTGAACGGGTGCTCACACAGGATTGTGAGGGGGCAGGATGTACGGCGATGTTGATGTGTTGACGGGGAGATCGAGCGCCGCGTGGGCGGCGCTCGATTTGCGCCCCAACACCAACCCCAAGGCATACACCCCGAACCCACCCGAAATCTCAGCCCTTTTCAGCCGGCGCCGAGCTCAATTCAAACGGGCTGCTGCTGCGACGCTGATTGCGGTCTTCCCGCGGCGTTGCACCAA  >  NC_002947/5094398‑5094696

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