Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,023,926 A→T 6.8% Y407* (TAT→TAA) ‡ opdP ← glycine‑glutamate dipeptide porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,023,9260AT 6.8% 217.9 / 5.2 74Y407* (TAT→TAA) ‡opdPglycine‑glutamate dipeptide porin
Reads supporting (aligned to +/- strand):  ref base A (30/39);  new base T (2/3);  total (32/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CGGAATGGTGGTGACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATAGGAAAG  >  NC_002947/1023782‑1024069
                                                                                                                                                |                                                                                                                                               
cGGAACGGGGGTGACCACGCGGAACTCACCCAGGCTGCCAGCGCCCTGGGCCTTGCTGGCGCGATGGGCGGTGTAGGTGGCGCGGACGCCGGTGTCCTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGt                                                                                                                                             <  8:227305/149‑1 (MQ=255)
 ggAATGGTGGTGACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTg                                                                                                                                            <  5:171429/149‑1 (MQ=255)
 cgactggtgtTCACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGAGGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTg                                                                                                                                            <  8:195249/145‑1 (MQ=255)
  gAATGGTGGTGACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTgg                                                                                                                                           <  6:66353/149‑1 (MQ=255)
      ggtggtGACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGttt                                                                                                                                       >  4:96907/1‑149 (MQ=255)
       gtggtgACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTc                                                                                                                                      <  2:207248/149‑1 (MQ=255)
       gtggtgACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTATGAGTggggtat                                                                                                                                      <  3:216213/149‑5 (MQ=255)
       gtggtgACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGGCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTc                                                                                                                                      <  7:298471/149‑1 (MQ=255)
          gtgACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGcc                                                                                                                                   <  6:202269/149‑1 (MQ=255)
               cACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATc                                                                                                                                >  7:277869/1‑147 (MQ=255)
                 cgcgGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTATGAGTgggg                                                                                                                                         <  1:61872/136‑2 (MQ=255)
                 cgcgGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTATGAGTgggg                                                                                                                                         >  2:61872/1‑135 (MQ=255)
                  gcgGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGc                                                                                                                           >  2:40310/1‑149 (MQ=255)
                    ggAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCAATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGTAACGCGTAGCTGGTGCCGATACCCCACTCATAGt                                                                                                                                             >  5:162836/1‑129 (MQ=255)
                    ggAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCAATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGTAACGCGTAGCTGGTGCCGATACCCCACTCATAGt                                                                                                                                             <  6:162836/129‑1 (MQ=255)
                         tCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACTCGTAGCTGGTGCCTATACCCCACTATGAGTggggta                                                                                                                                       <  5:331532/130‑4 (MQ=255)
                         tCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCAc                                                                                                                    >  8:23397/1‑149 (MQ=255)
                         tCATACAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAAAGCGTAGCTGGTGCCGATACCCCACTATGAGTggggta                                                                                                                                       >  6:331532/1‑127 (MQ=255)
                           aTCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCAc                                                                                                                    <  4:19378/147‑1 (MQ=255)
                              cAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCAcgtcg                                                                                                                >  2:298201/1‑148 (MQ=255)
                                ggCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAg                                                                                                             >  3:72799/1‑149 (MQ=255)
                                 gCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTggg                                                                                                                                          >  4:197057/1‑119 (MQ=255)
                                 gCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTggg                                                                                                                                          <  3:197057/119‑1 (MQ=255)
                                     ccATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGGTAGTGTCTTTCTGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAg                                                                                                         <  1:115579/148‑1 (MQ=255)
                                       aTCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGc                                                                                                      >  3:329701/1‑149 (MQ=255)
                                       aTCGCCCTGGGACTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGc                                                                                                      >  3:259638/1‑149 (MQ=255)
                                        tCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGATTCGCCAACCAGCCCATACACGTCGAAGGCAGTGcc                                                                                                     >  2:285903/1‑149 (MQ=255)
                                           cccTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCAcgtcgt                                                                                                               <  4:36225/136‑1 (MQ=255)
                                           cccTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCAcgtcgt                                                                                                               >  3:36225/1‑136 (MQ=255)
                                                ggCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGt                                                                                             <  4:272816/149‑1 (MQ=255)
                                                 gCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTg                                                                                            <  1:298201/149‑1 (MQ=255)
                                                  ccTTGCTTGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCGTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTgg                                                                                           <  4:324525/149‑1 (MQ=255)
                                                  ccTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTgg                                                                                           <  5:354818/149‑1 (MQ=255)
                                                       cTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTgg                                                                                                                                           >  1:158169/1‑96 (MQ=255)
                                                       cTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTgg                                                                                                                                           <  2:158169/96‑1 (MQ=255)
                                                           cgcgATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTc                                                                                   >  2:205677/1‑148 (MQ=255)
                                                              gATGGGTGGTGTAGGTGGCGCGGATGCTGGTGGCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATa                                                                               <  2:269775/149‑1 (MQ=255)
                                                                  ggtggtGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATAcccc                                                                           <  1:76454/149‑1 (MQ=255)
                                                                  ggtggtGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATAccc                                                                            >  3:268058/1‑148 (MQ=255)
                                                                      gtgtAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGcc                                                                       >  4:185768/1‑149 (MQ=255)
                                                                        gtAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATATCCCAGCcgc                                                                     <  4:128188/149‑1 (MQ=255)
                                                                            gTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCAt                                                                                                                                                <  8:297639/70‑1 (MQ=255)
                                                                            gTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCAt                                                                                                                                                >  7:297639/1‑70 (MQ=255)
                                                                               gcgcGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGtt                                                              <  1:285903/149‑1 (MQ=255)
                                                                                cgcgGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTg                                                             <  3:96907/149‑1 (MQ=255)
                                                                                cgcgGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTg                                                             <  7:23397/149‑1 (MQ=255)
                                                                                           gtgtCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGcc                                                                                      >  1:99701/1‑113 (MQ=255)
                                                                                           gtgtCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGcc                                                                                      <  2:99701/113‑1 (MQ=255)
                                                                                                tttCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACtt                                               >  4:326671/1‑147 (MQ=255)
                                                                                                tttCAGCGGCCCCCTCTGCACCGCGTAGCAGGTCCCCATACCCAACTCATAGTGGGTTTCGCCAACCAGGCCATTCACCTCGTAGGAAGAGCCGCGCTAGACGCTGCCGTCGATACCAAAGCCCCGGGCGATGTACAGGTTGAACtt                                               >  4:274275/1‑147 (MQ=255)
                                                                                                  tCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTca                                             <  3:185768/147‑1 (MQ=255)
                                                                                                    aGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCAATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTcaggc                                          >  8:185328/1‑148 (MQ=255)
                                                                                                      cGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTcaggccag                                       >  3:272990/1‑149 (MQ=255)
                                                                                                        gCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGAt                                                                                <  4:10811/106‑1 (MQ=255)
                                                                                                        gCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGAt                                                                                >  3:10811/1‑106 (MQ=255)
                                                                                                        gCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTcaggccaggc                                     <  4:268058/149‑1 (MQ=255)
                                                                                                        gCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTcaggccagg                                      >  4:108626/1‑148 (MQ=255)
                                                                                                           cGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGTCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACg                                  >  1:335536/1‑149 (MQ=255)
                                                                                                                  cAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAgg                           >  7:127441/1‑149 (MQ=255)
                                                                                                                   aaCGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAggg                          >  7:324502/1‑149 (MQ=255)
                                                                                                                        gTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTg                                                    <  1:83046/118‑1 (MQ=255)
                                                                                                                        gTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTg                                                    >  2:83046/1‑118 (MQ=255)
                                                                                                                            cTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGtt                 <  3:108626/149‑1 (MQ=255)
                                                                                                                               gTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCCCGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTc                <  3:274275/147‑1 (MQ=255)
                                                                                                                               gTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTAGAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCa               <  1:205677/148‑1 (MQ=255)
                                                                                                                               gTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCa               <  5:99331/148‑1 (MQ=255)
                                                                                                                                tGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGTCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGc             <  2:335536/149‑1 (MQ=255)
                                                                                                                                tGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGc             <  1:40310/149‑1 (MQ=255)
                                                                                                                                 gCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGca            >  5:68358/1‑149 (MQ=255)
                                                                                                                                     aTACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATa        <  8:277869/149‑1 (MQ=255)
                                                                                                                                     aTACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATa        <  3:326671/149‑1 (MQ=255)
                                                                                                                                        ccccACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATAGGa     <  8:324502/149‑1 (MQ=255)
                                                                                                                                         cccACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATAGGa     >  8:191366/1‑148 (MQ=255)
                                                                                                                                           cACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATAGGAAAg  <  8:127441/149‑1 (MQ=255)
                                                                                                                                                aTAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGcc                      >  6:168970/1‑124 (MQ=255)
                                                                                                                                                aTAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGcc                      <  5:168970/124‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
CGGAATGGTGGTGACCACGCGGAACTCATCCAGGCTGCCATCGCCCTGGGCCTTGCTGGCGCGATGGGTGGTGTAGGTGGCGCGGATGCTGGTGTCTTTCAGCGGCCCGCTCTGCAACGCGTAGCTGGTGCCGATACCCCACTCATAGTGGGTTTCGCCATCCAGGCCATTCACGTCGTAGGCAGTGCCGCGGTAGTGGGTGCCGTCGATACCCCAGCCGCGGGCGTTGTACAGGTTGAACTTCAGGCCAGGCACGCCATAGGGCGCCATGTTCAGCACATAGGAAAG  >  NC_002947/1023782‑1024069

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