Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,083,641 T→G 5.0% F99V (TTC→GTC)  PP_1862 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,083,6410TG 5.0% 252.4 / 3.5 80F99V (TTC→GTC) PP_1862membrane protein
Reads supporting (aligned to +/- strand):  ref base T (36/39);  new base G (2/2);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAGCTTTGCCGGCACCATGG  >  NC_002947/2083496‑2083788
                                                                                                                                                  |                                                                                                                                                   
gATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCa                                                                                                                                                  >  4:215000/1‑149 (MQ=255)
gATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCa                                                                                                                                                  >  3:168615/1‑149 (MQ=255)
gATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGGTTCGCCATCCGTTGCCGCTGTCGGCGTGCAGCTACTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGGGGGGCTTGA‑TGCGCATATct                                                                                                                                                  >  6:382698/1‑148 (MQ=255)
          ccGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTgg                                                                                                                                         <  7:79112/148‑1 (MQ=255)
            gAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGtt                                                                                                                                      >  2:289134/1‑149 (MQ=255)
            gAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATGAAATgcgcatcaag                                                                                                                                       >  8:85911/1‑142 (MQ=255)
            gAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATGAAATgcgcatcaag                                                                                                                                       <  7:85911/148‑7 (MQ=255)
             aaCGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGttt                                                                                                                                     >  6:2000/1‑149 (MQ=255)
                gTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGt                                                                                                                                  <  3:215000/149‑1 (MQ=255)
                 tgctgcTGGCGCTGTGGCGGTATGACCAGCGCATCCGCATCGCGTTTCGGCTTGCGTTGGCGCTGTCGGGGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTc                                                                                                                                 <  7:330868/149‑1 (MQ=255)
                    tgctgGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGc                                                                                                                              >  3:227676/1‑149 (MQ=255)
                      ctgGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCtt                                                                                                                            <  5:223774/149‑1 (MQ=255)
                       tgGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTc                                                                                                                           >  5:266233/1‑149 (MQ=255)
                         gcgcTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCt                                                                                                                         <  6:350809/149‑1 (MQ=255)
                           gcTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCCTTTCATGCGCCTGGGTTTTGTCGGCTTCCTgg                                                                                                                       >  6:260806/1‑149 (MQ=255)
                                   ggTATCACCAGCGCAACCGCATCGCGTTTGGCCATTCGTTGCCGGTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTGGATCTGGGCGGTGGGGCTTGA‑TGCGGCTTTCATGCGCCTGGGGTTTGTGGGCTTGCTGGTAGCGTTg                                                                                                               >  1:220044/1‑149 (MQ=255)
                                   ggTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGGTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTg                                                                                                               <  3:266281/149‑1 (MQ=255)
                                   ggTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTg                                                                                                               >  5:204605/1‑149 (MQ=255)
                                   ggTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTg                                                                                                               >  5:368080/1‑149 (MQ=255)
                                    gTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTgg                                                                                                              <  5:260806/149‑1 (MQ=255)
                                         aCCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGgtg                                                                                                         >  2:229274/1‑149 (MQ=255)
                                         aCCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATGAAATgcgca                                                                                                                                            >  7:138355/1‑113 (MQ=255)
                                         aCCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATGAAATgcgca                                                                                                                                            <  8:138355/114‑2 (MQ=255)
                                          ccAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGgtgt                                                                                                        <  2:243341/149‑1 (MQ=255)
                                              cgcAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCTGCTTCTTGTTCGCGTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGac                                                                                                    <  3:73117/149‑2 (MQ=255)
                                               gcAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAAc                                                                                                   >  7:271402/1‑149 (MQ=255)
                                               gcAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAAc                                                                                                   <  8:219286/149‑1 (MQ=255)
                                               gcAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGGCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAAc                                                                                                   <  1:303741/149‑1 (MQ=255)
                                                 aaCCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACgt                                                                                                 >  2:13483/1‑149 (MQ=255)
                                                    cGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACgtggt                                                                                              <  6:166858/149‑1 (MQ=255)
                                                     gCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCg                                                                                                                                <  7:141082/114‑1 (MQ=255)
                                                     gCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACgtggtg                                                                                             <  4:168615/149‑1 (MQ=255)
                                                          gcgTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGAt                                                                                        <  1:289134/149‑1 (MQ=255)
                                                            gTTTCGACATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGt                                                                                      <  7:130399/149‑1 (MQ=255)
                                                                cGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCGGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCa                                                                                   <  6:266233/148‑1 (MQ=255)
                                                                 gCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCAt                                                                                  <  1:229274/148‑1 (MQ=255)
                                                                 gCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCAt                                                                                  >  7:317398/1‑148 (MQ=255)
                                                                 gCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATg                                                                                 <  5:59614/149‑1 (MQ=255)
                                                                  ccATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGGCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGGATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGTTGTTCATg                                                                                 >  5:77282/1‑148 (MQ=255)
                                                                   cATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGaa                                                                               <  8:326976/149‑1 (MQ=255)
                                                                    aTGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAg                                                                              <  3:20728/149‑1 (MQ=255)
                                                                    aTGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTGGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAg                                                                              <  2:220044/149‑1 (MQ=255)
                                                                          tGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTcg                                                                        >  7:216062/1‑149 (MQ=255)
                                                                             cGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGGGGTGGTTAATCAAGTTGTcgcgg                                                                     >  5:325813/1‑147 (MQ=255)
                                                                             cGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTcgccg                                                                     >  5:85517/1‑149 (MQ=255)
                                                                                   cGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGt                                                               >  8:329025/1‑149 (MQ=255)
                                                                                     gCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGc                                                             >  8:150099/1‑149 (MQ=255)
                                                                                      cGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGcc                                                            <  6:204605/149‑1 (MQ=255)
                                                                                                cttGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCg                                                    <  3:207844/147‑1 (MQ=255)
                                                                                                  tGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCAt                                                >  7:152664/1‑149 (MQ=255)
                                                                                                   gTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTTGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCAt                                                <  3:100606/148‑1 (MQ=255)
                                                                                                   gTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCAt                                                <  7:210106/148‑1 (MQ=255)
                                                                                                    ttCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATg                                               <  7:150099/148‑1 (MQ=255)
                                                                                                    ttCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGc                                              <  6:279966/149‑1 (MQ=255)
                                                                                                    ttCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGc                                              <  8:271402/149‑1 (MQ=255)
                                                                                                    ttCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGc                                              >  6:201608/1‑149 (MQ=255)
                                                                                                    ttCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCAAGCTGCGTCATGc                                              >  7:193539/1‑149 (MQ=255)
                                                                                                      cGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCt                                            >  5:89128/1‑149 (MQ=255)
                                                                                                         ctggctggGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGc                                         >  6:297695/1‑149 (MQ=255)
                                                                                                          tggctggGGTCGATCTGGGCGGTGGGGCTTGACTGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGc                                         >  1:45641/1‑149 (MQ=255)
                                                                                                          tggctggGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGcc                                        >  6:178819/1‑149 (MQ=255)
                                                                                                            gctggGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCg                                       <  6:368080/148‑1 (MQ=255)
                                                                                                                 ggTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTg                                 >  2:135810/1‑149 (MQ=255)
                                                                                                                  gTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTgg                                <  6:85517/149‑1 (MQ=255)
                                                                                                                  gTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTgg                                >  5:109193/1‑149 (MQ=255)
                                                                                                                  gTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAAAGTGGTGAGGATGTTCATGAAGGTGTCGCCGCAGAGTGCCTAGCTGAGTCATGCCTGGCCGTGGCTgg                                >  5:177549/1‑149 (MQ=255)
                                                                                                                     gATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGgtgt                                                                                                        <  2:267526/74‑1 (MQ=255)
                                                                                                                      aTCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGgtgt                            >  1:300938/1‑149 (MQ=255)
                                                                                                                       tCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGtgtg                           >  7:140841/1‑149 (MQ=255)
                                                                                                                           ggCGGTGGGGCTTGA‑TGCGCATTTCTTGCGCCTGGGTTGTGGCGGCTTCCTGGTAGCGTGGGCGGTGTGGAACGTGGTGCGGAGGTGCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCt                         <  5:382698/147‑1 (MQ=255)
                                                                                                                           ggCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGa                                                          >  3:108202/1‑114 (MQ=255)
                                                                                                                           ggCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGa                                                          <  4:108202/114‑1 (MQ=255)
                                                                                                                           ggCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCt                         <  5:297695/147‑1 (MQ=255)
                                                                                                                           ggCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTgg                       >  2:286653/1‑149 (MQ=255)
                                                                                                                              ggTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCa                    <  8:317398/149‑1 (MQ=255)
                                                                                                                               gtggtgCTTGA‑TGCGCATTTCATGCGACGGGGTTTTGTCGGCTTCCTTGTAGCGTTGGCGGTGTGGAACGTGGTGCGTATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTTGGTGTGCTGGTCAg                   <  8:191944/149‑1 (MQ=255)
                                                                                                                               gTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAg                   <  7:329025/149‑1 (MQ=255)
                                                                                                                                         a‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAGCTTTGCCgg          <  6:109193/148‑1 (MQ=255)
                                                                                                                                            gcgcATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAGCTTTGCCGGCAc       >  5:339693/1‑149 (MQ=255)
                                                                                                                                              gcATTTCATGCGCCTCGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCAGGGTGTGCTGGGCAGCTTTGCCGGCGCCa     >  8:71242/1‑149 (MQ=255)
                                                                                                                                                 tttCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAGCTTTGCCGGCACCATgg  >  8:174564/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGA‑TGCGCATTTCATGCGCCTGGGTTTTGTCGGCTTCCTGGTAGCGTTGGCGGTGTGGAACGTGGTGCGGATGTTCATGAAGGTGTCGCCGCCGAGTGCCGAGCTGCGTCATGCCTGGCCGTGGCTGGGTGTGCTGGGCAGCTTTGCCGGCACCATGG  >  NC_002947/2083496‑2083788

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: ‑