Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,137,449 G→C 5.8% T34S (ACT→AGT)  PP_0998 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,137,4490GC 5.8% 332.9 / 9.6 104T34S (ACT→AGT) PP_0998hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (54/44);  new base C (3/3);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.72e-01

Only 100 of 105 total aligned reads displayed.

GTCCTTGGCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGCCGGCTTTGC  >  NC_002947/1137301‑1137598
                                                                                                                                                    |                                                                                                                                                     
gTCCTTGGCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGt                                                                                                                                                      <  2:617762/150‑1 (MQ=255)
      ggCCAGGAAGAACTTTTCTTCCACCGGCTTGCTTTCTACGCGGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTc                                                                                                                                                <  2:2000126/150‑1 (MQ=255)
       gCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCg                                                                                                                                               >  2:310758/1‑150 (MQ=255)
       gCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCACTGCGCtga                                                                                                                                               >  2:2344995/1‑148 (MQ=255)
       gCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCACTGCGCtga                                                                                                                                               >  2:2345474/1‑148 (MQ=255)
        ccAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGa                                                                                                                                              <  2:968781/150‑1 (MQ=255)
        ccAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGa                                                                                                                                              >  2:810151/1‑150 (MQ=255)
        ccAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCACTGCGCtgaa                                                                                                                                              <  1:2345474/150‑4 (MQ=255)
        ccAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCACTGCGCtgaa                                                                                                                                              <  1:2344995/150‑4 (MQ=255)
        ccAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCACTGCGCtgaa                                                                                                                                              <  1:1218206/150‑4 (MQ=255)
         cAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAt                                                                                                                                             >  2:1185153/1‑150 (MQ=255)
         cAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAt                                                                                                                                             >  2:2024849/1‑150 (MQ=255)
          aGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAtg                                                                                                                                            >  1:689295/1‑150 (MQ=255)
            gaagaaCTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAtgtg                                                                                                                                          >  1:990395/1‑150 (MQ=255)
            gaagaaCTTTTCTTCCACCGGCTTGCCTTCGATGCCGCCATACATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAtgtg                                                                                                                                          <  1:140505/150‑1 (MQ=255)
               gaaCTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTAc                                                                                                                                       <  2:1545048/150‑1 (MQ=255)
                aaCTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACa                                                                                                                                      <  1:1915288/150‑1 (MQ=255)
                   ttttCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGa                                                                                                                                   <  1:812477/150‑1 (MQ=255)
                   ttttCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGa                                                                                                                                   >  2:1903918/1‑150 (MQ=255)
                      tcttcCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCg                                                                                                                                <  1:1120176/150‑1 (MQ=255)
                        ttcCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGa                                                                                                                              <  1:1076199/150‑1 (MQ=255)
                        ttcCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGa                                                                                                                              >  2:1050369/1‑150 (MQ=255)
                         tcCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAg                                                                                                                             <  2:1623687/150‑1 (MQ=255)
                           cACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAtgt                                                                                                                                           <  2:1922180/134‑1 (MQ=255)
                           cACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGAtgt                                                                                                                                           >  1:1922180/1‑134 (MQ=255)
                           cACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTg                                                                                                                           >  1:905987/1‑150 (MQ=255)
                              cGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATa                                                                                                                        >  2:1153463/1‑150 (MQ=255)
                                  ttGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTgg                                                                                                                    >  1:1831604/1‑150 (MQ=255)
                                  ttGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTgg                                                                                                                    <  1:463182/150‑1 (MQ=255)
                                     ccTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGcc                                                                                                                 >  1:886949/1‑150 (MQ=255)
                                      cTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGt                                                                                                                                                      >  2:113472/1‑112 (MQ=255)
                                      cTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGt                                                                                                                                                      <  1:113472/112‑1 (MQ=255)
                                      cTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCa                                                                                                                >  2:440808/1‑150 (MQ=255)
                                       ttCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCAt                                                                                                               >  1:2164583/1‑150 (MQ=255)
                                       ttCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCAt                                                                                                               >  1:742868/1‑150 (MQ=255)
                                         cGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATgg                                                                                                             <  1:1185153/150‑1 (MQ=255)
                                             gctgcCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTAt                                                                                                         <  1:1903918/150‑1 (MQ=255)
                                                gcCATGCATGTGAAACTGAATCCCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTAtctc                                                                                                      >  1:1881952/1‑150 (MQ=255)
                                                gcCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGGCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTAtctc                                                                                                      >  1:1327530/1‑150 (MQ=255)
                                                    tgcatgTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTg                                                                                                                                                     <  2:1143981/99‑1 (MQ=255)
                                                    tgcatgTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTg                                                                                                                                                     >  1:1143981/1‑99 (MQ=255)
                                                       atgTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGa                                                                                               <  2:689295/150‑1 (MQ=255)
                                                        tgtgAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGaa                                                                                              <  2:2117078/150‑1 (MQ=255)
                                                             aaCTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATAATGGGCCATGGCTATCTCCTTTGGAAGGTCg                                                                                         >  1:2119597/1‑150 (MQ=255)
                                                                     aCCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGc                                                                                 <  1:440808/150‑1 (MQ=255)
                                                                     aCCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGc                                                                                 <  2:238257/150‑1 (MQ=255)
                                                                       ctctGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGa                                                                                                                              <  1:537997/103‑1 (MQ=255)
                                                                       ctctGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGa                                                                                                                              >  2:537997/1‑103 (MQ=255)
                                                                       ctctGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCgg                                                                               >  1:497658/1‑150 (MQ=255)
                                                                        tctGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATAATGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCggg                                                                              <  2:2119597/150‑1 (MQ=255)
                                                                          tGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGcc                                                                            >  1:323870/1‑150 (MQ=255)
                                                                            ccttccttGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCcgc                                                                          >  2:627379/1‑150 (MQ=255)
                                                                            ccgtccttGGGGTCCGCCACGGCCTCGTACGCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCcgc                                                                          <  1:1010129/147‑1 (MQ=255)
                                                                             cttccttGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCcgcg                                                                         >  1:1777924/1‑150 (MQ=255)
                                                                                 cttGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAg                                                                     >  2:2385118/1‑150 (MQ=255)
                                                                                  ttGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCGTAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTaa                                                                      >  2:1393735/1‑148 (MQ=255)
                                                                                  ttGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTaa                                                                      <  1:1393735/148‑1 (MQ=255)
                                                                                   tGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGcc                                                                   <  1:1050369/150‑1 (MQ=255)
                                                                                       gTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTc                                                                                          >  1:2292819/1‑123 (MQ=255)
                                                                                         ccGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGa                                                             >  1:1437615/1‑150 (MQ=255)
                                                                                            ccACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTGCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCTGCCGACTg                                                          >  2:1739622/1‑150 (MQ=255)
                                                                                              aCGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCa                                                        >  1:845061/1‑150 (MQ=255)
                                                                                                   cTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCaa                                                   <  2:886949/150‑1 (MQ=255)
                                                                                                   cTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCaa                                                   >  2:877002/1‑150 (MQ=255)
                                                                                                    tCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTggg                                                                                                                   <  1:1232153/85‑1 (MQ=255)
                                                                                                      gTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACgg                                                >  1:1885910/1‑150 (MQ=255)
                                                                                                      gTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACgg                                                <  1:2024849/150‑1 (MQ=255)
                                                                                                      gTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACgg                                                <  1:810151/150‑1 (MQ=255)
                                                                                                       tACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGAATTCAACgga                                               >  1:155326/1‑149 (MQ=255)
                                                                                                         cTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCg                                                                                >  2:4793/1‑115 (MQ=255)
                                                                                                         cTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCg                                                                                <  1:4793/115‑1 (MQ=255)
                                                                                                         cTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAg                                             <  1:877002/150‑1 (MQ=255)
                                                                                                          tCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGc                                            >  2:167346/1‑150 (MQ=255)
                                                                                                            gATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCAt                                                                                                               >  1:2283383/1‑81 (MQ=255)
                                                                                                            gATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCAt                                                                                                               <  2:2283383/81‑1 (MQ=255)
                                                                                                              tGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGaa                                                                                              <  2:974687/96‑1 (MQ=255)
                                                                                                              tGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGaa                                                                                              >  1:974687/1‑96 (MQ=255)
                                                                                                                ttCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCtt                                                                                                   >  1:106283/1‑89 (MQ=255)
                                                                                                                ttCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCtt                                                                                                   <  2:106283/89‑1 (MQ=255)
                                                                                                                   aGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGACCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCg                                   >  1:1230489/1‑150 (MQ=255)
                                                                                                                        aTAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGcact                              >  2:21254/1‑147 (MQ=255)
                                                                                                                         tAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTtct                             >  2:609904/1‑150 (MQ=255)
                                                                                                                            cTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCgg                          <  2:485072/150‑1 (MQ=255)
                                                                                                                             tgtgGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGt                         >  1:1508778/1‑150 (MQ=255)
                                                                                                                             tgtgGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGt                         <  1:2385118/150‑1 (MQ=255)
                                                                                                                              gtgGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTg                        <  2:1437615/150‑1 (MQ=255)
                                                                                                                              gtgGTCGGCGTTGTAAAGCGACGTGCGCTCGATGTATACATGACTGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTg                        <  2:1831604/150‑1 (MQ=255)
                                                                                                                               tgGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTgg                       >  1:893550/1‑150 (MQ=255)
                                                                                                                                    tgCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACGGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAAcc                  >  2:1650079/2‑150 (MQ=255)
                                                                                                                                    ggCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTggg                                                                                                                   <  2:319324/53‑1 (MQ=255)
                                                                                                                                    ggCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTggg                                                                                                                   >  1:319324/1‑53 (MQ=255)
                                                                                                                                          atcgaGCGCACTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTCCATTCAACGGTACCTGATCCCCGGCTTCCCCGTGGAAACCCGacag            >  2:1069826/5‑147 (MQ=255)
                                                                                                                                           ttCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGt                                               <  2:2125610/114‑1 (MQ=255)
                                                                                                                                           ttCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGt                                               >  1:2125610/1‑114 (MQ=255)
                                                                                                                                           ttCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGc           >  2:1583852/1‑150 (MQ=255)
                                                                                                                                                  cAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGCCGGCttt    <  2:1327530/150‑1 (MQ=255)
                                                                                                                                                   aGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGa                                                             <  2:2318557/92‑1 (MQ=255)
                                                                                                                                                   aGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGa                                                             >  1:2318557/1‑92 (MQ=255)
                                                                                                                                                   aGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGCCGGCtgtg   >  2:2184786/1‑150 (MQ=255)
                                                                                                                                                    gTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGCCGGCTTTGc  <  1:1583852/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCCTTGGCCAGGAAGAACTTTTCTTCCACCGGCTTGCCTTCGATGCTGCCATGCATGTGAAACTGAATACCTCTGCCTTCCTTGGGGTCCGCCACGGCCTCGTACTCGATGTTCAGGTCATAGCTGTGGTCGGCGTTGTTCAGCGCAGTGCGCTCGATGTGTACATGACCGGAGTGATACTGGGCCATGGCTATCTCCTTTGGAAGGTCGAAAATGGCGGGCCGCGTAAGCCCGCCGACTGCATTCAACGGTAGCTGATCCCCGGCTTCTCGGTGGAAACCCGAGTGCCGGCTTTGC  >  NC_002947/1137301‑1137598

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: