Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,391,824 G→T 8.1% L24L (CTC→CTA)  PP_5709 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,391,8240GT 8.1% 210.3 / 6.8 70L24L (CTC→CTA) PP_5709hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (34/30);  new base T (4/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

GCATCAGGGAGCTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAGTATCATT  >  NC_002947/5391681‑5391969
                                                                                                                                               |                                                                                                                                                 
gCATCAGGGAGCTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGagcag                                                                                                                                              <  2:280168/149‑1 (MQ=255)
     aGGGAGCTTCGTCAGCACTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCTAGCAGCCCtt                                                                                                                                         <  4:281519/149‑1 (MQ=255)
      gggAGCTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTc                                                                                                                                        <  6:235038/149‑1 (MQ=255)
       ggAGCTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCtt                                                                                                                                         >  8:144418/1‑147 (MQ=255)
           cTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGGAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATg                                                                                                                                   >  8:70903/1‑149 (MQ=255)
               gTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCtt                                                                                                                               >  3:40188/1‑149 (MQ=255)
                    cccTTCGTTCTGCTTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATGCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTTGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTg                                                                                                                            <  1:193666/147‑1 (MQ=255)
                    cacTTCGTTCTGCGTAAACTTGTAAGATAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTAAGATGccaattt                                                                                                                            >  2:193666/3‑142 (MQ=255)
                     ccTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGAATAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGAt                                                                                                                         >  4:289673/1‑149 (MQ=255)
                          gTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTc                                                                                                                     <  1:140809/148‑1 (MQ=255)
                           ttCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGc                                                                                                                   <  1:321787/149‑1 (MQ=255)
                                 gTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCTCGCAGccgtctc                                                                                                                                       <  6:100351/123‑6 (MQ=255)
                                 gTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCTCGCAGccgtctc                                                                                                                                       >  5:100351/1‑118 (MQ=255)
                                 gTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGt                                                                                                             >  6:288107/1‑149 (MQ=255)
                                      cTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCa                                                                                                        <  3:289673/149‑1 (MQ=255)
                                         gTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAg                                                                                                       <  3:51034/147‑1 (MQ=255)
                                          tAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGt                                                                                                     >  6:123357/1‑148 (MQ=255)
                                            aGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTCTCGCCTCGGTCTTCAGGtcaa                                                                                                  >  8:221298/1‑149 (MQ=255)
                                               tAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAgg                                                                                                      <  1:217265/142‑1 (MQ=255)
                                               tAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGtcaatc                                                                                                >  5:137833/1‑148 (MQ=255)
                                               tAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTAAgg                                                                                                      >  2:217265/1‑142 (MQ=255)
                                                    tGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAAtct                                                                                          >  4:151326/1‑149 (MQ=255)
                                                      gCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAAtcta                                                                                         >  5:242960/1‑147 (MQ=255)
                                                      gCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAAtc                                                                                           <  6:242960/146‑1 (MQ=255)
                                                      gCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTcg                                                                                        <  6:137833/149‑1 (MQ=255)
                                                      gCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTAAATCAAAtatct                                                                                        >  4:51063/1‑148 (MQ=255)
                                                      gCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGACAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTcg                                                                                        >  5:76757/1‑149 (MQ=255)
                                                                  tCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGGAGGTCGACTGCAACTGGGAGACGGCTGCGAGCACCCCTTCAGATGCCTTGTGTATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGct                                                                             >  2:270316/1‑147 (MQ=255)
                                                                    aCTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAg                                                                          >  5:58216/1‑149 (MQ=255)
                                                                        tCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCtt                                                                      >  6:326494/1‑149 (MQ=255)
                                                                        tCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCtt                                                                      <  5:326494/149‑1 (MQ=255)
                                                                         cACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTg                                                                     >  4:24201/1‑149 (MQ=255)
                                                                          aCCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGc                                                                    >  3:161576/1‑149 (MQ=255)
                                                                           ccGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGcc                                                                   <  6:93684/149‑1 (MQ=255)
                                                                           ccGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTAGCGATGCTAGCATAGCATTTcc                                                                   >  8:298601/1‑149 (MQ=255)
                                                                                  aaCGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGAc                                                            >  5:119215/1‑149 (MQ=255)
                                                                                  aaCGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGAc                                                            >  8:180215/1‑149 (MQ=255)
                                                                                     ggCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGagcagc                                                                                                                                             <  4:199505/65‑1 (MQ=255)
                                                                                     ggCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGagcagc                                                                                                                                             >  3:199505/1‑65 (MQ=255)
                                                                                     ggCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCgg                                                         <  7:144418/149‑1 (MQ=255)
                                                                                         aTCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTAc                                                     <  6:76757/149‑1 (MQ=255)
                                                                                            cAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATc                                                  <  6:58216/149‑1 (MQ=255)
                                                                                             agagTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCa                                                 <  5:288107/149‑1 (MQ=255)
                                                                                               agTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAg                                                >  5:208836/1‑148 (MQ=255)
                                                                                               agTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTCTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGa                                               <  7:221298/149‑1 (MQ=255)
                                                                                                    tcGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCtt                                          <  6:217197/148‑1 (MQ=255)
                                                                                                      ggTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGt                                        >  2:299476/1‑149 (MQ=255)
                                                                                                       gTTTGCCAGTGTATGTCGACTGCAACTGTGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGTTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGtt                                       <  7:298601/149‑1 (MQ=255)
                                                                                                        tttGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGttt                                      <  6:208836/149‑1 (MQ=255)
                                                                                                        tttGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGttt                                      >  8:253631/1‑149 (MQ=255)
                                                                                                          tGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGtttt                                     <  3:151326/148‑1 (MQ=255)
                                                                                                           gCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTc                                    <  7:180215/148‑1 (MQ=255)
                                                                                                            ccAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGc                                                                    <  1:42499/115‑1 (MQ=255)
                                                                                                            ccAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGc                                                                    >  2:42499/1‑115 (MQ=255)
                                                                                                               gggTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTc                               <  3:140337/149‑1 (MQ=255)
                                                                                                                        gACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCtt                      >  3:139605/1‑149 (MQ=255)
                                                                                                                         aCTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTg                     <  7:70903/149‑1 (MQ=255)
                                                                                                                             cAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCAt                 >  8:295606/1‑149 (MQ=255)
                                                                                                                                 tGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGAttt                                                                                                                       <  6:89381/43‑1 (MQ=255)
                                                                                                                                 tGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGAttt                                                                                                                       >  5:89381/1‑43 (MQ=255)
                                                                                                                                   ggAGACGGCTGCGAACAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAg            >  4:81628/1‑148 (MQ=255)
                                                                                                                                   ctgaagggCTGCTCGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAg            >  1:58119/7‑148 (MQ=255)
                                                                                                                                   ctgaagcgctgctCGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAg            >  7:109702/8‑148 (MQ=255)
                                                                                                                                    tgaagggCTGCTCGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAg            >  2:23604/6‑147 (MQ=255)
                                                                                                                                    gagaCGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAg            >  8:248113/1‑147 (MQ=255)
                                                                                                                                     gaagggCTGCTCGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATTAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAg         >  2:297398/5‑149 (MQ=255)
                                                                                                                                      gaCGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAg         >  6:317479/1‑148 (MQ=255)
                                                                                                                                         ggCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAGTATc     <  2:93337/149‑1 (MQ=255)
                                                                                                                                          gCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAGTATCa    <  5:123357/149‑1 (MQ=255)
                                                                                                                                            tGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGCTCTTCAGGTCAAGCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTATTCACTTCAATAACCTTGCCATGCAAGCAGTATCAtt  >  7:50644/1‑149 (MQ=255)
                                                                                                                                              cGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAGTATCAtt  <  1:299476/147‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GCATCAGGGAGCTTCGTCAGCCCTTCGTTCTGCGTCAACTTGTAAGCTAGGCTGGCAACCACTTCCTCACTATCACCGACTAAACGGCCATCCAGAGTGAACGGTTTGCCAGGGTAGGTCGACTGCAACTGGGAGACGGCTGCGAGCAGCCCTTCAGATGCCTTTTGGATTTTCGCCTCGGTCTTCAGGTCAATCAAATCTCGCGATGCTCGCATAGCCTTGCCCTATGACCGGCTACATCAGAATCTTGTTTTCACTTCAATAACCTTGCCATGCAAGCAGTATCATT  >  NC_002947/5391681‑5391969

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