Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 184,740 T→C 5.8% Q291Q (CAA→CAG)  PP_0161 ← putative transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947184,7400TC 5.8% 202.1 / 2.4 70Q291Q (CAA→CAG) PP_0161putative transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base C (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

AGTGTTGCTTACTGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCAGTC  >  NC_002947/184593‑184887
                                                                                                                                                   |                                                                                                                                                   
aGTGTTGCTTACTGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGtt                                                                                                                                                    <  3:209368/149‑1 (MQ=255)
      ccTTACTGTCCCACCCGCACCTGGCATCACACCCCGCCAGAACTCGCTACCGCCTTATCAGGCGATGCTCTCTACATGCCCCCCCCAGCCCGTACCCCGACGACGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCg                                                                                                                                              <  4:2459/148‑1 (MQ=255)
        ttACTGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCggg                                                                                                                                            >  5:21313/1‑149 (MQ=255)
          aCTGTCCTACCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGGGATGCTCTCTACATCCACCCACCAGCGTGTACGCCGACCAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCggg                                                                                                                                            <  1:176849/147‑1 (MQ=255)
           cTGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCt                                                                                                                                         >  7:292988/1‑149 (MQ=255)
            tGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCtt                                                                                                                                        >  8:210021/1‑149 (MQ=255)
            tGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCtt                                                                                                                                        >  6:300336/1‑149 (MQ=255)
              tCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTc                                                                                                                                       >  6:73188/1‑148 (MQ=255)
                cTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAg                                                                                                                                    <  5:300336/149‑1 (MQ=255)
                   gCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGc                                                                                                                                  <  5:95899/148‑1 (MQ=255)
                    ccGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCg                                                                                                                                 >  6:257640/1‑148 (MQ=255)
                        aCCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAaca                                                                                                                            >  5:349538/1‑149 (MQ=255)
                          cTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACAcc                                                                                                                          >  5:31713/1‑149 (MQ=255)
                             gATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACAccccc                                                                                                                       <  5:659/149‑1 (MQ=255)
                                  aaaaCCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAg                                                                                                                  >  7:81010/1‑149 (MQ=255)
                                        cGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAg                                                                                                            <  1:308588/149‑1 (MQ=255)
                                                         gCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATtgc                                                                                           >  1:250965/1‑149 (MQ=255)
                                                                 aGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGccc                                                                                   <  7:192232/149‑1 (MQ=255)
                                                                  ggTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGcc                                                                                    >  5:118400/1‑147 (MQ=255)
                                                                  ggTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGGCGGCAATCTCATTGCTgccgcc                                                                                    >  8:197790/1‑147 (MQ=255)
                                                                       tGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTaa                                                                             <  5:73188/149‑1 (MQ=255)
                                                                           ctctACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTACCCAt                                                                         >  8:333249/1‑149 (MQ=255)
                                                                              tACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGACCAGGTAACCATGGc                                                                      >  1:81592/1‑149 (MQ=255)
                                                                                 aTGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGa                                                                    >  7:177880/1‑148 (MQ=255)
                                                                                 aTGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGa                                                                    >  2:211104/1‑148 (MQ=255)
                                                                                       ccaccaGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAAcc                                                                           >  1:131699/1‑135 (MQ=255)
                                                                                       ccaccaGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAAcc                                                                           <  2:131699/135‑1 (MQ=255)
                                                                                        caccaGCGTGTACGCCGACGAAGTCGGACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCt                                                            >  7:132736/1‑149 (MQ=255)
                                                                                         accaGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTATCCATGGCGAGAGCGGCt                                                            >  5:256535/1‑148 (MQ=255)
                                                                                                   aCGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTgcgc                                                   >  2:186855/1‑147 (MQ=255)
                                                                                                    cGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGccc                                                                                   >  5:46057/1‑114 (MQ=255)
                                                                                                    cGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGccc                                                                                   <  6:46057/114‑1 (MQ=255)
                                                                                                     gCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAg                                                >  5:263730/1‑148 (MQ=255)
                                                                                                     gCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGAGCGAg                                                >  7:240159/1‑148 (MQ=255)
                                                                                                      ccGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAgg                                               >  5:160289/1‑148 (MQ=255)
                                                                                                        gacgaAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAg                                            >  1:115793/1‑149 (MQ=255)
                                                                                                          cgaAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGtt                                          <  8:81010/149‑1 (MQ=255)
                                                                                                          cgaAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGGAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGtt                                          >  8:201827/1‑149 (MQ=255)
                                                                                                           gaAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTg                                         <  7:210021/149‑1 (MQ=255)
                                                                                                             aGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGga                                       <  7:333249/149‑1 (MQ=255)
                                                                                                             aGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGga                                       <  6:21313/149‑1 (MQ=255)
                                                                                                               tCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCGCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGgaga                                     <  7:197790/149‑1 (MQ=255)
                                                                                                                 gTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGTGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCGAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGagag                                    <  2:250965/148‑1 (MQ=255)
                                                                                                                 gTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGagag                                    <  2:22814/148‑1 (MQ=255)
                                                                                                                 gTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAgg                                   <  1:186855/149‑1 (MQ=255)
                                                                                                                    cGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGgcg                                 <  7:201827/148‑1 (MQ=255)
                                                                                                                     gggCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCa                               >  8:184270/1‑149 (MQ=255)
                                                                                                                     gggCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCa                               >  2:337290/1‑149 (MQ=255)
                                                                                                                             tCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAg                                                <  6:282052/124‑1 (MQ=255)
                                                                                                                             tCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAg                                                >  5:282052/1‑124 (MQ=255)
                                                                                                                             tCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGa                       <  6:263730/149‑1 (MQ=255)
                                                                                                                             tCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTCCGCGAGGAAGGTGGAGCGGCGCAGGTCCt                         >  6:212490/1‑147 (MQ=255)
                                                                                                                              cTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAg                      <  8:292988/149‑1 (MQ=255)
                                                                                                                               tGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGt                     >  2:276050/1‑149 (MQ=255)
                                                                                                                                gCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGt                     >  5:20852/1‑148 (MQ=255)
                                                                                                                                 cGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGa                   >  2:329018/1‑149 (MQ=255)
                                                                                                                                 cGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGa                   >  6:307598/1‑149 (MQ=255)
                                                                                                                                 cGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGa                   >  3:119885/1‑149 (MQ=255)
                                                                                                                                        agcCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCaa            >  6:250161/1‑149 (MQ=255)
                                                                                                                                         gcCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAAc           <  6:118400/149‑1 (MQ=255)
                                                                                                                                          ccGGAGCAACTGCTTCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGgc                                  <  2:26242/125‑1 (MQ=255)
                                                                                                                                          ccGGAGCAACTGCTTCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGgc                                  >  1:26242/1‑125 (MQ=255)
                                                                                                                                          ccGGAGCAACTGCTTCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAAc           <  8:190999/148‑1 (MQ=255)
                                                                                                                                          ccGGAGCAACTGCTTCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAAc           >  7:190999/1‑148 (MQ=255)
                                                                                                                                          ccGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAAcc          <  5:250161/149‑1 (MQ=255)
                                                                                                                                          ccGAAGCAGGTGCTCCGGGTCTTCGAGGCTGAACACCCCCGACATGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAAcc          <  6:256535/149‑1 (MQ=255)
                                                                                                                                             aaGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGGCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCtt       <  5:257640/149‑1 (MQ=255)
                                                                                                                                               gCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCa     <  7:184270/149‑1 (MQ=255)
                                                                                                                                               gCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCa     <  6:349538/149‑1 (MQ=255)
                                                                                                                                                 aGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCAGt   >  6:40033/1‑149 (MQ=255)
                                                                                                                                                  gTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCAGTc  <  2:115793/149‑1 (MQ=255)
                                                                                                                                                   ttGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCAGTc  >  7:45967/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTGTTGCTTACTGTCCTAGCCGGACCTGGATTCAAAACCCGCCAGAAATCGATAACGCCTTATCAGGTGATGCTCTCTACATGCACCCACCAGCGTGTACGCCGACGAAGTCGTACGGGCAACGTCTGCGGTAGCAGCCGAAGCAGTTGCTCCGGGTCTTCGAGGCGGAACACCCCCGACAGGCGCAGGTCGGCAATCTCATTGCTGCAGCCCAGGTAACCATGGCGATAGCGGCTGACCTGCGCGAGGAAGTTGGAGAGGCGCATGTCCTGAGTGACGATCAACCCTTCAGTC  >  NC_002947/184593‑184887

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