Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,204,842 C→G 6.0% R993P (CGC→CCC)  PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,204,8420CG 6.0% 215.6 / 4.8 67R993P (CGC→CCC) PP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base C (38/25);  new base G (2/2);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GCTGCTTGAGCTCTACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTGG  >  NC_002947/4204693‑4204991
                                                                                                                                                      |                                                                                                                                                     
gctgctTGAGCTCTACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGAcgc                                                                                                                                                       >  2:1060197/1‑150 (MQ=255)
   gctTGAGCTCTACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCGTGACGCgag                                                                                                                                                    >  2:1019835/1‑150 (MQ=255)
      tGAGCTCTACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGagaag                                                                                                                                                 <  2:414898/150‑1 (MQ=255)
           tctACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGATATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCag                                                                                                                                            >  1:1143302/1‑150 (MQ=255)
           tctACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCag                                                                                                                                            >  2:1027771/1‑150 (MQ=255)
            ctactaTTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCGGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCaga                                                                                                                                           >  2:124459/1‑150 (MQ=255)
                      gCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGGGAGAAGACCAGAGCGGCGAAAg                                                                                                                                 >  1:250445/1‑150 (MQ=255)
                      gCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGGGAGAAGACCAGAGCGGCGAAAg                                                                                                                                 >  1:250433/1‑150 (MQ=255)
                      gCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAg                                                                                                                                 <  1:399094/150‑1 (MQ=255)
                       ccAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGc                                                                                                                                >  1:883642/1‑150 (MQ=255)
                         aaGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCtt                                                                                                                              >  1:1027666/1‑150 (MQ=255)
                           gCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTaa                                                                                                                            >  2:329346/1‑150 (MQ=255)
                             aGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACg                                                                                                                          >  2:1165214/1‑150 (MQ=255)
                               tGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGtt                                                                                                                        >  2:777013/1‑150 (MQ=255)
                                tacagcaacgtcgtaccacctTCGGTGAGCTCAGCCCGGCCTGGCCGGAGAGTG‑TCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTg                                                                                                                       <  2:658986‑M1/129‑1 (MQ=255)
                                 ccACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGa                                                                                                                      <  2:945705/150‑1 (MQ=255)
                                   aCCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGagaag                                                                                                                                                 >  1:983392/1‑121 (MQ=255)
                                   aCCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGagaag                                                                                                                                                 <  2:983392/121‑1 (MQ=255)
                                   aCCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAgc                                                                                                                    >  1:1253969/1‑150 (MQ=255)
                                        cTGCGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCg                                                                                                               <  1:1011517/150‑1 (MQ=255)
                                           aGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCGTGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGgtg                                                                                                            <  1:1019835/150‑1 (MQ=255)
                                              taccaccttcggtgagctcagcccggcctggccggagagtgtcaccgaAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAgc                                                                                                                    <  2:121785‑M2/92‑1 (MQ=255)
                                              taccaccttcggtgagctcagcccggcctggccggagagtgtcaccgaAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAgc                                                                                                                    >  1:121785‑M2/49‑140 (MQ=255)
                                               aGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGtgttgt                                                                                                        <  1:1060197/150‑1 (MQ=255)
                                               aGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGtgttgt                                                                                                        <  1:777013/150‑1 (MQ=255)
                                                   cATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAATGTGTTACTACGTTGCTGTAGCAACTCGCTCCCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCaa                                                                                                    <  1:286907/150‑1 (MQ=255)
                                                     TCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGtgtt                                                                                                          >  2:177734‑M1/1‑142 (MQ=255)
                                                     TCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGtgtt                                                                                                          <  1:177734‑M1/142‑1 (MQ=255)
                                                       GGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCa                                                                                                >  2:480530‑M1/1‑150 (MQ=255)
                                                        GTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGATCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCAt                                                                                               >  2:584501‑M1/1‑150 (MQ=255)
                                                         TGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATg                                                                                              >  2:1027666‑M1/1‑150 (MQ=255)
                                                              acTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGcaca                                                                                          <  2:1253969/150‑1 (MQ=255)
                                                                        ctggCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGa                                                                                                                      <  1:512759/110‑1 (MQ=255)
                                                                        ctggCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGa                                                                                                                      >  2:512759/3‑112 (MQ=255)
                                                                         aGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGGGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCa                                                                               <  2:250445/150‑1 (MQ=255)
                                                                         aGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGGGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCa                                                                               <  2:250433/150‑1 (MQ=255)
                                                                                tccgAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGTTACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTg                                                                        <  1:72905/147‑1 (MQ=255)
                                                                                 cTGAGCTCACCGAAGGTGGTAGGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTgg                                                                       <  2:525313/150‑1 (MQ=255)
                                                                                 cTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTgg                                                                       <  2:850243/150‑1 (MQ=255)
                                                                                  tGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACCGACATCGCTCGGCCGttcga                                                                      >  2:72905/1‑147 (MQ=255)
                                                                                  tGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTggg                                                                      >  1:800940/1‑150 (MQ=255)
                                                                                   gAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCa                                                                                                     >  2:878934/1‑118 (MQ=255)
                                                                                   gAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCa                                                                                                     <  1:878934/118‑1 (MQ=255)
                                                                                       tcaccgaAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAAt                                                                  <  1:1193557‑M2/142‑1 (MQ=255)
                                                                                       tcaccgaAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAAt                                                                  >  2:1193557‑M2/8‑149 (MQ=255)
                                                                                           cgaAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAAc                                                             >  2:580230‑M2/4‑150 (MQ=255)
                                                                                              aGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACggg                                                          <  1:1165214/150‑1 (MQ=255)
                                                                                                gtggtACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCa                                                        >  2:673062/1‑150 (MQ=255)
                                                                                                       gacgTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGcc                                                 >  2:1218797/1‑150 (MQ=255)
                                                                                                       gacgTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGcc                                                 >  1:1100711/1‑150 (MQ=255)
                                                                                                       gacgTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGcc                                                 >  1:1100887/1‑150 (MQ=255)
                                                                                                          gTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTa                                              <  2:1143302/150‑1 (MQ=255)
                                                                                                          gTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTa                                              >  2:1089555/1‑150 (MQ=255)
                                                                                                                     aaCTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAgcgc                                   >  2:562899/1‑150 (MQ=255)
                                                                                                                      aCTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGcc                                  >  2:872903/1‑150 (MQ=255)
                                                                                                                        tCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCaa                                >  1:1071033/1‑150 (MQ=255)
                                                                                                                        tCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCaa                                >  1:497573/1‑150 (MQ=255)
                                                                                                                          gATACCCTGCTTCGCATCGGCCATGACGCGAGAAGGTGAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATc                              >  1:622647/1‑150 (MQ=255)
                                                                                                                              cccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGt                          >  1:255102/1‑150 (MQ=255)
                                                                                                                               ccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCtt                            <  2:407318/147‑1 (MQ=255)
                                                                                                                               ccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCtt                            >  1:407318/1‑147 (MQ=255)
                                                                                                                               ccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTc                         <  2:1071033/150‑1 (MQ=255)
                                                                                                                               ccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTc                         >  2:591153/1‑150 (MQ=255)
                                                                                                                               ccTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTc                         >  1:374523/1‑150 (MQ=255)
                                                                                                                                 tGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCct                       <  2:374523/150‑1 (MQ=255)
                                                                                                                                          tCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATg                                                                 <  2:263924/99‑1 (MQ=255)
                                                                                                                                          tCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATg                                                                 >  1:263924/1‑99 (MQ=255)
                                                                                                                                               gaTGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAg         >  2:882467/2‑150 (MQ=255)
                                                                                                                                                     gcgAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTg   >  1:42684/1‑150 (MQ=255)
                                                                                                                                                     gcgAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTg   >  2:987317/1‑150 (MQ=255)
                                                                                                                                                      cgAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTgg  <  2:255102/150‑1 (MQ=255)
                                                                                                                                                      cgAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTgg  <  1:1272879/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTGCTTGAGCTCTACTATTTGGCCAAGCAGTGCCACCGCCTGAGCAAGGACATCGGTGA‑CACTCTCCGGCCAGGCCGGGCTGAGCTCACCGAAGGTGGTACGACGTTGCTGTAGCAACTCGATACCCTGCTTCGCATCGGCCATGACGCGAGAAGACCAGAGCGGCGAAAGCTTAACGTTGAGCGCTCGGTGTTGTTCAACCCATGCACAGCCATCGCTCAGCCGTTGGGAAATGGAACGGGCATCAGGCCGTAGCGCGAAGCGCCAATCTTGTCCTGCTGCAGCTGGCAGTGCCTGG  >  NC_002947/4204693‑4204991

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