Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,868,536 A→C 8.3% intergenic (+32/‑899) 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,5360AC 8.3% 185.2 / 10.0 73intergenic (+32/‑899)PP_2520/PP_2522hypothetical protein/transposase
Reads supporting (aligned to +/- strand):  ref base A (36/31);  new base C (3/3);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

AATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCGGCGAACC  >  NC_002947/2868387‑2868685
                                                                                                                                                     |                                                                                                                                                     
tctGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGACTGTACAAGCCCTCCTGCGCTGCTGGCGTTGTCAGCCTTCCTGGGGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCaa                                                                                                                                                       <  2:1653854/148‑1 (MQ=255)
    tCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCaagaga                                                                                                                                                   >  1:1653854/1‑148 (MQ=255)
        gTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTg                                                                                                                                               >  2:1532258/1‑150 (MQ=255)
        gTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTg                                                                                                                                               >  1:122210/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:929537/150‑1 (MQ=255)
                          gcagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTgg                                                                                                                             >  1:1833059/1‑150 (MQ=255)
                          gcagcaAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTgg                                                                                                                             <  2:820975/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                                                                                                                              <  1:367698/91‑1 (MQ=255)
                                                                                    ggTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTg                                                                                                                              >  2:367698/1‑91 (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                                                             <  2:777858/119‑1 (MQ=255)
                                                                                                                         tgtgATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCtt                                                             >  1:777858/1‑119 (MQ=255)
                                                                                                                           tgATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTcc                            >  2:1253955/1‑150 (MQ=255)
                                                                                                                                    ccGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCgctggc                   >  1:1407106/1‑150 (MQ=255)
                                                                                                                                    ccGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCgctggc                   >  2:424741/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                                                                    <  2:275439/91‑1 (MQ=255)
                                                                                                                                              ttCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTAc                                                                    >  1:275439/1‑91 (MQ=255)
                                                                                                                                              ttCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCg         >  1:1385552/1‑150 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGt                                                                  >  2:1779100/1‑87 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGt                                                                  <  1:1779100/87‑1 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGAc                                    <  1:952388/117‑1 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGAc                                    >  2:952388/1‑117 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGATGACTACATTCCCGCGCTGGCTGACGTGc           >  1:1555471/1‑142 (MQ=255)
                                                                                                                                                    aagaagACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGATGACTACATTCCCGCGCTGGCTGACGTGc           <  2:1555471/142‑1 (MQ=255)
                                                                                                                                                     agaagaCTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCGGCGAAcc  >  2:1279620/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATGTCGCGTCAAGCTGTGCCGCTCGGCAGCAAGCTGTTGGCCTGTAAAAGACCTCATGCGCTGCTGGCGTTGTCAGGCTTCCTGGTGATGATGACGCTTATGGTGATTTGTTTCTGAATGTGTGATCTCTCCCGCTCAGTCTTCGTCAAGAAGACTGCCCCTCGCCGTCGTCTGGGGATTGATACAGGGCCTTCCATGCAAACGCTGTTGAACGAGATTCTCGAAGAAGTACGTCCCTTGATCGGCAAGGGCAAAGTCGCTGACTACATTCCCGCGCTGGCTGACGTGCCGGCGAACC  >  NC_002947/2868387‑2868685

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: