Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,241,611 G→C 9.6% noncoding (35/76 nt) PP_t25 → tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,241,6110GC 9.6% 228.7 / 13.8 73noncoding (35/76 nt)PP_t25tRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base C (4/3);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

AATTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTT  >  NC_002947/2241466‑2241760
                                                                                                                                                 |                                                                                                                                                     
aaTTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTgcatggcatg                                                                                                                                                   >  2:2284794/1‑150 (MQ=255)
 aTTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGc                                                                                                                                                  <  1:1559166/150‑1 (MQ=255)
      tttGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCaggag                                                                                                                                             <  2:767500/150‑1 (MQ=255)
            aGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGc                                                                                                                                       <  1:134027/150‑1 (MQ=255)
            aGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGc                                                                                                                                       <  2:1440242/150‑1 (MQ=255)
             acccTGCAATCCGAACTTAATATGCGCCCCGTCCACAGAGTTGAACGAAAACGAAACGCCAGGCACCCCCACTCCCAGTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCg                                                                                                                                      <  1:540662/149‑1 (MQ=37)
                cTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAGAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGtt                                                                                                                                   <  2:176427/150‑1 (MQ=255)
                cTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGtt                                                                                                                                   <  1:2284794/150‑1 (MQ=255)
                  gCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGGGGTTCg                                                                                                                                 >  1:608486/1‑150 (MQ=255)
                      tCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATcc                                                                                                                             >  2:568321/1‑150 (MQ=255)
                          aaCTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCt                                                                                                                         <  2:2016451/150‑1 (MQ=255)
                                      gcCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaa                                                                                                             >  1:2252617/1‑150 (MQ=255)
                                      gcCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCaa                                                                                                             >  1:863066/1‑150 (MQ=255)
                                        cccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAAt                                                                                                           >  2:2344724/1‑150 (MQ=255)
                                         ccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAAtt                                                                                                          <  1:108783/150‑1 (MQ=255)
                                         ccGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGACAGAGATTCCGCTGATCGCGACATTGTGAGCCTATATCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAAtt                                                                                                          >  1:1499118/1‑150 (MQ=21)
                                            tCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCg                                                                                                       >  1:335366/1‑150 (MQ=255)
                                               acaGCGATGAACGAAAACTAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTTTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCtt                                                                                                    >  2:10514/1‑150 (MQ=255)
                                                        aacgaaaacgaaACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTaa                                                                                           <  2:866827/150‑1 (MQ=255)
                                                           gaaaacgaaaCGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                        >  1:1677137/1‑150 (MQ=255)
                                                              actgaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGa                                                                                     >  2:1030945/4‑150 (MQ=255)
                                                                 gaaacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATaa                                                                                  >  2:2101640/1‑150 (MQ=255)
                                                                 gaaacgACAGGCCCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATaa                                                                                  <  2:891328/150‑1 (MQ=255)
                                                                   aacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTa                                                                                                                       <  1:242335/111‑1 (MQ=255)
                                                                   aacgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTa                                                                                                                       >  2:242335/1‑111 (MQ=255)
                                                                     cgCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCaa                                                                              <  2:1517864/150‑1 (MQ=255)
                                                                      gCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGt                                                                                                                                    <  2:1683209/95‑1 (MQ=255)
                                                                      gCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGt                                                                                                                                    >  1:1683209/1‑95 (MQ=255)
                                                                         aGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGtt                                                                                                                                   >  1:1711505/1‑93 (MQ=255)
                                                                         aGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGtt                                                                                                                                   <  2:1711505/93‑1 (MQ=255)
                                                                         aGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGAcc                                                                          <  1:2101640/150‑1 (MQ=255)
                                                                           gCAGCCAGCCAGGGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATACCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                        <  2:1938216/150‑1 (MQ=255)
                                                                           gCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                        <  2:863066/150‑1 (MQ=255)
                                                                            cACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTTCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCggg                                                                       <  1:2344724/150‑1 (MQ=255)
                                                                             aCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                        >  2:789504/1‑148 (MQ=255)
                                                                             aCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCgg                                                                        <  1:789504/148‑1 (MQ=255)
                                                                             aCCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGt                                                                      <  2:335366/150‑1 (MQ=255)
                                                                               ccagccagCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                        >  1:1128496/1‑130 (MQ=255)
                                                                               ccagccagCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGc                                                                                        <  2:1128496/130‑1 (MQ=255)
                                                                                 agccagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAg                                                                  <  2:1499118/150‑1 (MQ=255)
                                                                                    cagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagcag                                                                                               <  2:906121/118‑1 (MQ=255)
                                                                                    cagcGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagcag                                                                                               >  1:906121/1‑118 (MQ=255)
                                                                                           tCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTc                                                        <  1:1030945/150‑1 (MQ=255)
                                                                                             aGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTc                                                                                                         >  2:1207741/1‑99 (MQ=255)
                                                                                             aGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTc                                                                                                         <  1:1207741/99‑1 (MQ=255)
                                                                                             aGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCCCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCtt                                                      >  2:74433/1‑150 (MQ=255)
                                                                                               ttCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCACGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGa                                                    >  2:2001080/1‑150 (MQ=255)
                                                                                               tgctaaGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagca                                                                                                <  1:740485/102‑1 (MQ=255)
                                                                                               tgctaaGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTagca                                                                                                >  2:740485/5‑106 (MQ=255)
                                                                                                          aTTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAgg                                         >  1:788185/1‑150 (MQ=255)
                                                                                                                         cTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCggg                                                          >  1:1901425/1‑118 (MQ=255)
                                                                                                                         cTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCggg                                                          <  2:1901425/118‑1 (MQ=255)
                                                                                                                               tGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTgg                    >  2:1673782/1‑150 (MQ=255)
                                                                                                                               tGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTgg                    >  1:1590462/1‑150 (MQ=255)
                                                                                                                                 ggAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGAcc                                                                          <  1:2036419/94‑1 (MQ=255)
                                                                                                                                 ggAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGAcc                                                                          >  2:2036419/1‑94 (MQ=255)
                                                                                                                                    gagCGCCTGCCTGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAg               <  1:1573461/150‑1 (MQ=255)
                                                                                                                                    gagCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAg               <  1:568321/150‑1 (MQ=255)
                                                                                                                                     agCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAg                       >  1:1094161/1‑141 (MQ=255)
                                                                                                                                     agCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAg                       <  2:1094161/141‑1 (MQ=255)
                                                                                                                                       ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCtt                                                      <  1:1201227/106‑1 (MQ=255)
                                                                                                                                       ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCtt                                                      >  2:1201227/3‑108 (MQ=255)
                                                                                                                                       ctccTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGac            >  2:627337/3‑150 (MQ=255)
                                                                                                                                       cgcCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTcccc                                <  1:799780/130‑1 (MQ=255)
                                                                                                                                       cgcCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTcccc                                >  2:799780/1‑130 (MQ=255)
                                                                                                                                        tccTGCATGCCATGCAGGCGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTaagaag                                          <  1:1119981/118‑1 (MQ=255)
                                                                                                                                        tccTGCATGCCATGCAGGCGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTaagaag                                          >  2:1119981/2‑119 (MQ=255)
                                                                                                                                         ccTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGacac          >  2:1620771/1‑150 (MQ=255)
                                                                                                                                          cTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAgg                                                                                         <  2:366842/70‑1 (MQ=255)
                                                                                                                                          cTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAgg                                                                                         >  1:366842/1‑70 (MQ=255)
                                                                                                                                          cTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACAcc         >  1:731961/1‑150 (MQ=255)
                                                                                                                                          cTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACAcc         <  1:1620771/150‑1 (MQ=255)
                                                                                                                                          cTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCa                                                                   <  1:732725/92‑1 (MQ=255)
                                                                                                                                          cTGCATGCCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCa                                                                   >  2:732725/1‑92 (MQ=255)
                                                                                                                                            gcatgccatgCAGGCGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGc       >  1:1577329/1‑150 (MQ=255)
                                                                                                                                                ggcatgCAGGAGGTCAGCGGTTCGCGCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCAGATAAGGAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCtt   >  2:2046290/1‑150 (MQ=255)
                                                                                                                                                 gcatgCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCttt  <  2:788185/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AATTAGTTTGACAGCCCTGCAATCCGAACTTAATATGCGCCCCGTCCACAGCGATGAACGAAAACGAAACGCCAGGCACCCAGCCAGCGATTCAGTTCAAAGCGACATTGTGGGGCTATAGCTCAGCTGGGAGAGCGCCTGCATGGCATGCAGGAGGTCAGCGGTTCGATCCCGCTTAGCTCCACCAAATTCCGCTTAGCAGCTAAGGCTGATAAGCAAGACCGGGTCCAGCAACCGGGTCTTGAAGGTAAGAAGGTTCGTCCCCTTCGTCTAGTGGCCTAGGACACCGCCCTTT  >  NC_002947/2241466‑2241760

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