Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,946 Δ1 bp 22.2% intergenic (+23/‑39) leuS → / → PP_4795 leucine‑‑tRNA ligase/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,458,9460T.22.2% 198.9 / 61.7 81intergenic (+23/‑39)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (32/31);  new base . (8/10);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

CCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGACGT  >  NC_002947/5458799‑5459094
                                                                                                                                                      |                                                                                                                                                    
ccAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGAtg                                                                                                                                                     <  2:960336/150‑2 (MQ=255)
 cAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGc                                                                                                                                                   <  1:215763/150‑1 (MQ=255)
   gCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCgg                                                                                                                                                 <  2:642673/150‑1 (MQ=255)
    ccGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCggg                                                                                                                                                >  1:947040/1‑150 (MQ=255)
      gcgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGtt                                                                                                                                              <  1:1342564/150‑2 (MQ=255)
       cgAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                             <  1:1740047/150‑1 (MQ=255)
              gTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTGGGCgtttt                                                                                                                                      <  1:538976/150‑5 (MQ=255)
              gTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAAc                                                                                                                                      >  1:1529022/1‑150 (MQ=255)
                 gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                                   >  2:864666/1‑150 (MQ=255)
                 gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                                   >  2:1459482/1‑150 (MQ=255)
                       gCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCaaa                                                                                                                             >  1:1898467/1‑150 (MQ=255)
                            ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                                   >  1:1512380/1‑139 (MQ=255)
                            ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                                   <  2:1512380/139‑1 (MQ=255)
                            ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                        <  2:1898467/150‑1 (MQ=255)
                            ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                        >  1:743383/1‑150 (MQ=255)
                             cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                                       <  1:1459482/150‑1 (MQ=255)
                             cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                                       <  2:206655/150‑1 (MQ=255)
                             cGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTAGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAg                                                                                                                       <  1:1581571/150‑1 (MQ=255)
                              gcagcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGc                                                                                                                      >  1:578391/1‑150 (MQ=255)
                                agcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAATCGCACACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCTCAAAATGGAGcaa                                                                                                                    <  2:1529022/150‑1 (MQ=255)
                                 gcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                             >  2:1154058/1‑124 (MQ=255)
                                 gcaACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                             <  1:1154058/124‑1 (MQ=255)
                                   aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGaa                                                                                                                                       <  2:745029/128‑1 (MQ=255)
                                   aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGaa                                                                                                                                       >  1:745029/1‑128 (MQ=255)
                                      gagaATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACAt                                                                                                              <  2:894288/150‑1 (MQ=255)
                                         aaTGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTgg                                                                                                                                             >  1:1745279/1‑116 (MQ=255)
                                         aaTGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTgg                                                                                                                                             <  2:1745279/116‑1 (MQ=255)
                                                tGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGc                                                                                                    >  2:579752/1‑150 (MQ=255)
                                                       cATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGc                                                                                             <  1:864666/150‑1 (MQ=255)
                                                                ccTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                             >  1:1636176/1‑93 (MQ=39)
                                                                ccTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                             <  2:1636176/93‑1 (MQ=39)
                                                                  tGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGc                                                                                  <  2:578391/150‑1 (MQ=255)
                                                                      cacccGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              >  1:1939896/4‑150 (MQ=255)
                                                                      cATCCGCGAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                                       <  2:1473857/125‑1 (MQ=255)
                                                                      cATCCGCGAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCaaa                                                                                                       >  1:1473857/1‑125 (MQ=255)
                                                                      cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTGGGCgtttt                                                                                                                                      >  2:1022174/1‑90 (MQ=37)
                                                                      cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTGGGCgtttt                                                                                                                                      <  1:1022174/94‑5 (MQ=37)
                                                                      cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              >  1:1353409/1‑150 (MQ=255)
                                                                      cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              <  2:947040/150‑1 (MQ=255)
                                                                             aaGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCAt                                                                       >  2:1057654/1‑150 (MQ=255)
                                                                               ggTCATCGTGTTACCGTGGAAGCTTGTCAACATCGTCGCCAACTGATGAAAACGCCCACCCACCCACCGA‑‑GGGCGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGc                                                                     <  1:611802/150‑1 (MQ=255)
                                                                                  cATCGTGGTACCGGGCAAGCTGGTCAATATCGTCGCCAACTTATGAAAACACCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCTCAAAATGGAGCAACAACATGATCAAACGCAATTCGCTGGTAATCGGCCTGGCGGTCATGCTCa                                                                  <  2:376475/150‑1 (MQ=255)
                                                                                    tCGTGGTACCTGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGCAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                                <  1:252378/150‑1 (MQ=255)
                                                                                    tCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAgc                                                                <  2:1939896/150‑1 (MQ=255)
                                                                                      gTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                                         <  1:1029715/91‑1 (MQ=255)
                                                                                      gTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAggg                                                                                                                         >  2:1029715/1‑91 (MQ=255)
                                                                                          tACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGc                                                          <  1:1057654/150‑1 (MQ=255)
                                                                                          tACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGc                                                          >  1:608742/1‑150 (MQ=255)
                                                                                            ccGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCAGCC‑‑CAAACGAATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACCACATGATCGCACACAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCgg                                                        >  1:1506029/1‑150 (MQ=255)
                                                                                               ggCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGttt                                                     >  2:202391/1‑150 (MQ=255)
                                                                                                gCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTc                                                    >  1:697422/1‑150 (MQ=255)
                                                                                                   aGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAg                                                 <  1:202391/150‑1 (MQ=255)
                                                                                                     cTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑TGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                                    <  1:1457974/129‑1 (MQ=255)
                                                                                                     cTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑TGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                                    >  2:1457974/1‑129 (MQ=255)
                                                                                                      tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTg                                              >  1:1279821/1‑150 (MQ=255)
                                                                                                      tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTg                                              >  1:376710/1‑150 (MQ=255)
                                                                                                        gTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcg                                            >  1:1457144/1‑150 (MQ=255)
                                                                                                         tCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcgc                                           >  2:445381/1‑150 (MQ=255)
                                                                                                         tCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATGTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTgcgc                                           >  2:218769/1‑150 (MQ=255)
                                                                                                          cAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGcgcg                                          <  2:376710/150‑1 (MQ=255)
                                                                                                             cATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCa                                       >  2:1534541/1‑150 (MQ=255)
                                                                                                                  tcgCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGgc                                  <  2:1457144/150‑1 (MQ=255)
                                                                                                                            atgaAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              >  1:652323/1‑96 (MQ=255)
                                                                                                                            atgaAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              <  2:652323/96‑1 (MQ=255)
                                                                                                                            atgaAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCGTGCGGTGTCCAGCTGCGCGGCACCGGCGCCCCTGAGc                        >  2:350083/1‑150 (MQ=255)
                                                                                                                               aaaaCGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCCCGATGTAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGa                     <  1:1180912/150‑1 (MQ=255)
                                                                                                                                ttccGCCCACC‑‑CGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGtt                                                      <  2:1340233/113‑1 (MQ=39)
                                                                                                                                ttccGCCCACC‑‑CGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGtt                                                      >  1:1340233/4‑116 (MQ=39)
                                                                                                                                  aCGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCg                  <  2:608742/150‑1 (MQ=255)
                                                                                                                                   cGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTTATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGt                 <  2:1252620/150‑1 (MQ=255)
                                                                                                                                   cGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGt                 <  1:218769/150‑1 (MQ=255)
                                                                                                                                     cccACC‑‑CGCAATCGATGCTGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATg                                                                      <  1:951642/95‑1 (MQ=255)
                                                                                                                                     cccACC‑‑CGCAATCGATGCTGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATg                                                                      >  2:951642/1‑95 (MQ=255)
                                                                                                                                     cccACC‑‑CGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGa               >  2:862284/1‑150 (MQ=255)
                                                                                                                                      ccACC‑‑CGCAATCGATGCGGGTGGGCGTAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCCGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCg                  >  2:458591/1‑146 (MQ=255)
                                                                                                                                      ccACC‑‑CGCAATCGATGCGGGTGGGCGTAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCCGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCg                  <  1:458591/146‑1 (MQ=255)
                                                                                                                                       cACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAg             >  2:147835/1‑150 (MQ=255)
                                                                                                                                         cc‑‑cGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              <  1:829253/83‑1 (MQ=39)
                                                                                                                                         cc‑‑cGCAATCGATGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                              >  2:829253/1‑83 (MQ=39)
                                                                                                                                              gCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAt        >  1:753116/1‑150 (MQ=255)
                                                                                                                                                   cGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGAc    <  2:1353409/150‑1 (MQ=255)
                                                                                                                                                     aTTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGACGt  >  1:359058/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
CCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAATGGACGT  >  NC_002947/5458799‑5459094

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: