Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,104,984 A→C 5.5% intergenic (+73/+80) PP_4488 → / ← phhR hypothetical protein/sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,104,9840AC 5.5% 316.0 / 3.0 114intergenic (+73/+80)PP_4488/phhRhypothetical protein/sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (33/72);  new base C (6/2);  total (39/75)
Fisher's exact test for biased strand distribution p-value = 1.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

TGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCA  >  NC_002947/5104835‑5105131
                                                                                                                                                     |                                                                                                                                                   
tGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCa                                                                                                                                                     <  2:4016685/150‑1 (MQ=255)
tGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCa                                                                                                                                                     <  2:2972107/150‑1 (MQ=255)
 gACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCag                                                                                                                                                    <  1:2840442/150‑1 (MQ=255)
  aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                   <  2:416645/150‑1 (MQ=255)
  aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                   <  1:5435405/150‑1 (MQ=255)
  aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                   <  2:624252/150‑1 (MQ=255)
  aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                   <  2:3364574/150‑1 (MQ=255)
  aCCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCaga                                                                                                                                                   <  1:5291324/150‑1 (MQ=255)
   ccTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCagag                                                                                                                                                  <  1:5370921/150‑1 (MQ=255)
    cTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACACATAATGCGCGCCGCCAGAGc                                                                                                                                                 >  2:1143239/1‑150 (MQ=255)
    cTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGGCACAAGGCCGCTCCCCCAGAAAAAGCGCTCAGCCAgcgc                                                                                                                                                 >  2:4358509/1‑150 (MQ=255)
     tACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGca                                                                                                                                                <  2:3285245/150‑1 (MQ=255)
     tACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGca                                                                                                                                                <  2:3285491/150‑1 (MQ=255)
       cGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGcaca                                                                                                                                              <  1:322853/150‑1 (MQ=255)
         aGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCCCAGATAAAGCGCTCAGCCCGAGCACAGc                                                                                                                                            >  2:2093444/1‑150 (MQ=255)
         aGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGc                                                                                                                                            <  2:1619075/150‑1 (MQ=255)
         aGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCAATCGCGACACAAGGCCGCTCCCACAGATAAAGCGCTCAGCCAGAGcacacc                                                                                                                                            >  2:2093512/1‑148 (MQ=255)
           gACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGt                                                                                                                                          <  1:5591254/150‑1 (MQ=255)
            aCTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTa                                                                                                                                         <  2:4080308/150‑1 (MQ=255)
             cTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTAt                                                                                                                                        <  1:3688252/150‑1 (MQ=255)
              tCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGGCCTGCGGGGCCCTATCGCGAAACAATGCCGCTACCACAGATAATTCGCTAAGGCAGAGCACAGCGTAtt                                                                                                                                       >  2:2563378/1‑150 (MQ=255)
              tCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGGCCTGCGGGCCCCTATCGCGTGAAAATGCCGCGCCCACAGATTATTGGCTCAGGCAGAGCACAGGGTAtt                                                                                                                                       >  2:2563815/1‑150 (MQ=255)
                ccAGGGCCAGAAGCACAGTCTCGAAGACCGCGTGAGGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTAc                                                                                                                                     <  2:5572330/150‑1 (MQ=255)
                  aGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTTATTCCTTTGAATCGGGGGCCCTGCGGGCCCCTATCGCGACCCAAGGCCGCTCCCACAGATTATGCGCTCAGCCAGAGCTCAGCGTATTACCt                                                                                                                                   <  1:2367553/150‑1 (MQ=255)
                  aGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCt                                                                                                                                   <  1:804985/150‑1 (MQ=255)
                  aGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCt                                                                                                                                   <  1:4704761/150‑1 (MQ=255)
                  aGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCt                                                                                                                                   >  2:779198/1‑150 (MQ=255)
                    ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                                                                                 <  1:4411745/150‑1 (MQ=255)
                    ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                                                                                 <  2:5003996/150‑1 (MQ=255)
                    ggCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgt                                                                                                                                 <  1:4411741/150‑1 (MQ=255)
                     gCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTCGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCtgtg                                                                                                                                <  1:3387696/150‑1 (MQ=255)
                      ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACCCAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                                                                               <  1:459183/150‑1 (MQ=255)
                      ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                                                                               <  2:2165448/150‑1 (MQ=255)
                      ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                                                                               <  2:570837/150‑1 (MQ=255)
                      ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                                                                               <  1:2628282/150‑1 (MQ=255)
                      ccAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTgg                                                                                                                               <  1:3221285/150‑1 (MQ=255)
                         gAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGGTAATACGCTGTGCTCTGGCTGAGCGCATTATCTGTGGGAg                                                                                                                            <  1:2602640/150‑1 (MQ=21)
                         gAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAg                                                                                                                            <  1:3383122/150‑1 (MQ=255)
                          aaGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGc                                                                                                                           <  1:3794102/150‑1 (MQ=255)
                          aaGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGc                                                                                                                           <  1:1020633/150‑1 (MQ=255)
                            gCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCgg                                                                                                                         <  1:4135012/150‑1 (MQ=255)
                              acaTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGAAACAAGGCCGCTCCCACCGATAATGCGCACAGCCAGAGCACAGCGTAGTACCAGTGGGAGCGGcc                                                                                                                       >  2:3427302/1‑150 (MQ=255)
                                aTTCTCGAATACCGCGTGATGGGGAACGGCTTCTCCATAGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCtt                                                                                                                     <  2:3158832/150‑1 (MQ=255)
                                aTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCtt                                                                                                                     <  1:1640889/150‑1 (MQ=255)
                                aTTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCtt                                                                                                                     <  1:4432952/150‑1 (MQ=255)
                                  tctcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACCCGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTtgt                                                                                                                   <  1:1835388/150‑1 (MQ=255)
                                   ctcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCTTTTTACCTGTGGGAGCGGCCTtgtg                                                                                                                  <  1:1677322/150‑1 (MQ=255)
                                    tcGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCAGTGGGAGCGGCCTTgtgt                                                                                                                 >  2:5511625/1‑150 (MQ=255)
                                       caTACCGCGTGCTGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACCAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTcgc                                                                                                              <  2:3275042/149‑1 (MQ=255)
                                       aaTACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTcgc                                                                                                              <  1:4764871/150‑1 (MQ=255)
                                           ccGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGaaa                                                                                                          <  2:500596/150‑1 (MQ=255)
                                            cgcgTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGGTAATACGCTGTGCTCTGGCTGAGCGCATTATCTGTGGGAGCGGCCTTGTGTCGCGATAg                                                                                                         <  2:4671106/150‑1 (MQ=17)
                                                       gAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGGCACAACGCCGCGCCCCCACATAATTCGCTCAGCCCGGGCACAGCGGATTAACAGTGGGGGGGGCCGTGTGTCGCGAAAGGGCTGCAACGc                                                                                              >  2:1920695/1‑150 (MQ=255)
                                                          cGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGt                                                                                           <  2:395790/150‑1 (MQ=255)
                                                            gctgctCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcc                                                                                         <  2:291970/150‑1 (MQ=255)
                                                             ctgctCCACCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCCCCCACATATAATGCGCTCAGCCAGACCACAGCGTATTACCTGGGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTccc                                                                                        <  1:4549981/150‑1 (MQ=39)
                                                             ctgctCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTccc                                                                                        <  2:4514771/150‑1 (MQ=255)
                                                                ctCCAACGGCTGACACAATGCAATCGGGGCCCCTGCGGGCCCCAAACCCCACCCACGCCCGCTCCCACAGATAATCCGCACACCCAGACCACACCGTATTACCTGTGGCAGCGCCCTTGTGTCGCGAAAGGACTCCAACGCAGTcctttc                                                                                     <  2:188767/150‑5 (MQ=255)
                                                                 tCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGTGACAAGGGCGCACCCACAGTTAATGCGGTCTGGCAGAGCACAGCGTGTTAGCTGTGGGAGGGGCCTTGTGTGGCGGAAGGGCGGCATCGCAGTCCCCggg                                                                                    >  2:2006009/1‑149 (MQ=255)
                                                                 gccACCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCGTTGCAGTcctttcg                                                                                    <  2:1794094/149‑6 (MQ=38)
                                                                  ccAACGGCTGATACCTTGGAATCGGGGGCCCGGCGGGCCCCAAACGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGGATTACAGGTGGGGGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGtt                                                                                   >  2:4425656/1‑150 (MQ=39)
                                                                     aCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCtt                                                                                <  1:5070380/150‑1 (MQ=255)
                                                                     aCGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCtt                                                                                >  2:4685451/1‑150 (MQ=255)
                                                                      cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATGGCGACACACGGCCGCTCCCCCAGATCATGCGCTCCGCCAGAGCACAGCGTATTACCTGGGGGAGCGGCCTTGTGGCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                               <  2:3520859/150‑1 (MQ=255)
                                                                      cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACCAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                               <  2:582590/150‑1 (MQ=255)
                                                                      cGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                               <  2:3099454/150‑1 (MQ=255)
                                                                      cGGCTGATACCTTGCAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCttt                                                                               <  2:1546893/150‑1 (MQ=255)
                                                                        gCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGa                                                                             <  2:4794394/150‑1 (MQ=255)
                                                                                cTTGGATTCGGGGGCCCTGGGGGCCCCTATCGCGACACACGGCGGCTCCCACCGATACTGCGCTCAGCCAGAGCACTGGGTATTACGTGTGGGCGGGGCGTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCaacaac                                                                     <  1:3624968/150‑1 (MQ=255)
                                                                                 ttGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGGATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACAGCAACGCCGTCCCAGGTTCTTTGACCAACAACt                                                                    >  2:5536453/1‑150 (MQ=255)
                                                                                        taggggCCCGCAGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATTCGCTCAGCCAGAGCACAGCGTATTAACAGTGGGAGCGGCCTTGGGTCGCGAAAGGACTGCAACGCAGGCCCAGGTTCTTTGAACAACAACTAATccgc                                                             >  1:2737076/3‑150 (MQ=37)
                                                                                        cGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATccgc                                                             <  2:4357466/150‑1 (MQ=255)
                                                                                             gCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACt                                                        >  2:3741701/1‑150 (MQ=255)
                                                                                             gCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACt                                                        >  2:3741742/1‑150 (MQ=255)
                                                                                              cccTGCGGGCCCCTATCGCGACACACGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGTGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                       <  1:5505513/150‑1 (MQ=255)
                                                                                              cccGGCGGGCCCATACCGCCACCCCCCGCCGCCCCCAAACATAAGGCGCTCAGCCAGAGCACAGCTTTTTACCTGTGGGAGCGGCATTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACtt                                                       <  1:2990987/150‑1 (MQ=255)
                                                                                                 tGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCACGTTCTTTGACCAACAACTAATCCGCCGACTTGtt                                                    >  2:3109266/1‑150 (MQ=255)
                                                                                                     ggCCCCTATCGCGACCCCCGGCCGCTCCCACAGATACTGCGCTCCGCCAGCGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                <  1:1748478/150‑1 (MQ=255)
                                                                                                     ggCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACCGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGc                                                <  2:5541835/150‑1 (MQ=255)
                                                                                                       ccccTACCGCGACCCAAGGCCGCCCCCACATATATTGCGCTCAGCCAGAGAAACGGGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGa                                              <  1:2639132/150‑1 (MQ=255)
                                                                                                        cccTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAg                                             <  2:5097496/150‑1 (MQ=255)
                                                                                                          cTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGGCCAACAACTAATCCGCCGACTTGTTAAGCGAGTa                                           >  2:2958645/1‑150 (MQ=255)
                                                                                                          cTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTa                                           >  2:5232367/1‑150 (MQ=255)
                                                                                                                   cacaCGGCCGCTCCCACAGATAATGCGCACAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCaa                                  <  1:4064425/150‑1 (MQ=255)
                                                                                                                      aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                               >  1:5163326/1‑150 (MQ=255)
                                                                                                                      aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                               >  1:5163345/1‑150 (MQ=255)
                                                                                                                      aaGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACtt                               >  2:4836123/1‑150 (MQ=255)
                                                                                                                         gCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACttgtt                            <  2:2368962/150‑1 (MQ=255)
                                                                                                                                ccACAGATACTGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                     <  2:2934913/150‑1 (MQ=255)
                                                                                                                                ccACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAAtt                     >  1:1963731/1‑150 (MQ=255)
                                                                                                                                 cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                    >  2:2824536/1‑150 (MQ=255)
                                                                                                                                 cacaGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTg                    >  1:922110/1‑150 (MQ=255)
                                                                                                                                     gATATTGCGCTCAGCCAGTGAACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgt                <  1:280093/150‑1 (MQ=255)
                                                                                                                                      aTAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCgtg               <  1:5511625/150‑1 (MQ=255)
                                                                                                                                        aaTGCGCACACCCACAGCACAGCTTATTACCTGTGGCAGCGGCCTTGTCTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGtgtg             <  1:2563815/150‑1 (MQ=255)
                                                                                                                                          tGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACcaac                                                                        <  1:875102/89‑1 (MQ=255)
                                                                                                                                               tCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAaca      >  2:4640163/1‑150 (MQ=255)
                                                                                                                                                 aGCCAGAGCACAGCGTTTTACGTGGGGGCGCGGCGTTGGGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACAcc    <  1:4335300/150‑1 (MQ=255)
                                                                                                                                                   ccAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa  >  1:1279924/1‑150 (MQ=255)
                                                                                                                                                   ccAGACCACACCGTATAACCTGTGGGACCCGCCTTGTGACGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCa  <  1:218940/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGACCTACGAGGACTCCCAGGGCCAGAAGCACATTCTCGAATACCGCGTGATGGGGAACGGCTGCTCCAACGGCTGATACCTTGGAATCGGGGGCCCTGCGGGCCCCTATCGCGACACAAGGCCGCTCCCACAGATAATGCGCTCAGCCAGAGCACAGCGTATTACCTGTGGGAGCGGCCTTGTGTCGCGAAAGGACTGCAACGCAGTCCCAGGTTCTTTGACCAACAACTAATCCGCCGACTTGTTAAGCGAGTAATCCCTCAACTTGTTGGCAATTGTCGTGTGCGAAACACCCA  >  NC_002947/5104835‑5105131

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