Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,168,438 T→A 6.6% intergenic (+338/+350) PP_1921 → / ← PP_1922 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,168,4380TA 6.6% 212.5 / 3.0 71intergenic (+338/+350)PP_1921/PP_1922hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/32);  new base A (2/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTGCAACTGAGGC  >  NC_002947/2168293‑2168585
                                                                                                                                                 |                                                                                                                                                   
gCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCg                                                                                                                                                  <  6:70739/149‑1 (MQ=255)
  aCCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGc                                                                                                                                                >  3:10975/1‑149 (MQ=255)
   ccGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGcc                                                                                                                                               <  8:27083/149‑1 (MQ=255)
   ccGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGcc                                                                                                                                               >  2:230592/1‑149 (MQ=255)
      ctgctCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGAATCGGcccag                                                                                                                                            <  5:152337/149‑4 (MQ=255)
        gctCACAAGGCTGACTTCAGTCCTGTTCCGCCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATg                                                                                                                                          >  4:66218/1‑149 (MQ=255)
         ctCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGAAAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATgg                                                                                                                                         <  4:70273/149‑1 (MQ=255)
           cacaAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGAATCGGcccag                                                                                                                                            >  6:53160/1‑141 (MQ=255)
           cacaAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGAATCGGcccag                                                                                                                                            <  5:53160/144‑4 (MQ=255)
            acaAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCg                                                                                                                                      >  6:302329/1‑149 (MQ=255)
                     aCTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGTACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCgg                                                                                                                               <  5:139188/147‑1 (MQ=255)
                        tCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGAATCGGcccaggg                                                                                                                                          >  4:47285/1‑128 (MQ=255)
                        tCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGAATCGGcccaggg                                                                                                                                          <  3:47285/133‑6 (MQ=255)
                          aGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAgg                                                                                                                        >  3:183130/1‑149 (MQ=255)
                          aGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAgg                                                                                                                        >  7:103500/1‑149 (MQ=255)
                          aGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAgg                                                                                                                        >  8:292923/1‑149 (MQ=255)
                           gTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGc                                                                                                                       >  5:340246/1‑149 (MQ=255)
                               ttttCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGa                                                                                                                          <  7:73760/142‑1 (MQ=255)
                               ttttCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGa                                                                                                                          >  8:73760/1‑142 (MQ=255)
                                     gTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTgg                                                                                                             >  1:37396/1‑149 (MQ=255)
                                     gTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTgg                                                                                                             >  7:69599/1‑149 (MQ=255)
                                     gTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTgg                                                                                                             >  5:122993/1‑149 (MQ=255)
                                     gTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAAAGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTg                                                                                                              >  3:276255/1‑148 (MQ=255)
                                      tCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCTACTgg                                                                                                             >  6:144725/1‑148 (MQ=255)
                                      tCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTgg                                                                                                             >  2:326444/1‑148 (MQ=255)
                                               cAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGa                                                                                                                          <  8:317270/126‑1 (MQ=255)
                                               cAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGa                                                                                                                          >  7:317270/1‑126 (MQ=255)
                                                    tAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCAtt                                                                                                         >  8:51171/1‑138 (MQ=255)
                                                    tAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCAtt                                                                                                         <  7:51171/138‑1 (MQ=255)
                                                      aTGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAgg                                                                                            <  4:10975/149‑1 (MQ=255)
                                                      aTGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAgg                                                                                            <  1:326444/149‑1 (MQ=255)
                                                      aTGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAgg                                                                                            <  3:66218/149‑1 (MQ=255)
                                                              cAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTAc                                                                                    <  7:115591/149‑1 (MQ=255)
                                                                 ttCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATTGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGaa                                                                                 <  1:115871/149‑1 (MQ=255)
                                                                 ttCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACg                                                                                   >  2:115871/1‑147 (MQ=255)
                                                                    cGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTc                                                                              >  5:277612/1‑149 (MQ=255)
                                                                    cGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTc                                                                              <  2:325766/149‑1 (MQ=255)
                                                                      ttGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTccc                                                                            <  6:340246/149‑1 (MQ=255)
                                                                         cTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGa                                                                          <  1:230592/148‑1 (MQ=255)
                                                                             aTGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCTACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAAc                                                                     <  5:144725/149‑1 (MQ=255)
                                                                               ggCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGGAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAAcc                                                                    <  4:183130/148‑1 (MQ=255)
                                                                                         ttCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTcc                                                                             >  8:42803/1‑129 (MQ=255)
                                                                                         ttCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTcc                                                                             <  7:42803/129‑1 (MQ=255)
                                                                                         ttCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGtctc                                                          <  4:276255/148‑1 (MQ=255)
                                                                                           cAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGcc                                                       >  3:79470/1‑149 (MQ=255)
                                                                                           cAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGcc                                                       >  7:234319/1‑149 (MQ=255)
                                                                                            aaGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACg                                                                                                    <  3:185925/103‑1 (MQ=255)
                                                                                            aaGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACg                                                                                                    >  4:185925/1‑103 (MQ=255)
                                                                                              gTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTcc                                                    >  8:185831/1‑149 (MQ=255)
                                                                                                    gTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCt                                              <  8:69599/149‑1 (MQ=255)
                                                                                                    gTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCt                                              <  2:37396/149‑1 (MQ=255)
                                                                                                         ggTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCATTCGTCTCGCCTCCAGCGCTGAGAt                                         <  5:302329/149‑1 (MQ=255)
                                                                                                          gTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGAtt                                        <  6:122993/149‑1 (MQ=255)
                                                                                                              aCGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGt                                    <  1:50930/149‑1 (MQ=255)
                                                                                                                  ggACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGa                                 >  3:275887/1‑148 (MQ=255)
                                                                                                                     cagcagGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCa                             >  5:187114/1‑149 (MQ=255)
                                                                                                                       gcagGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATg                           <  6:187114/149‑1 (MQ=255)
                                                                                                                       gcagGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATg                           <  8:103500/149‑1 (MQ=255)
                                                                                                                            cTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTtgct                      <  6:277612/149‑1 (MQ=255)
                                                                                                                             tgtgTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCCCCTCCAGCGCTGAGATTTTGTAGACTCATGTtgctg                     >  8:48859/1‑149 (MQ=255)
                                                                                                                                  aCAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGa                                                                          <  7:48659/91‑1 (MQ=255)
                                                                                                                                  aCAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGa                                                                          >  8:48659/1‑91 (MQ=255)
                                                                                                                                    aaCCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCAGACAAACCAAATCGTCTCGCATCCAGCGCTGAGATTTTGTAGACTCATGTTGCTCCGTGGAt              >  6:163062/1‑149 (MQ=255)
                                                                                                                                      cccTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTg            >  4:139771/1‑149 (MQ=255)
                                                                                                                                      cccTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTg            >  4:134579/1‑149 (MQ=255)
                                                                                                                                             ccGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATg                           <  4:283602/127‑1 (MQ=255)
                                                                                                                                             ccGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATg                           >  3:283602/1‑127 (MQ=255)
                                                                                                                                              cGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCAGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTTGATTGCAACTGAg    >  6:270986/1‑149 (MQ=255)
                                                                                                                                                aTTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTGCAACTGAGGc  <  4:275887/149‑1 (MQ=255)
                                                                                                                                                aTTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTGCAACTGAGGc  <  4:79470/149‑1 (MQ=255)
                                                                                                                                                aTTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTGCAACTGAGGc  >  2:235470/1‑149 (MQ=255)
                                                                                                                                                 ttCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGaa                                                                                 <  4:29237/69‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCACCGCTGCTCACAAGGCTGACTTCAGTCCTTTTCCGTCATCTCAGCAAGATAATGGGTAGCAATTCCGTTGCTGAATGGCGCTGGTGTTCAAGTCCCCGTATGGGTGAACGGGGACAGCAGGCTGTGTACAACCCTGGGCCGATTCGGCCTCATGGCCGGTTTCGGGTTGAGGCTGCCCACTGGCATTGGACGAACCTAGGCACTGTACGAAGTCCCGACAAACCCAATCGTCTCGCCTCCAGCGCTGAGATTTTGTAGACTCATGTTGCTGCGTGGATTGCAACTGAGGC  >  NC_002947/2168293‑2168585

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: