Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,654,136 G→C 10.9% P878R (CCG→CGG)  PP_1449 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,654,1360GC10.9% 184.7 / 11.2 64P878R (CCG→CGG) PP_1449hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base C (2/5);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-02

GACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGT  >  NC_002947/1653991‑1654284
                                                                                                                                                 |                                                                                                                                                     
gACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCgg                                                                                                                                                      <  5:327764/147‑1 (MQ=255)
 aCCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCgggg                                                                                                                                                    <  6:364063/148‑1 (MQ=255)
     cGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACTGCTCTGGTTTTCCAGGCGCTTGAGTTGTTCGACCGTCATGTCGGCCCCCACGCTGGTAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcc                                                                                                                                               <  8:169199/149‑1 (MQ=255)
      gACCTTTGCTGCGGACACCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcca                                                                                                                                              <  4:314556/149‑2 (MQ=255)
        ccTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTccagc                                                                                                                                            <  1:36389/149‑4 (MQ=255)
            ttctgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTccagcaa                                                                                                                                          <  1:369796/145‑6 (MQ=255)
            tgctgcGGACACCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGGCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCAccc                                                                                                                                         <  3:82732/148‑1 (MQ=255)
            tgctgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCAcc                                                                                                                                          <  1:179977/147‑1 (MQ=255)
               tgcGGACAACTCGCTGACCTGAGCAATTCCACGGCGCAGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCa                                                                                                                                     >  6:150720/1‑149 (MQ=255)
                 cGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCa                                                                                                                                     <  1:325246/147‑1 (MQ=255)
                      aaCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCg                                                                                                                              >  1:88969/1‑149 (MQ=255)
                      aaCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCg                                                                                                                              >  4:209286/1‑149 (MQ=255)
                       aCTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGGTTGAGGTGTTCGACCTTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCGTCCCCCAGCACCAGGCGt                                                                                                                             >  8:164320/1‑149 (MQ=255)
                        cTCGCTGACCTGAGCAATTCCCCGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTGTTCGGGCTAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                            <  1:77917/149‑1 (MQ=255)
                        cTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCTTGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                            >  5:303086/1‑149 (MQ=255)
                        cTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGtt                                                                                                                            <  5:346542/149‑1 (MQ=255)
                         tCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTc                                                                                                                           >  5:302464/1‑149 (MQ=255)
                         tCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTc                                                                                                                           >  5:202400/1‑149 (MQ=255)
                                   gAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCgg                                                                                                                 >  1:105670/1‑149 (MQ=255)
                                    aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCggg                                                                                                                >  6:203176/1‑149 (MQ=255)
                                    aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCg                                                                                                                  <  6:20571/147‑1 (MQ=255)
                                    aGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTAGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGGCTTCGGTCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCgg                                                                                                                 <  5:67682/148‑1 (MQ=255)
                                      cAATTCCTCGGCGCTGGTTCTCCAGGCGCTTGAGGTGTTCGACCGTCATGGCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGATGCTTGCTGGAAGCGCGGGGCTTCCTCACGCAGCACCAGGCGTGCGATAGgggggggc                                                                                                              <  8:290746/149‑3 (MQ=255)
                                      cAATTCCACGGCGCTGGTGTTCCAGGCGCGTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGGGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTc                                                                                                              <  8:167927/149‑1 (MQ=255)
                                       aagtCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTAGGCCCCCACGCTGGGAGGTTCCATGTCGTCTTCGGGCGAAGCTTGCTGGAAGCCCCGCGCTTcca                                                                                                                                              <  6:104648/113‑2 (MQ=255)
                                       aaTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCATCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGt                                                                                                               >  3:368471/1‑147 (MQ=255)
                                       aaTTCCACGGCGCTGGTTTTCAAGGCGGTTGAGGTGTTCGACCGTCATGTCGGCCCCCCCGCTGGGAGGTGCCATGGCGTCGTCGGGCGAAGATTGCTGGAAGCCCCGCGCTTcca                                                                                                                                              >  5:104648/1‑115 (MQ=255)
                                        aTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAg                                                                                                            <  5:85984/149‑1 (MQ=255)
                                        aTTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAg                                                                                                            <  6:302464/149‑1 (MQ=255)
                                          tCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGt                                                                                                          >  7:94992/1‑149 (MQ=255)
                                                gcgcTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCgg                                                                                                     >  2:225331/1‑148 (MQ=255)
                                                 cgcTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTa                                                                                                   >  8:277443/1‑149 (MQ=255)
                                                    tGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTc                                                                                                >  4:252112/1‑149 (MQ=255)
                                                          tCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGAGCGGGTCAGGTCGCGGTATTCGATAc                                                                                           <  1:10461/148‑1 (MQ=255)
                                                                 gcTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGccc                                                                                   <  2:105670/149‑1 (MQ=255)
                                                                           ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTAccc                                                                          <  6:9010/148‑1 (MQ=255)
                                                                           ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTAccc                                                                          >  4:259618/1‑148 (MQ=255)
                                                                           ttCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCg                                                                         >  3:78052/1‑149 (MQ=255)
                                                                              gACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCTGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                       <  3:209286/148‑1 (MQ=255)
                                                                              gACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                       <  2:88969/148‑1 (MQ=255)
                                                                               aCCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCAGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAc                                                                       >  6:160633/1‑147 (MQ=255)
                                                                                 cGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgt                                                                  >  2:204895/1‑149 (MQ=255)
                                                                                    cATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAgtcgt                                                               >  3:77659/1‑149 (MQ=255)
                                                                                          ggCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGc                                                         >  2:45192/1‑149 (MQ=255)
                                                                                            cccccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCtt                                                       <  6:303086/149‑1 (MQ=255)
                                                                                              cccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtg                                                     <  3:252112/149‑1 (MQ=255)
                                                                                               ccACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgt                                                    <  1:77065/149‑1 (MQ=255)
                                                                                                cACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                                   >  5:51403/1‑149 (MQ=255)
                                                                                                cACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgt                                                    >  5:37864/1‑148 (MQ=255)
                                                                                                   gCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTTGTCACGGCTTtgttcat                                                >  8:67128/1‑149 (MQ=255)
                                                                                                       ggAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgttcatgtt                                             <  3:259618/148‑1 (MQ=255)
                                                                                                           gTGCCATGGCGTCTTCTGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATgg                                        <  6:37864/149‑1 (MQ=255)
                                                                                                           gTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATgg                                        <  1:225331/149‑1 (MQ=255)
                                                                                                               cATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGACACAGTCGTCACGGCTTTGTTCATGTTCATGGGTa                                     <  1:45192/149‑1 (MQ=255)
                                                                                                                     gTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAg                              <  5:150720/149‑1 (MQ=255)
                                                                                                                           gggCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                                   >  6:64042/1‑122 (MQ=255)
                                                                                                                           gggCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTtgtt                                                   <  5:64042/122‑1 (MQ=255)
                                                                                                                             gCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACgg                      >  6:140834/1‑149 (MQ=255)
                                                                                                                              cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGt                     >  3:344415/1‑149 (MQ=255)
                                                                                                                              cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGC‑TTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt                    >  7:30997/1‑149 (MQ=255)
                                                                                                                              cGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGC‑TTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt                    <  8:30997/149‑1 (MQ=255)
                                                                                                                               gAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGtt                    <  4:344415/149‑1 (MQ=255)
                                                                                                                                       gaggAAGCCCCGCGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGc            >  7:332864/3‑149 (MQ=255)
                                                                                                                                          gAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGt          <  7:277443/148‑1 (MQ=255)
                                                                                                                                           aaGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGttt        <  2:222266/149‑1 (MQ=255)
                                                                                                                                           aaGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGttt        <  6:202400/149‑1 (MQ=255)
                                                                                                                                                 ggggCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGt  >  1:193861/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GACCTCGACCTTTGCTGCGGACAACTCGCTGACCTGAGCAATTCCACGGCGCTGGTTTTCCAGGCGCTTGAGGTGTTCGACCGTCATGTCGGCCCCCACGCTGGGAGGTGCCATGGCGTCTTCGGGCGAAGCTTGCTGGAAGCGCGGGGCTTCCTCACCCAGCACCAGGCGTTCGATAGGGCGGGTCAGGTCGCGGTATTCGATACTGGCGCCCAGGCTACCCGAC‑CAGTCGTCACGGCTTTGTTCATGTTCATGGGTATCCTGAGCTGCACGGTTGTCGATGCGTTTGGCGGT  >  NC_002947/1653991‑1654284

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