Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,890,847 A→C 7.2% intergenic (+79/+41) PP_5164 → / ← plpB penicillin amidase family protein/NLPA lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,890,8470AC 7.2% 191.2 / 7.1 69intergenic (+79/+41)PP_5164/plpBpenicillin amidase family protein/NLPA lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base C (2/3);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGCAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTG  >  NC_002947/5890698‑5890996
                                                                                                                                                     |                                                                                                                                                     
ggCAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGAcaca                                                                                                                                                       >  1:4025/1‑150 (MQ=255)
  cAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                     <  1:1176579/150‑1 (MQ=255)
  cAGTATCAGCGGGCGCAGATGGGCGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAg                                                                                                                                                     <  2:873254/150‑1 (MQ=255)
     tATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCGTCTACGGTGCTGTGCGGGCAGGGCCGAACACAACCGCCGACACAAgtgg                                                                                                                                                  >  2:1079935/1‑147 (MQ=255)
      aTCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGgtgt                                                                                                                                                 <  2:1463661/150‑1 (MQ=255)
            gggCGCAGATGGGGGTGATTCCGGCGCAAAGTCCCATAGGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGt                                                                                                                                           <  2:1610403/150‑1 (MQ=255)
            gggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGt                                                                                                                                           <  2:1631028/150‑1 (MQ=255)
             ggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                          <  2:1655124/150‑1 (MQ=255)
             ggCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                          <  1:599714/150‑1 (MQ=255)
              gcgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTg                                                                                                                                         <  2:1074135/150‑1 (MQ=255)
               cgcAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTgg                                                                                                                                        <  1:1771672/150‑1 (MQ=255)
                 cAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGca                                                                                                                                      <  1:617129/150‑1 (MQ=255)
                       ggggTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATc                                                                                                                                <  2:1755106/150‑1 (MQ=255)
                       ggggTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATc                                                                                                                                >  2:1762750/1‑150 (MQ=255)
                          gTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCacg                                                                                                                             <  2:566026/150‑3 (MQ=255)
                          gTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAgg                                                                                                                             <  1:1306642/150‑1 (MQ=255)
                           tGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGc                                                                                                                            <  2:431707/150‑1 (MQ=255)
                                  ggCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTa                                                                                                                     <  2:976506/150‑1 (MQ=255)
                                   gcgcAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAg                                                                                                                    <  2:1466961/150‑1 (MQ=255)
                                      cAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGa                                                                                                                 >  2:1272564/1‑150 (MQ=255)
                                       aaaGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGttttt                                                                                                                                       <  2:43969/127‑4 (MQ=255)
                                       aaaGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGttttt                                                                                                                                       >  1:43969/1‑124 (MQ=255)
                                           tACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCa                                                                                                            <  1:467438/150‑1 (MQ=255)
                                              cTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                         >  1:586562/1‑150 (MQ=255)
                                              cTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGaa                                                                                                         >  2:448849/1‑150 (MQ=255)
                                                   ggCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGttt                                                                                                                                         >  2:1725478/1‑112 (MQ=255)
                                                   ggCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACCTTGTGTCGGCGttt                                                                                                                                         <  1:1725478/113‑2 (MQ=255)
                                                        gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                                   <  2:1528391/130‑1 (MQ=255)
                                                        gTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGc                                                                                                                   >  1:1528391/1‑130 (MQ=255)
                                                         tACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTa                                                                                              <  2:123609/150‑1 (MQ=255)
                                                             cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                          <  2:1105479/102‑1 (MQ=255)
                                                             cccGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGtt                                                                                                                                          >  1:1105479/1‑102 (MQ=255)
                                                                     aaGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGc                                                                                                                                             <  2:115620/91‑1 (MQ=255)
                                                                     aaGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGc                                                                                                                                             >  1:115620/1‑91 (MQ=255)
                                                                         ggCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAg                                                                                                                              <  1:1071381/102‑1 (MQ=255)
                                                                         ggCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAg                                                                                                                              >  2:1071381/1‑102 (MQ=255)
                                                                             gctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGc                                                                                                                                       >  2:1266129/1‑89 (MQ=255)
                                                                             gctTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGc                                                                                                                                       <  1:1266129/89‑1 (MQ=255)
                                                                                   cGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCaatg                                                                    >  1:205994/1‑150 (MQ=255)
                                                                                          ttCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCt                                                             >  1:409323/1‑150 (MQ=255)
                                                                                             gcagcaAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGACGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGa                                                          >  1:1852524/1‑150 (MQ=255)
                                                                                                 caAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGt                                                                                               >  2:373517/1‑109 (MQ=255)
                                                                                                 caAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGt                                                                                               <  1:373517/109‑1 (MQ=255)
                                                                                                   aaaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCg                                                    >  1:809281/1‑150 (MQ=255)
                                                                                                    aaGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCgg                                                   >  1:1833901/1‑150 (MQ=255)
                                                                                                     aGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGGGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCggg                                                  <  1:1099713/150‑1 (MQ=255)
                                                                                                        cTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTg                                               <  2:1833901/150‑1 (MQ=255)
                                                                                                         ttCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTTGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTgg                                              <  2:4025/150‑1 (MQ=255)
                                                                                                              cGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAgg                                         <  2:1852524/150‑1 (MQ=255)
                                                                                                                    tgtgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACGTGTGGGAAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTg                                   >  1:1583303/1‑150 (MQ=255)
                                                                                                                      tgCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGc                                 >  1:403682/1‑150 (MQ=255)
                                                                                                                          ggCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTAGTCAGGGCCTTGGCCAGt                             >  1:1161071/1‑150 (MQ=255)
                                                                                                                             aGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTc                          >  2:1093141/1‑150 (MQ=255)
                                                                                                                             aGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTc                          >  1:1330492/1‑150 (MQ=255)
                                                                                                                              gggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCt                         <  1:1272564/150‑1 (MQ=255)
                                                                                                                               ggCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTg                        >  2:1533926/1‑150 (MQ=255)
                                                                                                                                  tGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGAt                     >  1:1747427/1‑150 (MQ=255)
                                                                                                                                    aagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTTGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGtt                            <  1:555356/141‑1 (MQ=255)
                                                                                                                                    aagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTTGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGtt                            >  2:555356/1‑141 (MQ=255)
                                                                                                                                    aagaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATgg                   >  1:844760/1‑150 (MQ=255)
                                                                                                                                     agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc                  <  1:448849/150‑1 (MQ=255)
                                                                                                                                     agaaAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGc                  <  1:1093141/150‑1 (MQ=255)
                                                                                                                                        aaaaCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCtt               <  2:1330492/150‑2 (MQ=255)
                                                                                                                                           aCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCg                                                                                                                  >  2:1034992/1‑48 (MQ=255)
                                                                                                                                           aCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCg                                                                                                                  <  1:1034992/48‑1 (MQ=255)
                                                                                                                                           aCGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACt            >  1:637891/1‑150 (MQ=255)
                                                                                                                                            cGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTAGACTTCTGGGCTGGTCAGGGCCTTGGCCCGTTTCTGGATGGCCTCACTg           >  2:707183/1‑150 (MQ=255)
                                                                                                                                               cGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACt            <  2:1903497/146‑1 (MQ=255)
                                                                                                                                               cGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACt            >  1:1903497/1‑146 (MQ=255)
                                                                                                                                                 acacAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCtt      <  2:403682/150‑1 (MQ=255)
                                                                                                                                                  cacaAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCttg     <  2:1747427/150‑1 (MQ=255)
                                                                                                                                                     aaGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCttgttg  >  1:1276777/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCAGTATCAGCGGGCGCAGATGGGGGTGATTCCGGCGCAAAGTACCTTACGGCTGGTACCCCCGAAGTAAGAGGCTGCTTTGCGGCCGCTTCGCAGCAAAAGGCTTCTACGGTGCTGTGCGGGCAGGGCTGAAGAAAAACGCCGACACAAGGTGTCGGCGTTGGCACATATCAGGCGTTGCTAGCGAAATCAGAACGCCGGCAGTACGGCACCTTGGTACTTCTTGGCAATGAATGCCTTGACTTCCGGGCTGGTCAGGGCCTTGGCCAGTTTCTGGATGGCCTCACTGTCCTTGTTG  >  NC_002947/5890698‑5890996

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