Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,672,522 C→G 6.0% G76A (GGT→GCT)  PP_4135 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,672,5220CG 6.0% 213.7 / 3.0 67G76A (GGT→GCT) PP_4135hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/27);  new base G (2/2);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CAGGTTGTCAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCAGGTGCTCCC  >  NC_002947/4672379‑4672669
                                                                                                                                               |                                                                                                                                                   
cAGGTTGTCAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGCAAGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCCGc                                                                                                                                                 >  8:208213/1‑148 (MQ=255)
cAGGTTGTCAAAAACTGATAGATCAGGCATGCACTTTCCCTCCGCTCGTAACAGTCATACAGGCGGTTGAACCCAGGCTATACATGACCGCCGATTCGACAAAACGGATGTAGAACACCCGCAAGGTCCGCCAATGATCATGACCGGc                                                                                                                                                 >  8:202844/1‑148 (MQ=255)
    ttGTCAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGGTCATGATCATGGGCGGcccttg                                                                                                                                            <  4:95630/149‑5 (MQ=255)
      gTCAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATcc                                                                                                                                          >  3:315167/1‑149 (MQ=255)
        cAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTc                                                                                                                                        >  5:96016/1‑149 (MQ=255)
         aaaaaCTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTc                                                                                                                                        >  7:297620/1‑148 (MQ=255)
            aaCTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAgg                                                                                                                                    <  1:9180/149‑1 (MQ=255)
             aCTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGgc                                                                                                                                   >  1:234478/1‑149 (MQ=255)
                 aTAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCtt                                                                                                                               >  5:163934/1‑149 (MQ=255)
                         ggCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGcccc                                                                                                                       <  4:315167/149‑1 (MQ=255)
                               cACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCtctc                                                                                                                 >  4:181653/1‑149 (MQ=255)
                                  ttGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGcc                                                                                                                         >  4:264938/1‑138 (MQ=255)
                                  ttGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGcc                                                                                                                         <  3:264938/138‑1 (MQ=255)
                                     cccTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTa                                                                                                           >  4:223443/1‑149 (MQ=255)
                                        tCCGCGCGGAACAGGCAGACAGGCGGTTGAACACAAGCTATACATTACCGACGCTTCTACAAAACGCATGTACCACACCCGCAAGGTCCGCCCAACATCATGACAGGCCACATCCTAGAGCCACATGTAGCCCCGCTCTCGCAGTACt                                                                                                         >  2:234626/1‑148 (MQ=255)
                                           gcgcgGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGAcc                                                                                                     >  7:132636/1‑149 (MQ=255)
                                             gcgGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGc                                                                                                   <  6:163934/149‑1 (MQ=255)
                                                   cAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCgg                                                                                             <  8:297620/149‑1 (MQ=255)
                                                      gCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACAACCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGc                                                                                           >  8:102841/1‑148 (MQ=255)
                                                         gACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAAt                                                                                       <  2:234478/149‑1 (MQ=255)
                                                             ggCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTcgc                                                                                   <  6:96016/149‑1 (MQ=255)
                                                              gCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTcgcg                                                                                  <  8:132636/149‑1 (MQ=255)
                                                                      aCCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCt                                                                          >  7:138091/1‑149 (MQ=255)
                                                                          aGGCAATACATGACCGCCGATTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCAAGATCATGACCGGCCACATCCTCGAGGAGCTTGTAGCCCCGCTCTCGCACTACTGACCGCCAGCGGGCGAATTACCCTGCTGCCTCGAt                                                                      >  4:212917/1‑149 (MQ=255)
                                                                            gCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATgg                                                                    <  5:135789/149‑1 (MQ=255)
                                                                                  cATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAgcg                                                               >  8:166404/1‑148 (MQ=255)
                                                                                       ccgccgCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCa                                                                  <  2:91843/140‑1 (MQ=255)
                                                                                       ccgccgCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCa                                                                  >  1:91843/1‑140 (MQ=255)
                                                                                        cgccgcTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAAcgcg                                                        >  6:59251/1‑149 (MQ=255)
                                                                                           cgcTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCgg                                                     <  8:138091/149‑1 (MQ=255)
                                                                                             cTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGc                                                    >  3:333474/1‑148 (MQ=255)
                                                                                                  aCAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTc                                              >  4:261532/1‑149 (MQ=255)
                                                                                                  aCAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTc                                              <  3:223443/149‑1 (MQ=255)
                                                                                                    aaaaCGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTc                                                                                     <  1:238323/108‑1 (MQ=255)
                                                                                                    aaaaCGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTc                                                                                     >  2:238323/1‑108 (MQ=255)
                                                                                                    aaaaCGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTcc                                             >  4:2838/1‑148 (MQ=255)
                                                                                                      aaCGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAAt                                          >  8:306500/1‑149 (MQ=255)
                                                                                                      aaCGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCACTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAAt                                          >  7:228470/1‑149 (MQ=255)
                                                                                                         ggATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCATGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGAAGAGCAACGAGCTGCAACTTCCAATg                                         >  1:97043/1‑147 (MQ=255)
                                                                                                            tGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCt                                    >  4:328278/1‑149 (MQ=255)
                                                                                                              tAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCa                                                  <  5:46651/133‑1 (MQ=255)
                                                                                                              tAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCa                                                  >  6:46651/1‑133 (MQ=255)
                                                                                                              tAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCtt                                   <  3:328278/148‑1 (MQ=255)
                                                                                                                gCACACCCGCACCGTCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCa                                 <  3:212917/148‑1 (MQ=255)
                                                                                                                gCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCa                                 <  4:333474/148‑1 (MQ=255)
                                                                                                                gCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCAc                                <  3:181653/149‑1 (MQ=255)
                                                                                                                  acacCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTCCTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACt                               >  2:249914/1‑148 (MQ=255)
                                                                                                                  acacCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACt                               >  2:46220/1‑148 (MQ=255)
                                                                                                                   cacCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGa                             <  5:59251/149‑1 (MQ=255)
                                                                                                                    acCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCg                                                                                                                <  3:11768/65‑1 (MQ=255)
                                                                                                                    acCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCg                                                                                                                >  4:11768/1‑65 (MQ=255)
                                                                                                                    acCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGa                             >  4:314064/1‑148 (MQ=255)
                                                                                                                     cccGCAAGGGCCTCCCATGATCATGACAGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCTAATTCGCGTGCTGCCTCGATGGCAGCGCAACGAGCGGCCACTTCCAATGCAGCTTCACTGACa                           <  7:202844/149‑1 (MQ=255)
                                                                                                                     cccGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACa                           <  7:208213/149‑1 (MQ=255)
                                                                                                                      ccGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTcgc                                                                                   >  3:84347/1‑92 (MQ=255)
                                                                                                                      ccGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTcgc                                                                                   <  4:84347/92‑1 (MQ=255)
                                                                                                                              gccggtcATGATCATGGGCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCgg                   >  8:312432/7‑148 (MQ=255)
                                                                                                                              gccggtcATGATCATGGGCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCgg                   >  4:206531/7‑148 (MQ=255)
                                                                                                                               ccgccCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCGGCTTCACTGACAGATCGCgg                   >  7:148118/1‑147 (MQ=255)
                                                                                                                                   ccATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTg             >  1:207036/1‑149 (MQ=255)
                                                                                                                                   ccATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTg             >  5:265555/1‑149 (MQ=255)
                                                                                                                                    cATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGc            <  6:265555/149‑1 (MQ=255)
                                                                                                                                     aTGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCa           <  3:314064/149‑1 (MQ=255)
                                                                                                                                      tGATCATGGGCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCAg          <  3:206531/149‑1 (MQ=255)
                                                                                                                                      tGATCATGACCGGCCACAACCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCAg          <  7:102841/149‑1 (MQ=255)
                                                                                                                                             gACCGGCCACATCCTCGAGGCGCTTGGAGCCCCGCTCTCGGAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAGCGCGCGGCCGCTTCCATGGGGGCTTCGCTGACAGCGCGCGGCAGTTGCAGGTGCTcc   >  4:66842/1‑149 (MQ=255)
                                                                                                                                              aCCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCAGGTGCTccc  <  1:46220/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CAGGTTGTCAAAAACTGATAGCTCAGGCATGCACTTGCCCTCCGCGCGGAACAGGCAGACAGGCGGTTGAACCCAGGCTATACATGACCGCCGCTTCGACAAAACGGATGTAGCACACCCGCAAGGGCCGCCCATGATCATGACCGGCCACATCCTCGAGGCGCTTGTAGCCCCGCTCTCGCAGTACTGACCGCCAGCGGGCGAATTCGCGTGCTGCCTCGATGGCAGCGCAACGCGCGGCCACTTCCAATGCAGCTTCACTGACAGATCGCGGCAGTTGCAGGTGCTCCC  >  NC_002947/4672379‑4672669

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