Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,458,940:1 +A 22.2% intergenic (+17/‑45) 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,9401.A22.2% 198.7 / 78.1 81intergenic (+17/‑45)leuS/PP_4795leucine‑‑tRNA ligase/lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (31/32);  new base A (8/10);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 7.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

GCCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458792‑5459088
                                                                                                                                                       |                                                                                                                                                    
gccGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCa                                                                                                                                                     >  1:206655/1‑150 (MQ=255)
    cTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGa                                                                                                                                                 <  2:708835/150‑1 (MQ=255)
       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:1459482/1‑150 (MQ=255)
                        gAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAgg                                                                                                                             >  2:864666/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                                                                                                                  >  1:743383/1‑150 (MQ=255)
                                   ccGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGa                                                                                                                  <  2:1898467/150‑1 (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                                                                                                                                 >  1:745029/1‑128 (MQ=255)
                                          aaCGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGaa                                                                                                                                 <  2:745029/128‑1 (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                                                                                                                                       <  2:1636176/93‑1 (MQ=39)
                                                                       ccTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTgg                                                                                                                                       >  1:1636176/1‑93 (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                                                                                                                                <  1:1022174/94‑5 (MQ=37)
                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCAATCGATGCGGGTGGGCgtttt                                                                                                                                >  2:1022174/1‑90 (MQ=37)
                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                        <  2:947040/150‑1 (MQ=255)
                                                                             cATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                        >  1:1353409/1‑150 (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                                                              >  2:1457974/1‑129 (MQ=255)
                                                                                                            cTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑TGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCt                                                              <  1:1457974/129‑1 (MQ=255)
                                                                                                             tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTg                                        >  1:376710/1‑150 (MQ=255)
                                                                                                             tGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTg                                        >  1:1279821/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                                                                        <  2:652323/96‑1 (MQ=255)
                                                                                                                                   atgaAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTgg                                                                        >  1:652323/1‑96 (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)
                                                                                                                                                       |                                                                                                                                                    
GCCGCTGCCAGCCGCGAAGAAGTCGAAGCCGCTGCCCGCAGCAACGAGAATGTCCTGCGCTTCATCGATGGCCTGACCATCCGCAAGGTCATCGTGGTACCGGGCAAGCTGGTCAACATCGTCGCCAACTGATGAAAACGCCCACC‑‑CGCATCGATTGCGGGTGGGCGGAACAGGCCCAAAAGGGAGCAACAACATGATCAAACGCAATTTGCTGGTAATCGGCCTGGCGGTCATGCTCAGCGCCTGCGGTTTCCAGCTGCGCGGCACCGGCGCCACCGAGCTGAGCGTGAAGGAAAT  >  NC_002947/5458792‑5459088

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: