Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,532 G→T 7.9% intergenic (+28/‑903) PP_2520 → / → PP_2522 hypothetical protein/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,868,5320GT 7.9% 235.9 / 6.8 77intergenic (+28/‑903)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base G (39/31);  new base T (3/3);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

GCTAATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCG  >  NC_002947/2868384‑2868678
                                                                                                                                                    |                                                                                                                                                  
gcTAATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGt                                                                                                                                                   >  2:379075/1‑150 (MQ=255)
 cTAATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTc                                                                                                                                                  >  2:374259/1‑150 (MQ=255)
   tctGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGACTGTACAAGCCCTCCTGCGCTGCTGGCGTTGTCAGCCTTCCTGGGGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCaa                                                                                                                                                <  2:1653854/148‑1 (MQ=255)
       tCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCaagaga                                                                                                                                            >  1:1653854/1‑148 (MQ=255)
           gTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTg                                                                                                                                        >  1:122210/1‑150 (MQ=255)
           gTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTg                                                                                                                                        >  2:1532258/1‑150 (MQ=255)
            tCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGc                                                                                                                                       <  1:374259/150‑1 (MQ=255)
             cAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCTTGACGAAGACTgag                                                                                                                                      <  2:1520235/150‑3 (MQ=255)
                gCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCt                                                                                                                                   >  1:820975/1‑150 (MQ=255)
                 cTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTAAGTCTTCGTCAAGAAGACTGCCCCTc                                                                                                                                  <  2:1641249/150‑1 (MQ=255)
                          tCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTAAGGTGATTTTTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCgtcgtc                                                                                                                         >  1:703167/1‑150 (MQ=255)
                           cGGCAGCAAGCTGTTGGCCTGTAAAAGACCTAATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTATTTCTGAATGTGTAATCTCACCCGCTCAGTCTTCCTCAATGAGACTTCCCCTCGCCGTCGTCt                                                                                                                        <  2:498445/150‑1 (MQ=255)
                             gcagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTgg                                                                                                                      <  2:820975/150‑1 (MQ=255)
                             gcagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTgg                                                                                                                      >  1:1833059/1‑150 (MQ=255)
                             gcagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTgg                                                                                                                      <  2:929537/150‑1 (MQ=255)
                              cagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTggg                                                                                                                     <  1:379075/150‑1 (MQ=255)
                              cagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTggg                                                                                                                     <  2:1718995/150‑1 (MQ=255)
                                 caAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGAt                                                                                                                  >  1:1636405/1‑150 (MQ=255)
                                  aaGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCCGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGAtt                                                                                                                 <  2:1633938/150‑1 (MQ=255)
                                      tGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATa                                                                                                             >  1:11445/1‑150 (MQ=255)
                                       gTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATAc                                                                                                            >  1:1501174/1‑150 (MQ=255)
                                             cTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGcc                                                                                                      >  2:308573/1‑150 (MQ=255)
                                             cTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCATATGGTGATTTGTTTCTGAATGTGTGATCTCGCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGcc                                                                                                      <  2:1438869/150‑1 (MQ=255)
                                              tGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCt                                                                                                     >  1:1461459/1‑150 (MQ=255)
                                              tGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCt                                                                                                     >  1:867439/1‑150 (MQ=255)
                                                    aGACCTCATGCGCTACTGGCGTTGTCAGGCTTCCTGGTGATGATGACGTTTATGGTTATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGTATTGATACAGGGCCTTCCATg                                                                                               >  1:1279810/1‑150 (MQ=255)
                                                      aCCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATg                                                                                               >  1:1001335/1‑148 (MQ=255)
                                                      aCCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCCTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATg                                                                                               <  2:1001335/148‑1 (MQ=255)
                                                           aTGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGc                                                                                        <  2:1501174/150‑1 (MQ=255)
                                                           aTGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGc                                                                                        >  1:1547890/1‑150 (MQ=255)
                                                              cgcTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACATAAATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGt                                                                                     <  1:1841650/150‑1 (MQ=255)
                                                                  gctgGCGTTGTCAGGCTTCCTGGTGCTGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCATTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGCCCTTCCATGCAAACGCTGTTGaa                                                                                 <  2:948665/150‑1 (MQ=255)
                                                                     ggCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACga                                                                              <  1:308573/150‑1 (MQ=255)
                                                                      gCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACgag                                                                             <  1:1619676/150‑1 (MQ=255)
                                                                      gCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCAATGCAAACGCTGTTGAACgag                                                                             <  2:122210/150‑1 (MQ=255)
                                                                        gTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAAc                                                                                >  2:178500/1‑145 (MQ=255)
                                                                        gTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAAc                                                                                <  1:178500/145‑1 (MQ=255)
                                                                              aGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCga                                                                     >  2:1301918/1‑150 (MQ=255)
                                                                                gCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCgaag                                                                   <  1:1532258/150‑1 (MQ=255)
                                                                                      tGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCATCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGCGATTCTCGAAGAAGTAc                                                             >  2:584204/1‑150 (MQ=255)
                                                                                       ggTGCTGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCTTGACGAAGACTgag                                                                                                                                      <  1:133140/76‑3 (MQ=37)
                                                                                       ggTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCTTGACGAAGACTgag                                                                                                                                      >  2:133140/1‑74 (MQ=38)
                                                                                       ggTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTg                                                                                                                       >  2:367698/1‑91 (MQ=255)
                                                                                       ggTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTg                                                                                                                       <  1:367698/91‑1 (MQ=255)
                                                                                          gatgatgaCGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTcc                                                         >  2:1216069/1‑150 (MQ=255)
                                                                                              atgaCGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCt                                                                                       >  2:1854196/1‑116 (MQ=255)
                                                                                              atgaCGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCt                                                                                       <  1:1854196/116‑1 (MQ=255)
                                                                                                     ttATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGAt                                                   >  2:957284/1‑145 (MQ=255)
                                                                                                     ttATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGAt                                                   <  1:957284/145‑1 (MQ=255)
                                                                                                        tGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAgg                                           <  1:1544519/150‑1 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCtt                                                      <  2:600087/137‑1 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCtt                                                      >  1:600087/1‑137 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGc                                         >  1:132663/1‑150 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGc                                         >  2:233455/1‑150 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCACCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGc                                         >  2:261015/1‑150 (MQ=255)
                                                                                                          gTGATTTGTTTCTGAATGTGTGATCTCTCCCGATCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGc                                         >  1:1095987/1‑150 (MQ=255)
                                                                                                           tGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAgg                                                                                                         <  2:331732/85‑1 (MQ=255)
                                                                                                           tGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAgg                                                                                                         >  1:331732/1‑85 (MQ=255)
                                                                                                             aTTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCaaa                                      <  1:1539314/150‑1 (MQ=255)
                                                                                                                   ttCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCt                                <  2:1267190/150‑1 (MQ=255)
                                                                                                                   ttCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCt                                <  2:132663/150‑1 (MQ=255)
                                                                                                                    tCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTg                               >  2:532362/1‑150 (MQ=255)
                                                                                                                     cTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCATCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGa                              <  1:584204/150‑1 (MQ=255)
                                                                                                                        aaTGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGAATCTCGAAGAAGGACGTCCCTTGATCGGCAAGGGCACAGTAGCTGACTa                           >  2:312244/1‑150 (MQ=255)
                                                                                                                            tgtgATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCtt                                                      >  1:777858/1‑119 (MQ=255)
                                                                                                                            tgtgATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCtt                                                      <  2:777858/119‑1 (MQ=255)
                                                                                                                              tgATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTcc                     >  2:1253955/1‑150 (MQ=255)
                                                                                                                                       ccGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCgctggc            >  2:424741/1‑150 (MQ=255)
                                                                                                                                       ccGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCgctggc            >  1:1407106/1‑150 (MQ=255)
                                                                                                                                           tCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGCCGAGGTACTTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGAc        <  1:233455/150‑1 (MQ=255)
                                                                                                                                           gcAGTCTTCTTGACGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGt                                                                                     >  1:917462/2‑73 (MQ=38)
                                                                                                                                           gcAGTCTTCTTGACGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGt                                                                                     <  2:917462/72‑1 (MQ=38)
                                                                                                                                           gcAGTCTTCTTGACGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGAc        >  2:956102/2‑150 (MQ=255)
                                                                                                                                              gTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTg     >  1:264150/1‑150 (MQ=255)
                                                                                                                                                 ttCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTAc                                                             >  1:275439/1‑91 (MQ=255)
                                                                                                                                                 ttCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTAc                                                             <  2:275439/91‑1 (MQ=255)
                                                                                                                                                 ttCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCg  >  1:1385552/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCTAATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCG  >  NC_002947/2868384‑2868678

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