Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,251,601 T→C 5.1% M83V (ATG→GTG)  PP_1094 ← electron transport complex protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,251,6010TC 5.1% 634.9 / 13.3 196M83V (ATG→GTG) PP_1094electron transport complex protein
Reads supporting (aligned to +/- strand):  ref base T (89/97);  new base C (2/8);  total (91/105)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 197 total aligned reads displayed.

CGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGCGA  >  NC_002947/1251453‑1251750
                                                                                                                                                     |                                                                                                                                                     
cGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATa                                                                                                                                                      <  2:1780533/150‑1 (MQ=255)
 gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAg                                                                                                                                                     <  2:4888298/150‑1 (MQ=255)
 gATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAg                                                                                                                                                     <  1:2533175/150‑1 (MQ=255)
     gCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGc                                                                                                                                                 >  1:4983407/1‑150 (MQ=255)
     gCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGc                                                                                                                                                 <  2:5087405/150‑1 (MQ=255)
         aGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACgg                                                                                                                                             >  2:624862/1‑150 (MQ=255)
          gCCTGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGt                                                                                                                                            >  1:1063776/1‑150 (MQ=255)
             cGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCGGTGTGAGCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAg                                                                                                                                         >  2:3439919/1‑150 (MQ=255)
              ggggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAgc                                                                                                                                        <  1:4589180/150‑1 (MQ=255)
                ggCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAgccg                                                                                                                                      <  2:2796198/150‑1 (MQ=255)
                 gCTTTCCTGTTGTTCGACATCCGGGCTGCCGCTCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCGGGGCCTGTGTGATCAGGAGAATGGCGCTGGGCTTACCGTGCCACTATGGCAGGTAAgccaa                                                                                                                                     <  2:3505219/150‑3 (MQ=255)
                 gCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCGTGCCACTATGGCAGGTAAgccaa                                                                                                                                     <  1:2903167/150‑3 (MQ=255)
                    ttCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAg                                                                                                                                  <  1:98110/150‑1 (MQ=255)
                       cTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGAt                                                                                                                               <  1:4274918/150‑1 (MQ=255)
                         gttgttCGACAACCTGGCTGCCTTTCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGGTTTCCTGCCTTAGTGGCACGGTAAGCCGCCAGGATAc                                                                                                                             <  1:3946429/150‑1 (MQ=255)
                          ttgttCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCCCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACg                                                                                                                            >  1:1109291/1‑150 (MQ=255)
                                gACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTAtt                                                                                                                      >  2:591607/1‑150 (MQ=255)
                                 aCAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTg                                                                                                                     >  1:2736635/1‑150 (MQ=255)
                                 aCAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTg                                                                                                                     >  1:3150520/1‑150 (MQ=255)
                                  cAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTgg                                                                                                                    >  1:1564365/1‑150 (MQ=255)
                                   aaCCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTggg                                                                                                                   >  1:3689356/1‑150 (MQ=255)
                                      cTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCct                                                                                                                >  2:3918361/1‑150 (MQ=255)
                                            gCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCctct                                                                                                              >  2:4745594/1‑146 (MQ=255)
                                            gCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAACGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgc                                                                                                          <  2:800475/150‑1 (MQ=255)
                                               gATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGGCTCTGTgcggc                                                                                                       >  2:2076091/1‑150 (MQ=255)
                                               gATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggc                                                                                                       >  1:1237493/1‑150 (MQ=255)
                                               gATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggc                                                                                                       >  2:484716/1‑150 (MQ=255)
                                                aTCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggcg                                                                                                      <  1:1544868/150‑1 (MQ=255)
                                                aTCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTgcggcg                                                                                                      >  2:5018819/1‑150 (MQ=255)
                                                 tCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGcggcgg                                                                                                     <  1:2374637/150‑1 (MQ=255)
                                                    gccggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGcc                                                                                                  >  1:4882405/1‑150 (MQ=255)
                                                      cggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCtt                                                                                                <  1:349352/150‑1 (MQ=255)
                                                      cggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCtt                                                                                                <  1:3244077/150‑1 (MQ=255)
                                                       ggccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTg                                                                                               <  1:2625568/150‑1 (MQ=255)
                                                        gccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTgg                                                                                              <  2:4584189/150‑1 (MQ=255)
                                                        gccgTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTgg                                                                                              <  2:2184623/150‑1 (MQ=255)
                                                              cgtGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCgg                                                                                        <  2:3281533/150‑1 (MQ=255)
                                                                   tgATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCttt                                                                                   >  1:1700215/1‑150 (MQ=255)
                                                                    gATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTc                                                                                  >  1:4204140/1‑150 (MQ=255)
                                                                       ggCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAg                                                                               >  2:3266882/1‑150 (MQ=255)
                                                                       ggCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAg                                                                               <  2:1935323/150‑1 (MQ=255)
                                                                       ggCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAg                                                                               <  1:3543653/150‑1 (MQ=255)
                                                                        gCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAgg                                                                              >  1:1448317/1‑150 (MQ=255)
                                                                        gCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAgg                                                                              <  1:1201020/150‑1 (MQ=255)
                                                                         cGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAggg                                                                             <  2:2533264/150‑1 (MQ=255)
                                                                              cTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTgg                                                                        <  1:3396983/150‑1 (MQ=255)
                                                                              cTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTgg                                                                        <  2:3236812/150‑1 (MQ=255)
                                                                               tCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGt                                                                       <  1:4649997/150‑1 (MQ=255)
                                                                                cAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTTTTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGTGCTTTCCAGGGGCTGGtt                                                                      <  2:4099935/150‑1 (MQ=255)
                                                                                cAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGtt                                                                      <  1:1992417/150‑1 (MQ=255)
                                                                                cAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGtt                                                                      <  2:2969836/150‑1 (MQ=255)
                                                                                 aGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATATGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTg                                                                     >  2:1674589/1‑150 (MQ=255)
                                                                                  gCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTCTTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGt                                                                    >  2:4206520/1‑150 (MQ=255)
                                                                                  gCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGt                                                                    >  2:3700638/1‑150 (MQ=255)
                                                                                  gCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGt                                                                    >  2:4206263/1‑150 (MQ=255)
                                                                                    aCGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                  <  1:1810081/150‑1 (MQ=255)
                                                                                    aCGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                  <  1:1809776/150‑1 (MQ=255)
                                                                                    aCGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCa                                                                  <  2:3712986/150‑1 (MQ=255)
                                                                                           gCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAc                                                           >  1:5306434/1‑150 (MQ=255)
                                                                                           gCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAc                                                           >  1:4088295/1‑150 (MQ=255)
                                                                                            cccGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTAcc                                                          >  1:4144354/1‑150 (MQ=255)
                                                                                             ccGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCg                                                         >  1:220043/1‑150 (MQ=255)
                                                                                              cGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTAGTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGGTGTCATAGCTACCGt                                                        >  1:4012714/1‑150 (MQ=255)
                                                                                               ggCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATGGcgcggcta                                                       >  1:1324723/1‑143 (MQ=255)
                                                                                               ggCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATGGCTACCGTc                                                       >  1:2352134/1‑150 (MQ=255)
                                                                                                   tACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgc                                                   >  1:5305217/1‑150 (MQ=255)
                                                                                                    aCCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgca                                                  <  2:2674225/150‑1 (MQ=255)
                                                                                                      cccTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgcagc                                                >  2:885546/1‑150 (MQ=255)
                                                                                                        cTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCAc                                              >  2:3301117/1‑150 (MQ=255)
                                                                                                         tGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACg                                             <  1:885546/150‑1 (MQ=255)
                                                                                                            gCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGcc                                          <  1:73205/150‑1 (MQ=255)
                                                                                                               tgtgtgATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGc                                       >  1:1279355/1‑150 (MQ=255)
                                                                                                               tgtgtgATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGc                                       <  1:4521169/150‑1 (MQ=255)
                                                                                                                    gATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTg                                  >  2:3006856/1‑150 (MQ=255)
                                                                                                                         ggATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTggcggc                             <  1:1431175/150‑1 (MQ=255)
                                                                                                                          gATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCt                            <  2:2756503/150‑1 (MQ=255)
                                                                                                                          gATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCt                            >  1:2720112/1‑150 (MQ=255)
                                                                                                                           aTAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCtt                           <  2:2574727/150‑1 (MQ=255)
                                                                                                                            tAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCttt                                                                                   >  1:2959367/1‑93 (MQ=255)
                                                                                                                            tAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCttt                          <  1:3801892/150‑1 (MQ=255)
                                                                                                                              aTGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGc                        <  2:1448317/150‑1 (MQ=255)
                                                                                                                              aTGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGc                        <  1:3424872/150‑1 (MQ=255)
                                                                                                                                ggCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGgc                                                   >  1:648653/1‑121 (MQ=255)
                                                                                                                                 gcgcTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAgc                     >  2:1723738/1‑150 (MQ=255)
                                                                                                                                  cgcTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGGCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAgct                    <  1:2076091/150‑1 (MQ=255)
                                                                                                                                    ggcggCTTACCGTGCCACTATGGCAGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCggg                                                                                       >  2:1397568/4‑81 (MQ=21)
                                                                                                                                        gCTTACCGTGCCACTATGGCAGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGtt                                                                      <  2:3626159/94‑1 (MQ=37)
                                                                                                                                        gCTTACCGTGCCACTATGGCAGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGtt                                                                      >  1:3626159/1‑94 (MQ=37)
                                                                                                                                         cTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTc             <  2:724037/150‑1 (MQ=255)
                                                                                                                                          ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa            >  2:661399/1‑150 (MQ=255)
                                                                                                                                          ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa            <  2:2941412/150‑1 (MQ=255)
                                                                                                                                          ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa            <  2:2482024/150‑1 (MQ=255)
                                                                                                                                          ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa            <  2:2977961/150‑1 (MQ=255)
                                                                                                                                          ttACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCa            <  2:2736635/150‑1 (MQ=255)
                                                                                                                                           gaCCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAg           <  2:2412522/149‑1 (MQ=255)
                                                                                                                                              cTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAgcgc        >  2:5367579/1‑150 (MQ=255)
                                                                                                                                                tGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGcc       >  2:815342/1‑150 (MQ=255)
                                                                                                                                                  ccATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCgg     <  1:5143924/150‑1 (MQ=255)
                                                                                                                                                    aTAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGCg   <  1:3071433/150‑1 (MQ=255)
                                                                                                                                                    aTAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGCg   <  1:3734895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATCGCCCAGCCCGGGGCTTTCCTGTTGTTCGACAACCTGGCTGCCGATCAGCCGGCCGTCCGTGGTGATGGCGACGCTCAGCACGATCGGCCCGGCATACCCCTGGGCCTGTGTGATCAGGATAATGGCGCTTGGCTTACCTGCCATAGTGGCACGGTAAGCCGCCAGGATACGGCTATTGGGCCTCTGTGCGGCGGGCCTTGGCAGCGGGCTTTCCAGGGGCTGGTTGTCATAGCTACCGTCGGGCAGCACGGCCAGCCGCTGGCGGCTTTGCAGCTGCTTTTCAGCGCCGGCGA  >  NC_002947/1251453‑1251750

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