Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,058 A→T 13.8% intergenic (+38/+65) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0580AT13.8% 275.4 / 21.5 122intergenic (+38/+65)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base T (2/14);  total (55/67)
Fisher's exact test for biased strand distribution p-value = 2.05e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

AGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541911‑5542206
                                                                                                                                                   |                                                                                                                                                    
aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  2:909726/150‑1 (MQ=255)
aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:1839322/150‑1 (MQ=255)
aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:4582879/150‑1 (MQ=255)
aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  1:5157414/150‑1 (MQ=255)
aGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAgg                                                                                                                                                    <  2:390334/150‑1 (MQ=255)
  aTCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGt                                                                                                                                                  <  1:3986991/150‑1 (MQ=255)
     gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                               <  1:1676365/150‑1 (MQ=255)
     gATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGaa                                                                                                                                               <  1:2819310/150‑1 (MQ=255)
      aTGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAc                                                                                                                                              <  1:965358/150‑1 (MQ=255)
       tGCCTTCCTGCATGGGCCGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAcc                                                                                                                                             <  1:222127/150‑1 (MQ=255)
         ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCGGCGAATACGAGGGTGAACCCa                                                                                                                                           >  2:3037593/1‑150 (MQ=255)
         ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCGGCGAATACGAGGGTGAACCCa                                                                                                                                           >  2:3038077/1‑150 (MQ=255)
         ccTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCa                                                                                                                                           <  1:2077067/150‑1 (MQ=255)
          cTTCCTGCATGGGCAGCCGGTTGTCCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAc                                                                                                                                          <  2:637945/150‑1 (MQ=255)
           ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCGCGAATACGATGGTGAACCCagc                                                                                                                                         >  2:4161076/1‑148 (MQ=255)
           ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                         >  2:3300004/1‑150 (MQ=255)
           ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGAGCCCCGCAAAAGAGGGGGTACACAcc                                                                                                                                         >  1:1116608/1‑150 (MQ=255)
           ttCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACACAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCAcc                                                                                                                                         <  2:3285406/150‑1 (MQ=255)
              cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                      >  2:5487005/1‑150 (MQ=255)
              cTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGAc                                                                                                                                      >  2:5487073/1‑150 (MQ=255)
               tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCACGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg                                                                                                                                     <  1:4935399/150‑1 (MQ=255)
               tGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg                                                                                                                                     <  1:3073708/150‑1 (MQ=255)
               tGCATGGGCAGCCGGTGGCCCAGGCAGCGTTCGACACCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGGGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACg                                                                                                                                     <  1:3073445/150‑1 (MQ=255)
                gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCCCGAATACGAGGGTGTACCCACCCACGc                                                                                                                                    >  2:2116987/1‑150 (MQ=255)
                   tGGGCAGCCGGTGGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCAtt                                                                                                                                 <  1:193373/150‑1 (MQ=255)
                      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCGCGAAAAAGAGGGTGTAACCGCCGACGCAttagg                                                                                                                              >  2:1407211/1‑147 (MQ=255)
                      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                              <  1:3766627/150‑1 (MQ=255)
                       cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgcg                                                                                                                             <  1:1542484/150‑1 (MQ=255)
                       cAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAAACCACCGACGCATTcgcg                                                                                                                             >  2:4476792/1‑150 (MQ=255)
                        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  1:1039449/150‑1 (MQ=25)
                        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                            <  1:605955/150‑1 (MQ=25)
                        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  1:4611341/150‑1 (MQ=255)
                        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:3345262/150‑1 (MQ=255)
                        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCACCACAAGCGCGAACTGCCGTTCGCTCCTTGATAACTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:3345407/150‑1 (MQ=255)
                          ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCCCGAATAACAGGGTGAACCCACCGACCCAATAGCGGGt                                                                                                                          >  2:4390191/1‑150 (MQ=255)
                             gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGaa                                                                                                                       >  2:4145889/1‑150 (MQ=255)
                              ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGAACCCCCAAAAACAGGGGGAAACCACCCACCCAATAGCCGGTGAAc                                                                                                                      >  1:4931295/1‑150 (MQ=255)
                              ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                      <  1:1460192/150‑1 (MQ=25)
                              ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAAc                                                                                                                      <  1:1568255/150‑1 (MQ=25)
                              ttGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAAc                                                                                                                      <  2:1928662/150‑1 (MQ=255)
                                  ccAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTTCACCCGc                                                                                                                  <  2:4701742/150‑1 (MQ=255)
                                   cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGGATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCt                                                                                                                 >  2:1789231/1‑150 (MQ=255)
                                            ttCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCCggg                                                                                                        >  2:5044502/1‑150 (MQ=255)
                                               gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTc                                                                                                     >  2:4215708/1‑150 (MQ=255)
                                                  atcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                  >  2:2716692/1‑150 (MQ=255)
                                                  aacatcGTCGCCACCTCCCGCACCACCCACCACACGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCTGGAGCGGGTTCACCCGCGAATACGTGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGcc                                                                                                  <  2:2686957/148‑1 (MQ=37)
                                                            ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGAACCCCCGATAACAGGGTGTACCCCCCGCCGCCATCGCCGGTGTAACCCCGCCTGCAGGGTTGGCCGGCCGCCgtt                                                                                        >  2:1477022/1‑150 (MQ=255)
                                                              accgcccGCACGCCCCAGCCCACGCGCTACCGGCGGTTCGCTCCGTGATAAGGGGCCTGCAGGAGGGGGTTCACCGGCGAATACGAGGGTGACCCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGa                                                                                      <  2:3018209/145‑1 (MQ=25)
                                                              aCCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCAACCGCGAATACGAGGGTGAACCCACCGACCCAATCGCGGGTGAACCCGCTCCTGCAGGGTTAGCCGTCAGCCTTTGa                                                                                      >  2:5055149/1‑150 (MQ=255)
                                                                  cccGAAAGCCCCAGCACACGCGCTAACTGCCGTTCGCCCCTTGAAAAGTGGCCTGCAGGCGCGGGTTAACCCGCGATTACGAGGGTGAACCAACCGACGCTTTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTAc                                                                                  <  2:570887/149‑1 (MQ=255)
                                                                    cgcaccaCCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGTTATGTGGCCTGCAGGCGGGGGTTCACCCGCTATTCCGCGGGTTAACCCACCGACGCATTGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACtt                                                                                <  2:1052564/144‑1 (MQ=255)
                                                                     gCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGAACCCACCGACGCAATCGCGGGTGAAACCGCACCTGCAGGGTTCGCCGGCAGCCTTTGATTACTTg                                                                               >  2:4140136/1‑150 (MQ=255)
                                                                       aacaCCCCGAACACGCGCGAACTGCCGTTCGCTCCTTCATACGGGGCCTGCAGGCGCGGGTACACCCGCTAATACGCGGGTGACCCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                             <  2:2835113/147‑1 (MQ=37)
                                                                          ccccAGCACACGCGCGAACTGCGGTTCGCTCCTTTAAACGTGGCCTGCAGGTGCGGTTTCACCCGCTAATACGAGGGTGAACCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                          <  2:2106707/149‑1 (MQ=255)
                                                                          aCCCCGCCCACGCGCGAACTGCCGTTCGCTCATTTATACTTGGCCTCCAGGCGGGGTTTCACCCGCTAATACGCGGGTGAACCCACCTACGCATGCGCGGGTTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                          <  2:2106762/150‑1 (MQ=255)
                                                                          aCCCAGCACACGCGCTAACTGCCGTTCGCTCATTGATATGTGGCCTGCTGGAGGGGTTTCACCCGCGAATACGCGGGTTAACCCACCTACGCTTTGGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAt                                                                          <  2:2106754/150‑1 (MQ=255)
                                                                            ccAGCACCAGCGCTAAGTGCCGTTCGCTCCTTGATAGGTGGCCTGCTGGAGCGTGTTCACCCGCTAATACGTGGGTCACCCCACCGCCGCTTTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGa                                                                        <  2:637311/150‑1 (MQ=255)
                                                                              aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            <  2:4113068/96‑1 (MQ=255)
                                                                              aGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                            >  1:4113068/1‑96 (MQ=255)
                                                                                       gccaagtGCGGTTCGCTCGTTGATATGTGGCCTGCTGGCGGGGGTTCCCCCGCTAATCGGTGGTTCACCCCACCGACGCTTTGGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                             <  1:1287917/144‑1 (MQ=255)
                                                                                         gAACTGCCGTTCGCTCCTTTATTAGTGGCCTGCAGGAGCGGTTTCACCCGCGAATACGAGGGTGAACCCACCGACGCTTTCGCGGTTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcat                                                           <  2:837233/150‑1 (MQ=255)
                                                                                          aaCTGCCGTTCGCTCCTTGTTATGTGGCCTGCTGGAGCGGTTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtcatc                                                          <  1:5113883/150‑1 (MQ=38)
                                                                                             tGCGGTTCGCTCCTTGATATGTGGCCTGCTGGCGCGGTTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                       <  1:4136881/150‑1 (MQ=38)
                                                                                             tGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAAAACGAGGGTGAACCCACCCACGCATTCGCGGGTGAACCCCCTCCAGCACGGTTCGCCCTCAGCCTTTGATTACTTGACGATGACCGTTCCGTTAatcatc                                                       >  2:3738588/1‑150 (MQ=38)
                                                                                             tGCCGTGCGCTCTTTTATACGTGGCCTGCCGGCGCGGGTCCCCCCGAGAATCCGCGGGACACCCCACCGCCGCTTTGGGGGGCTAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                       <  1:3370679/150‑1 (MQ=255)
                                                                                              gCCGTTCGCTCCTTGATAAGTGGCCTGCTGGAGGGGGTTCACCCGCTAATACGTGGTTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                      <  1:4679971/150‑1 (MQ=39)
                                                                                              gCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                      <  2:3743065/150‑1 (MQ=255)
                                                                                                              tctgTGGCCGGCTGGCGGGGTTCCCCCCGCTATTCCGGGGTTCAACCCCCCGCCGCATTCGGGGGTCACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGaa                                      <  1:2083588/147‑1 (MQ=255)
                                                                                                                 gTGGCCTGCTGGCGGGTTGTCACCCCCGATTACGTGGGTAAACCCCCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGt                                   <  1:1407211/150‑1 (MQ=255)
                                                                                                                   ggCCTGCAGGAGCGGTTCCACCCGCTAATCCGCGGGTCAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGc                                 <  2:1121030/150‑1 (MQ=255)
                                                                                                                         ccggcgCGGGTTCACCCGCGAATACGCGGGTGAACCCCCCGACGCATTCGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                           <  2:4223301/145‑1 (MQ=255)
                                                                                                                         cAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGTCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCagaaga                           <  2:4620492/150‑1 (MQ=255)
                                                                                                                           ggtgCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                         <  2:330540/147‑1 (MQ=255)
                                                                                                                            gcgCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGc                        <  2:1631683/148‑1 (MQ=255)
                                                                                                                                ggtttCCCCGGCTAATCCGTGGGTTAACCCACCGCCGCTTGCGGGGTTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                    <  1:2396575/147‑1 (MQ=255)
                                                                                                                                gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                    >  1:2044872/1‑150 (MQ=255)
                                                                                                                                  gTTCACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                  >  2:1398056/1‑150 (MQ=255)
                                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                  >  1:5431820/1‑150 (MQ=255)
                                                                                                                                   ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                 >  1:997555/1‑150 (MQ=255)
                                                                                                                                   ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTc                 >  1:997557/1‑150 (MQ=255)
                                                                                                                                     cACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTcac               >  1:4760007/1‑150 (MQ=25)
                                                                                                                                       cccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGCTGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca             >  2:2495275/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCGGGGACGGTGCAGAAGTAGCTGTAGTcaccac            >  2:3515020/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:1023051/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:349447/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:2244828/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:5073629/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:495465/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:2173043/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:1596616/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:625855/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:1467427/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:1290767/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGAGTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  1:2829943/1‑150 (MQ=255)
                                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCAGTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac            >  2:2406062/1‑150 (MQ=255)
                                                                                                                                         cgcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCaccacc           >  2:5325079/1‑150 (MQ=255)
                                                                                                                                               tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc     <  1:1789231/150‑1 (MQ=255)
                                                                                                                                                aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc    >  1:4399320/1‑150 (MQ=255)
                                                                                                                                                aCGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTcc    >  2:5493535/1‑150 (MQ=255)
                                                                                                                                                  gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg  >  1:830839/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGATCGATGCCTTCCTGCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAG  >  NC_002947/5541911‑5542206

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