Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 693,590 C→G 8.2% G151A (GGC→GCC)  PP_0595 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947693,5900CG 8.2% 208.0 / 9.7 73G151A (GGC→GCC) PP_0595LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (31/36);  new base G (3/3);  total (34/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.16e-01

TCTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGC  >  NC_002947/693443‑693727
                                                                                                                                                    |                                                                                                                                           
tcTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGcc                                                                                                                                            <  5:117592/149‑1 (MQ=255)
        ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGATAAACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTc                                                                                                                                      <  1:245602/148‑1 (MQ=255)
        ctgGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGc                                                                                                                                        <  6:64000/145‑1 (MQ=255)
        agggccgtcgGCTTACAACCGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGGTCGCTGTACAAGGGTTGATACTCAAGCCCCGAT‑AACGGGCTGTTCTGAGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCt                                                                                                                                       <  3:243276/141‑1 (MQ=255)
          ggTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGAcg                                                                                                                                  <  5:183938/149‑2 (MQ=255)
          ggTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGAcc                                                                                                                                  >  5:185489/1‑149 (MQ=255)
                  gCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTggg                                                                                                                          >  3:12810/1‑149 (MQ=255)
                     tAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGc                                                                                                                       >  6:223239/1‑149 (MQ=255)
                       gAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCg                                                                                                                      <  3:186695/148‑1 (MQ=255)
                       gAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCg                                                                                                                      >  4:9371/1‑148 (MQ=255)
                       gAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCg                                                                                                                      <  3:286342/148‑1 (MQ=255)
                       gAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGa                                                                                                                     <  6:185489/149‑1 (MQ=255)
                        aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGc                                                                                                                                        <  1:269328/129‑1 (MQ=255)
                        aaCAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGc                                                                                                                                        >  2:269328/1‑129 (MQ=255)
                          cAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCa                                                                                                                  >  3:228791/1‑149 (MQ=255)
                          cAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCGCGACCTGGCTGGGGGCGATCa                                                                                                                  >  8:103916/1‑149 (MQ=255)
                           aGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGAcg                                                                                                                                  >  2:53082/1‑131 (MQ=255)
                           aGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGAcg                                                                                                                                  <  1:53082/132‑2 (MQ=255)
                           aGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑t                                                                                                                <  4:12810/149‑1 (MQ=255)
                                   gCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCt                                                                                                        >  2:64306/1‑149 (MQ=255)
                                      aaCCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                    >  5:249501/1‑119 (MQ=255)
                                      aaCCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGGCGTGCTCGa                                                                                                                                    <  6:249501/119‑1 (MQ=255)
                                          gcgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATgcg                                                                                                 >  1:284625/1‑149 (MQ=255)
                                           cgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGa                                                                                                           >  6:9723/1‑138 (MQ=255)
                                           cgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGa                                                                                                           <  5:9723/138‑1 (MQ=255)
                                           cgcAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATgcg                                                                                                 >  5:147361/1‑148 (MQ=255)
                                                   agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCTACGTGCAGGCTGCCGTCGAGAACTCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAgg                                                                                        <  5:223239/149‑1 (MQ=255)
                                                   agcagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAgg                                                                                        <  1:8124/149‑1 (MQ=255)
                                                     cagcGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGc                                                                                      >  2:265021/1‑149 (MQ=255)
                                                       gcgagcgTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGACCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGccc                                                                                    >  3:304206/1‑149 (MQ=255)
                                                        cgagcgTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGcccc                                                                                   <  7:58857/149‑1 (MQ=255)
                                                             gTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGTTCTGTGGTACGACGCCAACTTGCAGGCTGCCGTCGAGCACGCCGTGCTAGACCTTGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATcc                                                                             >  3:88971/1‑149 (MQ=255)
                                                             gTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGAACTGGATGCGCACATCACGGCCCC‑GATcc                                                                             >  3:93383/1‑149 (MQ=255)
                                                                cGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           <  5:224267/148‑1 (MQ=255)
                                                                cGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           <  5:2543/148‑1 (MQ=255)
                                                                cGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           <  6:147361/148‑1 (MQ=255)
                                                                cGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           <  1:265021/148‑1 (MQ=255)
                                                                        aaCGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGGTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GCTCCTTGAGTTACGc                                                                  >  5:206242/1‑149 (MQ=255)
                                                                              tGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           <  6:100500/134‑1 (MQ=255)
                                                                              tGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCtt                                                                           >  5:100500/1‑134 (MQ=255)
                                                                               gATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGc                                                            <  2:284625/148‑1 (MQ=255)
                                                                                    tCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGt                                                      <  3:172087/149‑1 (MQ=255)
                                                                                    tCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGt                                                      >  4:172087/1‑149 (MQ=255)
                                                                                      aaGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGa                                                    <  1:8310/149‑1 (MQ=255)
                                                                                       aGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGAt                                                   <  7:103916/149‑1 (MQ=255)
                                                                                        gCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATa                                                  <  5:117734/149‑1 (MQ=255)
                                                                                           ccGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAgg                                                              >  5:278550/1‑134 (MQ=255)
                                                                                           ccGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAgg                                                              <  6:278550/134‑1 (MQ=255)
                                                                                               t‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCAtgt                                           <  1:64306/149‑1 (MQ=255)
                                                                                                  aCGGGCTGGTCTGTGGTACGACGCCAACTTGCAGGCTGCCTTCGAGCACGCCGTTCTCGAGCTGGCTGGGGGCGATGA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGc                                         <  4:304206/149‑1 (MQ=255)
                                                                                                    gggCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCgg                                       >  4:94926/1‑149 (MQ=255)
                                                                                                     ggCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGc                                      <  6:132368/149‑1 (MQ=255)
                                                                                                           tCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTc                                >  6:3743/1‑149 (MQ=255)
                                                                                                              gtgGTACGACGCCTACGTGCTGGCTGCCGTCGTGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAc                              <  8:174978/148‑1 (MQ=255)
                                                                                                              gtgGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATAAGGGCCCC‑GAGCCTTGAGTTCCGCCACGGCATGGGTGATACGCATGTGCGGCAGGGTCAc                              >  7:174978/1‑148 (MQ=255)
                                                                                                               tgGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAc                              >  8:293014/1‑147 (MQ=255)
                                                                                                               tgGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAc                              <  5:3743/147‑1 (MQ=255)
                                                                                                                ggTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCCTGATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCAc                              >  5:269838/1‑147 (MQ=255)
                                                                                                                     gacgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                   <  8:143035/136‑1 (MQ=255)
                                                                                                                     gacgCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAgg                                   >  7:143035/1‑136 (MQ=255)
                                                                                                                         ccAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGt                  <  3:9371/149‑1 (MQ=255)
                                                                                                                          cAACGTGCAGGCTGCCGTCGCGCACGCCGTGCTCGCCCGGGGTGGGGGCGATCA‑AGCGGATATGGCTGCGGACCTCAGGGCCCC‑GGTCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGt                  >  5:198624/1‑148 (MQ=255)
                                                                                                                               tGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGcc            <  7:195995/149‑1 (MQ=255)
                                                                                                                                gCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATAc                                                 >  5:163426/1‑111 (MQ=255)
                                                                                                                                gCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATAc                                                 <  6:163426/111‑1 (MQ=255)
                                                                                                                                    gCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGt       >  2:96032/1‑149 (MQ=255)
                                                                                                                                      ccaggTCGAGCACGGCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGtt      >  2:14827/5‑148 (MQ=255)
                                                                                                                                        ccGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTc     <  1:96032/147‑1 (MQ=255)
                                                                                                                                        ccGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑AGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGGTGTCGGTCAGGCCGATGTTCAg   <  4:53444/149‑1 (MQ=255)
                                                                                                                                         cGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGc  >  3:264194/1‑149 (MQ=255)
                                                                                                                                          gTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGc  >  6:91675/1‑148 (MQ=255)
                                                                                                                                          gTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGc  <  5:91675/148‑1 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCAGTGCGGATCTGGATgcgc                                                                                                >  1:184644/1‑47 (MQ=255)
                                                                                                                                                   gCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATgcgc                                                                                                <  2:184644/46‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
TCTATCTGCTGGTCGTCGGCGTAGAACAGCGGGTGGCCAACCGCGCAGTACAGCAGCGAGCGTTCGCTGTACAACGGTTGATACTCAAGCCCCGAT‑AACGGGCTGGTCTGTGGTACGACGCCAACGTGCAGGCTGCCGTCGAGCACGCCGTGCTCGACCTGGCTGGGGGCGATCA‑TGCGGATCTGGATGCGCACATCAGGGCCCC‑GATCCTTGAGTTCCGCCAGGGCATGGGTGATACGCATGTGCGGCAGGGTCACCAGGTTGTCGGTCAGGCCGATGTTCAGC  >  NC_002947/693443‑693727

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑