Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,671,822 A→G 6.0% I51V (ATA→GTA)  PP_4133 → transcriptional regulator‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,671,8220AG 6.0% 193.4 / 2.5 66I51V (ATA→GTA) PP_4133transcriptional regulator‑like protein
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base G (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACTCACATACTTGCTGAGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGCCAGGCCATCCT  >  NC_002947/4671678‑4671970
                                                                                                                                                |                                                                                                                                                    
aCTCACATACTTGCTGAGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCa                                                                                                                                                      >  5:264677/1‑145 (MQ=255)
    acaTACTTGCTGAGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATc                                                                                                                                              >  3:159008/1‑149 (MQ=255)
     caTACTTGCTGAGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCg                                                                                                                                             <  6:264677/149‑1 (MQ=255)
                aGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAgcat                                                                                                                                  <  4:159008/149‑1 (MQ=255)
                aGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAgcat                                                                                                                                  >  1:214134/1‑149 (MQ=255)
                  gTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAg                                                                                                                                                   >  4:191778/1‑130 (MQ=255)
                  gTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAg                                                                                                                                                   <  3:191778/130‑1 (MQ=255)
                        aCCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCGTAGATATCGCCAAAAAGCATGCATGCt                                                                                                                           >  6:66870/1‑148 (MQ=255)
                        aCCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCGTAGATATCGCCAAAAAGCATGCATGCt                                                                                                                           >  2:127990/1‑148 (MQ=255)
                          cctccTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCac                                                                                                                        <  4:166587/146‑1 (MQ=255)
                               ttAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGa                                                                                                                                                  >  6:322045/1‑118 (MQ=255)
                               ttAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGa                                                                                                                                                  <  5:322045/118‑1 (MQ=255)
                                 aaGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCactggactg                                                                                                                 >  3:279811/1‑149 (MQ=255)
                                   gCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGc                                                                                                               <  2:214134/149‑1 (MQ=255)
                                       cttACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTa                                                                                                             <  4:279811/147‑1 (MQ=255)
                                           cTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCC‑AAAAGCATGCATGCTCACTGGACTGGCTATTGCTgg                                                                                                      <  6:137863/149‑1 (MQ=255)
                                            tGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGTAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGc                                                                                                               <  1:30391/140‑1 (MQ=255)
                                            tGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGc                                                                                                               >  2:30391/1‑140 (MQ=255)
                                             ggCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTggg                                                                                                     >  1:166284/1‑149 (MQ=255)
                                              gcgcAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTgg                                                                                                      <  1:18712/147‑1 (MQ=255)
                                              gcgcAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGc                                                                                                    >  7:10008/1‑149 (MQ=255)
                                                 cAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCGTAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCaa                                                                                                 <  5:66870/149‑1 (MQ=255)
                                                 cAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCGTAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCaa                                                                                                 <  1:127990/149‑1 (MQ=255)
                                                  aaaCACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCactggac                                                                                                                   >  2:107357/1‑130 (MQ=255)
                                                  aaaCACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCactggac                                                                                                                   <  1:107357/130‑1 (MQ=255)
                                                      acTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGa                                                                                            <  8:282235/149‑1 (MQ=255)
                                                        tGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTgg                                                                                                      <  1:9269/137‑1 (MQ=255)
                                                        tGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTgg                                                                                                      >  2:9269/1‑137 (MQ=255)
                                                         gATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCAAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGa                                                                                            >  8:136191/1‑146 (MQ=255)
                                                         gATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCAAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGa                                                                                            <  7:136191/146‑1 (MQ=255)
                                                           tGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATACATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGcac                                                                                        >  8:102508/1‑148 (MQ=255)
                                                                  cTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAAcgccgc                                                                                >  5:69320/1‑149 (MQ=255)
                                                                    gTCACGTGCAGTGTCGGTCAGCCCGCAAACGATTAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGc                                                                                                    <  6:338534/127‑1 (MQ=255)
                                                                    gTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTATCCATAGATATCGCCAAAAAGCATGCATGCTAACTGGACTCGCTATTGATGGGc                                                                                                    >  5:338534/1‑127 (MQ=255)
                                                                     tCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACgccgcc                                                                               <  8:10008/147‑1 (MQ=255)
                                                                       aCGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCaa                                                                           <  2:166284/149‑1 (MQ=255)
                                                                       aCGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCaa                                                                           <  2:292639/149‑1 (MQ=255)
                                                                         gTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTAGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGt                                                                         >  3:298848/1‑149 (MQ=255)
                                                                            cACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGa                                                                       <  6:69320/148‑1 (MQ=255)
                                                                               tGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAg                                                                                           >  3:97629/1‑125 (MQ=255)
                                                                               tGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAg                                                                                           <  4:97629/125‑1 (MQ=255)
                                                                                    gTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTg                                                                        <  3:284022/139‑1 (MQ=255)
                                                                                    gTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTg                                                                        >  4:284022/1‑139 (MQ=255)
                                                                                    gTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTg                                                              >  2:203412/1‑149 (MQ=255)
                                                                                            gCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGc                                                       >  3:284469/1‑148 (MQ=255)
                                                                                            gCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGc                                                       >  3:278111/1‑148 (MQ=255)
                                                                                            gCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGc                                                       <  4:284469/148‑1 (MQ=255)
                                                                                                   cTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCt                                               <  4:298848/149‑1 (MQ=255)
                                                                                                        cagcTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATACATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAAc                                          <  7:102508/149‑1 (MQ=255)
                                                                                                                   ttCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACg                               >  5:207503/1‑149 (MQ=255)
                                                                                                                   ttCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACg                               >  5:186478/1‑149 (MQ=255)
                                                                                                                      gcgAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGctc                            <  4:278111/149‑1 (MQ=255)
                                                                                                                              aTTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAg                    >  6:174458/1‑149 (MQ=255)
                                                                                                                               ttCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGc                                                <  1:295747/120‑1 (MQ=255)
                                                                                                                               ttCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGc                                                >  2:295747/1‑120 (MQ=255)
                                                                                                                                 ccTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCAc                                <  7:104429/134‑1 (MQ=255)
                                                                                                                                 ccTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCAc                                >  8:104429/1‑134 (MQ=255)
                                                                                                                                 ccTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGa                 >  2:6456/1‑149 (MQ=255)
                                                                                                                                   ttACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGAtg                                                                     >  6:224546/1‑95 (MQ=255)
                                                                                                                                   ttACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGAtg                                                                     <  5:224546/95‑1 (MQ=255)
                                                                                                                                     aCTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATc               <  1:6456/147‑1 (MQ=255)
                                                                                                                                      cTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGcc            <  5:174458/149‑1 (MQ=255)
                                                                                                                                         aCTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAg                    <  6:52803/138‑1 (MQ=255)
                                                                                                                                         aCTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAg                    >  5:52803/1‑138 (MQ=255)
                                                                                                                                         aCTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGCCAgg         <  6:207503/149‑1 (MQ=255)
                                                                                                                                           tttGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGCCAGGc        >  5:80157/1‑148 (MQ=255)
                                                                                                                                               cATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGCCAGGCCATc    <  6:186478/148‑1 (MQ=255)
                                                                                                                                                aTAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCGCGCTCTCGCACAGCGATCGCCAGGCCATCct  >  8:154011/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACTCACATACTTGCTGAGGTGCTCACCCGCCTTAAGCTTCTTACTGGCGCAAACACTGATGCCGAGCTGTCACGTGCACTGTCGGTCAGCCCGCAAACGCTGAGCAGCTGGAAGGTTCGCGAAAGCATTCCTTACTCACTTTGCATAGATATCGCCAAAAAGCATGCATGCTCACTGGACTGGCTATTGCTGGGCCAAGCCGAGCACAACGCCGCCCCAAGTGATGTGGGTTGGGAGTGCGACATGCTCGAACGGCTGCGCACGCTCTCGCACAGCGATCGCCAGGCCATCCT  >  NC_002947/4671678‑4671970

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: ‑