Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 165,932 A→C 8.0% D325E (GAT→GAG pntAA ← pyridine nucleotide transhydrogenase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947165,9320AC 8.0% 274.5 / 11.1 96D325E (GAT→GAGpntAApyridine nucleotide transhydrogenase subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (47/41);  new base C (4/4);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTC  >  NC_002947/165784‑166070
                                                                                                                                                    |                                                                                                                                          
cTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCAt                                                                                                                                           >  1:1218012/1‑150 (MQ=255)
cTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCAt                                                                                                                                           >  1:338173/1‑150 (MQ=255)
  tAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCg                                                                                                                                         <  2:1218012/150‑1 (MQ=255)
     ccGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACCGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                      <  2:190385/150‑1 (MQ=255)
     ccGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                      <  1:837926/150‑1 (MQ=255)
     ccGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCa                                                                                                                                      <  1:2008207/150‑1 (MQ=255)
       gTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTAGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACt                                                                                                                                    >  2:1735530/1‑150 (MQ=255)
             gCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTg                                                                                                                              <  2:1505485/150‑1 (MQ=255)
              cGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTgg                                                                                                                             <  1:1060197/150‑1 (MQ=255)
                   gacCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGCGCCgg                                                                                                                        >  1:2288374/1‑150 (MQ=255)
                    acCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                        >  1:1563608/1‑149 (MQ=255)
                    acCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCgg                                                                                                                        <  2:1563608/149‑1 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                      >  1:1717092/1‑150 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                      <  2:2099533/150‑1 (MQ=255)
                     ccTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGCGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCa                                                                                                                      <  1:2297412/150‑1 (MQ=255)
                      cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACtg                                                                                                                                   <  2:71848/136‑2 (MQ=255)
                      cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACtg                                                                                                                                   >  1:71848/1‑135 (MQ=255)
                      cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                     >  2:1426567/1‑150 (MQ=255)
                      cTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAg                                                                                                                     >  2:435801/1‑150 (MQ=255)
                       tGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAgg                                                                                                                    >  1:1920353/1‑150 (MQ=255)
                        gagcATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGt                                                                                                                   <  1:1062188/147‑1 (MQ=255)
                        gACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGt                                                                                                                   <  1:1544285/150‑1 (MQ=255)
                          ccATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTTCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTg                                                                                                                 <  1:1024841/150‑1 (MQ=255)
                           cATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGc                                                                                                                                       >  1:1887315/1‑127 (MQ=255)
                           cATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGc                                                                                                                                       <  2:1887315/127‑1 (MQ=255)
                             tCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTc                                                                                                              >  1:1952712/1‑150 (MQ=255)
                               aCGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCgg                                                                                                            <  2:513701/150‑1 (MQ=255)
                               aCGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCgg                                                                                                            >  1:327529/1‑150 (MQ=255)
                                cGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCggg                                                                                                           <  2:327529/150‑1 (MQ=255)
                                 ggCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGc                                                                                                          <  1:435801/150‑1 (MQ=255)
                                    acaTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCAc                                                                                                                                     <  1:1981880/120‑1 (MQ=255)
                                    acaTCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCAc                                                                                                                                     >  2:1981880/1‑120 (MQ=255)
                                          ggcaggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatga                                                                                                 >  1:1193473/1‑150 (MQ=255)
                                            caggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATc                                                                                               <  2:1952712/150‑1 (MQ=255)
                                             aggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGc                                                                                                          >  2:594775/1‑138 (MQ=255)
                                             aggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGc                                                                                                          <  1:594775/138‑1 (MQ=255)
                                             aggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCa                                                                                              >  2:969136/1‑150 (MQ=255)
                                              ggcCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCAc                                                                                             <  1:969136/150‑1 (MQ=255)
                                                    gactaTGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGTTCGGGCCGACGATGATCCCGCCGtt                                                                                       >  2:754535/1‑150 (MQ=255)
                                                        aTGTCGTCTTCGAGGTTGAAGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTcc                                                                                   <  2:1590455/150‑1 (MQ=255)
                                                          gtcgtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCTGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTg                                                                                 >  1:701198/1‑150 (MQ=255)
                                                             gtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCgg                                                                                                                                        >  1:505473/1‑92 (MQ=255)
                                                             gtcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCgg                                                                                                                                        <  2:505473/92‑1 (MQ=255)
                                                              tcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCa                                                                                              >  2:2098578/1‑133 (MQ=255)
                                                              tcTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCa                                                                                              <  1:2098578/133‑1 (MQ=255)
                                                                ttCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGcaccac                                                                           >  1:801661/1‑150 (MQ=255)
                                                                 tCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCaccacc                                                                          >  2:1411853/1‑150 (MQ=255)
                                                                    aGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTGCAGCAGGTTGCGTGCGGACCGTGCCGAGGCATCGGCCCTCCCCTGGGCCGGCAGGTTGGTCGGGGCGACGAGGATCACGCCGTTCTCCTGCACCACCTgg                                                                       >  2:320914/1‑150 (MQ=255)
                                                                    aGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                <  1:1862282/141‑1 (MQ=255)
                                                                    aGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGc                                                                                >  2:1862282/1‑141 (MQ=255)
                                                                      gTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGCGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTc                                                                     >  1:1994153/1‑150 (MQ=255)
                                                                       ttGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCt                                                                    >  2:1264630/1‑150 (MQ=255)
                                                                       ttGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGCGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCt                                                                    <  2:1994153/150‑1 (MQ=255)
                                                                            gatCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGAc                                                               <  2:386591/150‑1 (MQ=255)
                                                                             ttccgcgcgcCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACt                                                                                                                                    <  1:1912855/76‑1 (MQ=255)
                                                                             ttccgcgcgcCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGATGCCTCGGCACt                                                                                                                                    >  2:1912855/5‑80 (MQ=255)
                                                                             atCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACg                                                              >  2:1414399/1‑150 (MQ=255)
                                                                              tCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACgg                                                             >  1:554690/1‑150 (MQ=255)
                                                                                aGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTc                                                           >  1:1383678/1‑150 (MQ=255)
                                                                                     gcCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCgg                                                      >  2:70426/1‑150 (MQ=255)
                                                                                       cATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGc                                                    >  1:796416/1‑150 (MQ=255)
                                                                                       cATCCTTGGCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGc                                                    >  2:93904/1‑150 (MQ=255)
                                                                                            ttGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                <  2:815217/149‑1 (MQ=255)
                                                                                            ttGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                >  1:815217/1‑149 (MQ=255)
                                                                                            ttGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTg                                               >  1:866541/1‑150 (MQ=255)
                                                                                              gTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgcc                                             <  2:866979/150‑1 (MQ=255)
                                                                                              gTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgcc                                             <  2:866541/150‑1 (MQ=255)
                                                                                               tCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgccg                                            >  1:1294379/1‑150 (MQ=255)
                                                                                                cGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                >  2:1279991/1‑145 (MQ=255)
                                                                                                cGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGtt                                                <  1:1279991/145‑1 (MQ=255)
                                                                                                cGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTgccgc                                           >  2:1753796/1‑150 (MQ=255)
                                                                                                 gAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCg                                                                                                             >  1:911538/1‑83 (MQ=255)
                                                                                                 gAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCg                                                                                                             <  2:911538/83‑1 (MQ=255)
                                                                                                   aCAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCtt                                        >  2:1279747/1‑150 (MQ=255)
                                                                                                     agcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTgg                                      <  2:1383678/150‑1 (MQ=255)
                                                                                                      gcagcTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTggg                                     <  1:1279747/150‑1 (MQ=255)
                                                                                                            tCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAg                               <  1:93904/150‑1 (MQ=255)
                                                                                                              aTGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTTCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGcc                             >  1:1574151/1‑150 (MQ=255)
                                                                                                              aTGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGcc                             >  2:5848/1‑150 (MQ=255)
                                                                                                              aTGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCAACTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGcc                             >  2:278743/1‑150 (MQ=255)
                                                                                                               tGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatga                                                                                                 >  2:1963217/1‑81 (MQ=255)
                                                                                                               tGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACgatga                                                                                                 <  1:1963217/81‑1 (MQ=255)
                                                                                                                gAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAg                           >  1:1285166/1‑150 (MQ=255)
                                                                                                                gAAGTCCAGCAGGTTGCGTGCATACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGCTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAg                           >  2:598718/1‑150 (MQ=255)
                                                                                                                 aaGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAgg                          >  2:1776294/1‑150 (MQ=255)
                                                                                                                 aaGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAgg                          >  1:90515/1‑150 (MQ=255)
                                                                                                                  aGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCTGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGt                         <  2:701198/150‑1 (MQ=255)
                                                                                                                       agcagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACAATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATg                    <  2:2029654/150‑1 (MQ=255)
                                                                                                                         cagGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGac                  <  2:1285166/150‑1 (MQ=255)
                                                                                                                            gTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGaccac               <  1:5848/150‑1 (MQ=255)
                                                                                                                            gTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGaccac               >  2:2085065/1‑150 (MQ=255)
                                                                                                                                 gtgcgtACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                       >  2:1220159/1‑121 (MQ=255)
                                                                                                                                 gtgcgtACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTg                                       <  1:1220159/121‑1 (MQ=255)
                                                                                                                                        cAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCtt   <  1:1776294/150‑1 (MQ=255)
                                                                                                                                         aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTGTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCTTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTc  <  1:1748275/150‑1 (MQ=255)
                                                                                                                                         aGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTc  <  1:1414399/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
CTTAGCCGTTTTTGCGGACGACCTGACCATCACGGCACATCAGGCAGGCCGCGACGATGTCGTCTTCGAGGTTGATGATCAGCGCGCCATCCTTGTCGAACAGCAGCTTCATGAAGTCCAGCAGGTTGCGTGCGTACAGTGCCGAGGCATCGGCACTCACCTGGGCCGGCAGGTTGGTCGGGCCGACGATGATCACGCCGTTCTCCTGCACCACCTGGTCTGCGACGGTCAGCGGGCAGTTGCCGCCTTGGGCTGCAGCCAGGTCGATGACCACCGAGCCGGGCTTC  >  NC_002947/165784‑166070

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