Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,706,643 T→A 7.1% L202Q (CTG→CAG)  PP_5009 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,706,6430TA 7.1% 287.9 / 11.9 99L202Q (CTG→CAG) PP_5009hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (49/43);  new base A (4/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

GGGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGAC  >  NC_002947/5706495‑5706791
                                                                                                                                                    |                                                                                                                                                    
gggCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTg                                                                                                                                                     <  1:2401301/150‑1 (MQ=255)
    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCACTGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTg                                                                                                                                                 >  2:42404/1‑150 (MQ=255)
    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTg                                                                                                                                                 >  1:1257435/1‑150 (MQ=255)
        cgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCCAGGTATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGa                                                                                                                                             <  1:392055/150‑1 (MQ=255)
         gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgtg                                                                                                                                            <  2:110639/150‑3 (MQ=255)
           aGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                          <  1:2065028/150‑1 (MQ=255)
           aGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                          <  1:2439734/150‑1 (MQ=255)
               aGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                                    <  2:678090/136‑1 (MQ=255)
               aGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                                    >  1:678090/1‑136 (MQ=255)
               aGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCCAGGCGATa                                                                                                                                      >  2:372208/1‑150 (MQ=255)
                   aaGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAg                                                                                                                                  >  1:2046852/1‑150 (MQ=255)
                     gTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTc                                                                                                                                <  1:1500853/150‑1 (MQ=255)
                        tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACg                                                                                                                             >  1:753257/1‑150 (MQ=255)
                        tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAAATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACg                                                                                                                             >  2:864744/1‑150 (MQ=255)
                            gcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCg                                                                                                                         >  1:2185761/1‑150 (MQ=255)
                              aTGCTCTGGCTGGAGCACACCATCAAGTCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGa                                                                                                                       >  1:568458/1‑150 (MQ=255)
                                 gtctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACg                                                                                                                    >  1:1757846/2‑150 (MQ=255)
                                 ctctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACg                                                                                                                    >  1:821242/1‑150 (MQ=255)
                                 ctctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACg                                                                                                                    >  2:1302358/1‑150 (MQ=255)
                                   cgggccggAGCACACCATCAAGCCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACANGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                                  <  1:2037522/148‑1 (MQ=255)
                                    tggctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGAt                                                                                                                 <  1:2074166/150‑1 (MQ=255)
                                    tggctcgATCACAACATCAAGGCTGCATTTCTCGATAACCTCCCTCGCCAAGGCCTCGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGAt                                                                                                                 <  2:1124720/150‑1 (MQ=255)
                                        tggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                             >  2:987013/1‑150 (MQ=255)
                                         ggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgg                                                                                                            <  2:192655/150‑1 (MQ=255)
                                         ggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgg                                                                                                            <  2:821242/150‑1 (MQ=255)
                                          gAGCACACCATCAAGGCTGCATTGGTCGAGACCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                           <  2:1081974/150‑1 (MQ=255)
                                          gAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCCAGACCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                           <  1:258515/150‑1 (MQ=255)
                                             cacaCCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGca                                                                                                        <  1:239528/150‑1 (MQ=255)
                                              acacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcac                                                                                                       >  2:2450081/1‑150 (MQ=255)
                                              acacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGATATTAAGGACGGGGTACACGACAGCCTGGCTGTCGAGGCGAAACAAGGCTCGCTCGGAACGGATGGCGGGGcac                                                                                                       >  1:413599/1‑150 (MQ=255)
                                               cacCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcaca                                                                                                      <  1:2450081/150‑1 (MQ=255)
                                                 ccANCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcacaca                                                                                                    >  2:1498829/1‑150 (MQ=255)
                                                      aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                                    >  1:2326906/1‑97 (MQ=255)
                                                      aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                                    <  2:2326906/97‑1 (MQ=255)
                                                      aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGc                                                                                                                                                   <  1:1389243/98‑1 (MQ=255)
                                                      aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGc                                                                                                                                                   >  2:1389243/1‑98 (MQ=255)
                                                       aGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTa                                                                                              <  2:1757846/150‑1 (MQ=255)
                                                        ggCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAc                                                                                             <  1:864744/150‑1 (MQ=255)
                                                        ggCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAc                                                                                             >  2:2347822/1‑150 (MQ=255)
                                                        ggCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAc                                                                                             >  2:40489/1‑150 (MQ=255)
                                                              aTTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTAcc                                                                                       >  1:307892/1‑150 (MQ=255)
                                                                   tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                                  >  1:27766/1‑150 (MQ=255)
                                                                   tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                                  >  1:2269213/1‑150 (MQ=255)
                                                                   tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                                  >  1:1020058/1‑150 (MQ=255)
                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                                  >  2:1122530/1‑116 (MQ=255)
                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                                  <  1:1122530/116‑1 (MQ=255)
                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTa                                                                                >  1:1581924/1‑150 (MQ=255)
                                                                         aCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGc                                                                            <  2:2185761/150‑1 (MQ=255)
                                                                          cctccctcGCCAAGGCCTGGCAGGTCGGGCAAAAAGTGCAAGGGACGCCCGGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                           >  1:1233615/1‑150 (MQ=255)
                                                                           tcgcctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                        <  1:430795/85‑1 (MQ=255)
                                                                           tcgcctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                        >  2:430795/4‑88 (MQ=255)
                                                                           ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtc                                                                          <  1:987013/150‑1 (MQ=255)
                                                                           ctccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtc                                                                          >  1:2246478/1‑150 (MQ=255)
                                                                             ccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcat                                                                        <  2:1827364/150‑1 (MQ=255)
                                                                                  gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAc                                                                                             >  2:235773/1‑124 (MQ=255)
                                                                                  gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAc                                                                                             <  1:235773/124‑1 (MQ=255)
                                                                                  gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                                   >  2:1660279/1‑150 (MQ=255)
                                                                                    cAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTa                                                                                              >  1:1650998/1‑121 (MQ=255)
                                                                                    cAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTa                                                                                              <  2:1650998/121‑1 (MQ=255)
                                                                                        gCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                             <  1:1302358/150‑1 (MQ=255)
                                                                                                  tgatgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaa                                                   <  2:753257/150‑1 (MQ=255)
                                                                                                       cAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTgg                                              >  1:1282114/1‑150 (MQ=255)
                                                                                                         agaGCTGCCATGTCCGCCCGGCTTGACGAGATACACAACAGCCTGGCCGTCGAGGCTATTGCTGTTACGCTCGGAACGGATGGCGGGGCACACACATTCACCGCAACGTTTCTATTGCTCATCCCCACAGGTAGCCTCCGCGACCTGGtt                                            <  2:568458/148‑1 (MQ=255)
                                                                                                          aaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACCAACCGGTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTc                                           >  1:2232101/1‑150 (MQ=255)
                                                                                                           aGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcac                                                                     <  2:1139283/123‑1 (MQ=255)
                                                                                                           aGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcac                                                                     >  1:1139283/1‑123 (MQ=255)
                                                                                                              tGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgca                                       >  2:1572659/1‑150 (MQ=255)
                                                                                                                aagtGTTCGCCCTGCTCGACGGGATAGACGACTGCCTGGCGGCCGAGGCGATACAACCCACGCTCGAAACGGATGGCTGGGCACACCCCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCTCAACCGGGTTCtgcaat                                     <  2:814272/147‑1 (MQ=255)
                                                                                                                        gCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGccc                             <  2:1581924/150‑1 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                          >  2:601221/1‑55 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                          <  1:601221/55‑1 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                     <  2:421656/140‑1 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaat                                     >  1:421656/1‑140 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGc                               >  1:1611784/1‑146 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGc                               <  2:1611784/146‑1 (MQ=255)
                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTg                           >  2:2349603/1‑150 (MQ=255)
                                                                                                                            tGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcg                         <  1:40489/150‑1 (MQ=255)
                                                                                                                            tGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcg                         <  1:2349603/150‑1 (MQ=255)
                                                                                                                              cTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCCATGCAAGCCCTGcggt                       >  1:181945/1‑150 (MQ=255)
                                                                                                                               tCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGtggtt                      <  2:2269213/150‑1 (MQ=255)
                                                                                                                                 gACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                             >  2:2311882/1‑61 (MQ=255)
                                                                                                                                 gACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                             <  1:2311882/61‑1 (MQ=255)
                                                                                                                                 gACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACAACGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcg                    >  2:865482/1‑150 (MQ=255)
                                                                                                                                  aCGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcgg                   <  1:1660279/150‑1 (MQ=255)
                                                                                                                                      ggTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAGCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGgttcggttcg               >  1:1340026/1‑150 (MQ=255)
                                                                                                                                        tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc             >  1:399432/1‑150 (MQ=255)
                                                                                                                                        tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc             <  1:2347822/150‑1 (MQ=255)
                                                                                                                                        tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc             <  2:413599/150‑1 (MQ=255)
                                                                                                                                        tACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcctc             >  2:329736/1‑150 (MQ=255)
                                                                                                                                        gcctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc             >  1:89539/5‑150 (MQ=255)
                                                                                                                                         cctcGACAGCCAGGCTGTCGTGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg            >  2:966529/4‑150 (MQ=255)
                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                             >  2:206534/3‑100 (MQ=255)
                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                             <  1:206534/98‑1 (MQ=255)
                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc           >  2:1917161/3‑150 (MQ=255)
                                                                                                                                             gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                             <  1:2054613/97‑1 (MQ=255)
                                                                                                                                             gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                             >  2:2054613/1‑97 (MQ=255)
                                                                                                                                                aGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc     <  2:89539/150‑1 (MQ=255)
                                                                                                                                                  ccTGGCTGTCGAGGCGATACAACTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCCACCGCCACCTGGTTCTGCAATGCAGGCCCGGCGGTTCGGTTCGTCGCCATGTCGa   >  2:828052/1‑150 (MQ=255)
                                                                                                                                                   cTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGAc  >  1:1219632/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTCGAC  >  NC_002947/5706495‑5706791

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