Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,030 C→G 5.1% V479L (GTG→CTG)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,0300CG 5.1% 257.9 / 3.2 79V479L (GTG→CTG) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base C (39/36);  new base G (2/2);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGA  >  NC_002947/460885‑461174
                                                                                                                                                 |                                                                                                                                                   
ggCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgc                                                                                                                                                  >  6:54171/1‑149 (MQ=255)
   aGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgctgc                                                                                                                                               >  1:111921/1‑149 (MQ=255)
   aGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgctgc                                                                                                                                               <  6:6858/149‑1 (MQ=255)
        aCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTgg                                                                                                                                            >  8:31289/1‑147 (MQ=255)
         ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGc                                                                                                                                           <  6:281613/147‑1 (MQ=255)
         ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                                                                                                                         <  3:63192/149‑1 (MQ=255)
         ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                                                                                                                         <  2:111921/149‑1 (MQ=255)
         ccccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAACCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGc                                                                                                                                           <  1:15267/147‑1 (MQ=255)
          cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAgg                                                                                                                                        >  4:359054/1‑149 (MQ=255)
          cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAgg                                                                                                                                        >  3:39102/1‑149 (MQ=255)
          cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCAgt                                                                                                                                        <  4:307971/149‑2 (MQ=255)
            cGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGcc                                                                                                                                      >  2:178619/1‑149 (MQ=255)
               aCCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgc                                                                                                                                   >  3:287013/1‑149 (MQ=255)
                 cgcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGgcg                                                                                                                                  >  8:131797/1‑148 (MQ=255)
                   cgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAg                                                                                                                               <  5:54171/149‑1 (MQ=255)
                     tCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTg                                                                                                                             >  7:215757/1‑149 (MQ=255)
                      cTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTg                                                                                                                             >  7:336388/1‑148 (MQ=255)
                             ttGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                                                                                                                      >  8:266863/1‑148 (MQ=255)
                               gCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa                                                                                                                    <  7:129579/148‑1 (MQ=255)
                               gCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAg                                                                                                                   >  6:213624/1‑149 (MQ=255)
                                cTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGcagta                                                                                                                  >  1:370518/1‑145 (MQ=255)
                                cTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAg                                                                                                                   <  7:131797/148‑1 (MQ=255)
                                   gAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTAtt                                                                                                               >  3:175063/1‑149 (MQ=255)
                                     gCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTcc                                                                                                             >  8:338059/1‑149 (MQ=255)
                                      cAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTcc‑‑c                                                                                                          <  4:175063/149‑1 (MQ=255)
                                       aGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACt                                                                                                                        >  1:226193/1‑136 (MQ=255)
                                       aGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACt                                                                                                                        <  2:226193/136‑1 (MQ=255)
                                         ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgctg                                                                                                                                                <  5:99842/110‑1 (MQ=255)
                                         ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgctg                                                                                                                                                >  6:99842/1‑110 (MQ=255)
                                         ggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCgg                                                                                                       >  5:66382/1‑149 (MQ=255)
                                           cccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCa                                                                                                    <  2:372486/149‑1 (MQ=255)
                                            ccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGc                                                                                                                            >  7:185189/1‑127 (MQ=255)
                                            ccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGc                                                                                                                            <  8:185189/127‑1 (MQ=255)
                                            ccTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCAc                                                                                                   >  1:197162/1‑149 (MQ=255)
                                                     gATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACt                                                                                          <  7:338059/149‑1 (MQ=255)
                                                      aTCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTg                                                                                         <  7:266863/149‑1 (MQ=255)
                                                        cAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCCGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCActg                                                                                                 >  4:272768/1‑139 (MQ=255)
                                                        cAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCCGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCActg                                                                                                 <  3:272768/139‑1 (MQ=255)
                                                         aGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCa                                                                                       <  8:215757/148‑1 (MQ=255)
                                                         aGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCa                                                                                       <  4:287013/148‑1 (MQ=255)
                                                         aGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCAGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCa                                                                                       <  6:90317/148‑1 (MQ=255)
                                                           ggCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGtt                                                                                    >  8:116536/1‑149 (MQ=255)
                                                              cGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCa                                                                                 <  2:197162/149‑1 (MQ=255)
                                                               gggggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCAc                                                                                <  7:31289/149‑1 (MQ=255)
                                                                ggggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACg                                                                               >  8:43552/1‑149 (MQ=255)
                                                                ggggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCAGTACAGCAGtgccggg                                                                               <  8:336388/149‑7 (MQ=255)
                                                                 gggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACg                                                                               >  2:107396/1‑148 (MQ=255)
                                                                   gCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAg                                                                            <  6:66382/149‑1 (MQ=255)
                                                                          cGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGt                                                                     >  2:324887/1‑149 (MQ=255)
                                                                           gTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGt                                                                     >  6:82845/1‑148 (MQ=255)
                                                                             gCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTcac                                                                  <  6:354278/149‑1 (MQ=255)
                                                                             gCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTca                                                                   >  6:100717/1‑148 (MQ=255)
                                                                              cAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTcac                                                                  >  6:128159/1‑148 (MQ=255)
                                                                              cAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTcac                                                                  <  4:39102/148‑1 (MQ=255)
                                                                                  gCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGt                                                             >  3:84952/1‑149 (MQ=255)
                                                                                  gCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGt                                                             >  7:380160/1‑149 (MQ=255)
                                                                                    agccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTc                                                           <  7:43552/149‑1 (MQ=255)
                                                                                    agccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTc                                                           <  5:34929/149‑1 (MQ=255)
                                                                                       cagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                        <  5:82845/149‑1 (MQ=255)
                                                                                                    gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAgg                                           >  3:346952/1‑149 (MQ=255)
                                                                                                       gcgtGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGc                                         <  3:304542/145‑1 (MQ=255)
                                                                                                       gCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGc                                         <  6:188815/148‑1 (MQ=255)
                                                                                                       gCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGc                                         <  5:128159/148‑1 (MQ=255)
                                                                                                            aGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTc                                   <  5:100717/149‑1 (MQ=255)
                                                                                                             gCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTcc                                  <  5:213624/149‑1 (MQ=255)
                                                                                                               aTCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAg                                <  1:178619/149‑1 (MQ=255)
                                                                                                                tCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGc                               <  8:36800/149‑1 (MQ=255)
                                                                                                                tCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCAGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGc                               >  3:192714/1‑149 (MQ=255)
                                                                                                                tCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCACTGCAGCTGGCAGGCCCGCGCACTCCAACAGCCAGTAATCC‑‑CCCGCACTGCTGAACAGCAGTAGCACGGAGGCGCGGTCACAACGTTCCTCCAGTGGCGCCACGGCCTGTTCCAg                                >  1:258654/1‑148 (MQ=255)
                                                                                                                  aaggccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGca                              >  6:53629/1‑148 (MQ=255)
                                                                                                                       caaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCTGGAGGACAATGGTTTACGGCAGTTGCCCGGAGGCGCGGACAAAACTTGCGGGCAGTGCCCCCAGGGCATGTTccgccaccc                           >  1:125880/1‑142 (MQ=255)
                                                                                                                       caaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                        >  3:27327/1‑149 (MQ=255)
                                                                                                                       caaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                        >  6:412677/1‑149 (MQ=255)
                                                                                                                       caaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGCGCCAGGGCCTGTTCCAGCACACGCTg                       >  6:94290/1‑149 (MQ=255)
                                                                                                                        aaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTg                       >  1:249113/1‑149 (MQ=255)
                                                                                                                           gcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTaa                    <  4:346952/149‑1 (MQ=255)
                                                                                                                            cTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCAGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAAc                   >  7:101691/1‑149 (MQ=255)
                                                                                                                                        cccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCAt        <  2:249113/148‑1 (MQ=255)
                                                                                                                                        cccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTAC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGgcg           >  4:225993/1‑145 (MQ=255)
                                                                                                                                        ccagcagcaGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCa         >  5:295987/4‑147 (MQ=255)
                                                                                                                                        ccagcagcaGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCa         >  5:377325/4‑147 (MQ=255)
                                                                                                                                         ccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCActgctg                                                                                              >  5:196436/1‑61 (MQ=255)
                                                                                                                                         ccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCActgctg                                                                                              <  6:196436/61‑1 (MQ=255)
                                                                                                                                             agcaGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACTGCTGGCGCATAGGTGa  <  6:377325/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCC‑‑CCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGA  >  NC_002947/460885‑461174

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: