Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 841,360 G→C 12.4% L564V (CTC→GTC)  PP_0725 ← TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947841,3600GC12.4% 179.7 / 13.0 64L564V (CTC→GTC) PP_0725TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base C (4/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCACTGTTGCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAGGCCTG  >  NC_002947/841215‑841504
                                                                                                                                                 |                                                                                                                                                
gcCACTGTTGCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGt                                                                                                                                               >  6:159177/1‑149 (MQ=255)
gcCACTGTTGCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGt                                                                                                                                               >  5:310307/1‑149 (MQ=255)
 ccACTGTTGCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTc                                                                                                                                              <  5:10879/149‑1 (MQ=255)
         gCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCa                                                                                                                                       <  1:89064/148‑1 (MQ=255)
          cGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAg                                                                                                                                      <  5:159177/148‑1 (MQ=255)
            ggTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGACCCTCTTgcgccagat                                                                                                                                   >  5:113663/1‑144 (MQ=255)
            ggTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGACCCTCTTgcgc                                                                                                                                        <  6:113663/144‑1 (MQ=255)
                gCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTg                                                                                                                                           <  4:110260/137‑1 (MQ=255)
                gCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTg                                                                                                                                           >  3:110260/1‑137 (MQ=255)
                gCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAggg                                                                                                                               <  7:335456/149‑1 (MQ=255)
                        ccTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGc                                                                                                                       >  5:291029/1‑149 (MQ=255)
                                 ccTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAGGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTc                                                                                                              >  5:47549/1‑149 (MQ=255)
                                          tACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCtt                                                                                                     >  8:43878/1‑149 (MQ=255)
                                          tACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCtt                                                                                                     >  8:38571/1‑149 (MQ=255)
                                            ccGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTg                                                                                                    <  1:2098/148‑1 (MQ=255)
                                            ccGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTg                                                                                                    <  1:310241/148‑1 (MQ=255)
                                               aaGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAggg                                                                                                                               >  3:301166/1‑118 (MQ=255)
                                               aaGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAggg                                                                                                                               <  4:301166/118‑1 (MQ=255)
                                                      cgcAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTg                                                                                         <  1:60574/149‑1 (MQ=255)
                                                         aGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCg                                                                                      >  3:44619/1‑149 (MQ=255)
                                                          gCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCgg                                                                                     <  6:310307/149‑1 (MQ=255)
                                                           ccGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGcc                                                                                                        <  7:249858/129‑1 (MQ=255)
                                                           ccGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGcc                                                                                                        >  8:249858/1‑129 (MQ=255)
                                                               ccAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTgg                                                                                                   <  8:104328/130‑1 (MQ=255)
                                                               ccAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTgg                                                                                                   >  7:104328/1‑130 (MQ=255)
                                                                cagcagGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCgg                                                                                     <  5:147633/143‑1 (MQ=255)
                                                                cagcagGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCgg                                                                                     >  6:147633/1‑143 (MQ=255)
                                                                 agcagGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTcgc                                                                              >  6:125693/1‑149 (MQ=255)
                                                                  gcagGGCGATCAAGGTAAGGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTcgcg                                                                             >  7:173330/1‑149 (MQ=255)
                                                                  gcagGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTcgcg                                                                             <  1:80418/149‑1 (MQ=255)
                                                                      ggCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGaa                                                                                                           <  1:330461/115‑1 (MQ=255)
                                                                      ggCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGaa                                                                                                           >  2:330461/1‑115 (MQ=255)
                                                                               ggTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAggg                                                                                                                               <  4:42744/86‑1 (MQ=255)
                                                                                      gTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCcaca                                                         <  6:291029/149‑1 (MQ=255)
                                                                                           gCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTgg                                                               <  1:298768/138‑1 (MQ=255)
                                                                                           gCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTgg                                                               >  2:298768/1‑138 (MQ=255)
                                                                                           gCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACt                                                    >  8:162510/1‑149 (MQ=255)
                                                                                            cAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACt                                                    <  8:173330/148‑1 (MQ=255)
                                                                                              aTGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCg                                                 >  2:130720/1‑149 (MQ=255)
                                                                                                   gcAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCggg                                                                                    >  5:193341/1‑109 (MQ=255)
                                                                                                   gcAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCggg                                                                                    <  6:193341/109‑1 (MQ=255)
                                                                                                   gcAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAg                                            <  6:47549/149‑1 (MQ=255)
                                                                                                   gcAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAg                                            >  3:69774/1‑149 (MQ=255)
                                                                                                            tGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGggcggc                                   <  4:65579/149‑1 (MQ=255)
                                                                                                            tGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGggcggc                                   <  5:125693/149‑1 (MQ=255)
                                                                                                             ggCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGggcggc                                   >  3:266754/1‑148 (MQ=255)
                                                                                                              gCGTTAAAGAGTCTCGGCTCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGGCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCACCTGggcgaag                                  <  6:261626/148‑4 (MQ=255)
                                                                                                                     agagTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTg                          <  4:69774/149‑1 (MQ=255)
                                                                                                                      gagTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTgg                         >  8:230642/1‑149 (MQ=255)
                                                                                                                           tcGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCggg                    >  6:198033/1‑149 (MQ=255)
                                                                                                                             ggATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCgggg                   <  4:266754/148‑1 (MQ=255)
                                                                                                                               aTCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAg                >  7:323091/1‑149 (MQ=255)
                                                                                                                                  cgttCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCGAGGTGGGCGGCCACTTCGTGGTCGGGGAAGtt              <  7:43878/145‑1 (MQ=255)
                                                                                                                                  cGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGtt              <  7:38571/148‑1 (MQ=255)
                                                                                                                                   ggTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTgg            <  4:294592/149‑1 (MQ=255)
                                                                                                                                       ctgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGc                                                              <  6:108834/93‑1 (MQ=255)
                                                                                                                                       ctgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGt                           >  1:10467/3‑130 (MQ=255)
                                                                                                                                       ctgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGt                           <  2:10467/128‑1 (MQ=255)
                                                                                                                                       ctgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGtt              <  4:226649/141‑1 (MQ=255)
                                                                                                                                       ctgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGtt              >  3:226649/3‑143 (MQ=255)
                                                                                                                                        tgcgcAAGACCCTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAg       >  2:142955/2‑149 (MQ=255)
                                                                                                                                         gcgcAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAgg      <  5:198033/149‑1 (MQ=255)
                                                                                                                                             aaGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAGGCCTg  >  4:108444/1‑149 (MQ=255)
                                                                                                                                             aaGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGAGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAGGCCTg  >  2:252603/1‑149 (MQ=255)
                                                                                                                                               gagGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAgg      <  4:206729/143‑1 (MQ=255)
                                                                                                                                               gagGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAgg      >  3:206729/1‑143 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GCCACTGTTGCGGGTCGCCGTGGCCCTGCAGGGCCTCACGGCTACCGAAGCCGGCGCAGCCGGCCAGCAGGGCGATCAAGGTAAAGGTGATGCAATGGCGCAGGAACATGGCGTTAAAGAGTCTCGGATCCGGTCAGGCGCAAGAGGGTCTTGCGCAGGATAGGGCTGTCGGCCTGGGCTTCGAAGCCCTTGGCCCATACCTGGCGGGCTTCGCGGCGCTTGCCGTTGGCCCACAGCACTTCGCCCAGGTGGGCGGCCACTTCGTGGTCGGGGAAGTTGGCAAAGGCCTG  >  NC_002947/841215‑841504

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