Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,709,561 (C)6→7 coding (915/1692 nt) eptB ← KDO phosphoethanolamine transferase, Ca(2+)‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,709,5551.C95.7% 234.6 / 2.0 69coding (921/1692 nt)eptBKDO phosphoethanolamine transferase, Ca(2+)‑inducible
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base C (37/29);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCC  >  NC_000913/3709425‑3709677
                                                                                                                                    |                                                                                                                          
tgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccct                                                                                                                    <  2:770449/139‑3 (MQ=255)
       cATCGTGTTGCTGTAGAACCACATTTCGCTCGGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTAc                                                                                                              <  2:368128/139‑1 (MQ=255)
                   gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  1:1024893/1‑139 (MQ=255)
                           acaTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTgaga                                                                                          >  1:1083036/1‑139 (MQ=255)
                            caTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         >  1:378148/1‑139 (MQ=255)
                             aTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   >  2:1231936/1‑108 (MQ=255)
                             aTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   <  1:1231936/111‑4 (MQ=255)
                                    tctGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAGAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCCGTGAGAGTTTGGTTg                                                                                 <  2:1024893/139‑1 (MQ=255)
                                         aTAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTa                                                                            <  2:751455/139‑1 (MQ=255)
                                           ctgaagcagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  2:868696/8‑139 (MQ=255)
                                                agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccct                                                                                                                    >  2:466111/1‑89 (MQ=255)
                                                agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccct                                                                                                                    <  1:466111/91‑3 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgc‑g                                                                                                                             <  2:346701/81‑1 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgc‑g                                                                                                                             >  1:346701/1‑81 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAg                                                                   >  1:1502155/1‑139 (MQ=255)
                                                   ggTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                  <  1:1151908/139‑1 (MQ=255)
                                                       aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAAcc                                                              >  1:1254559/1‑139 (MQ=255)
                                                       aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAAcc                                                              >  2:762465/1‑139 (MQ=255)
                                                         aaCTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCAc                                                            >  1:400033/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           <  1:783797/103‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           >  2:783797/1‑103 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       >  1:247903/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       <  1:868696/139‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       >  2:740026/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑CCCCCCTg                                                                                                                    <  1:5416/77‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑CCCCCCTg                                                                                                                    >  2:5416/1‑77 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                  <  2:1254559/139‑1 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                  >  1:559240/1‑139 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   >  1:1106449/1‑70 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   <  2:1106449/73‑4 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         <  2:949715/95‑1 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   >  1:259938/1‑65 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccctg                                                                                                                   <  2:259938/68‑4 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCATTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         >  1:949715/1‑95 (MQ=255)
                                                                         cTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑cccc                                                                                                                        <  1:344168/62‑1 (MQ=255)
                                                                         cTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑cccc                                                                                                                        >  2:344168/1‑62 (MQ=255)
                                                                           tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccc                                                                                                                         <  1:333423/59‑1 (MQ=255)
                                                                           tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccc                                                                                                                         >  2:333423/1‑59 (MQ=255)
                                                                                aCCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     <  2:378148/139‑1 (MQ=255)
                                                                                 ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgc‑g                                                                                                                             <  1:706478/50‑1 (MQ=255)
                                                                                 ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgc‑g                                                                                                                             >  2:706478/1‑50 (MQ=255)
                                                                                   gcgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                   <  2:640583/90‑1 (MQ=255)
                                                                                   gcgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                   >  1:640583/1‑90 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCaa                                                                                                    <  1:820454/72‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCaa                                                                                                    >  2:820454/1‑72 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACg                                                                                                  <  2:1266046/74‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  1:1266046/1‑74 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         <  1:565379/99‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         >  2:565379/1‑99 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  2:98826/1‑98 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          <  1:98826/98‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     <  1:1195238/135‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     >  2:1195238/1‑135 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTggg                                 >  1:1326564/1‑139 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccct                                                                                                                    <  1:1510755/55‑3 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccct                                                                                                                    >  2:1510755/1‑53 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              >  1:471276/1‑89 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              <  2:471276/89‑1 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccc                                                                                                                     <  2:464859/49‑2 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑ccccccc                                                                                                                     >  1:464859/1‑48 (MQ=255)
                                                                                                gttctTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                     <  2:1502155/139‑1 (MQ=255)
                                                                                                gttctTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                     <  2:247903/139‑1 (MQ=255)
                                                                                                 ttctTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  1:989073/1‑85 (MQ=255)
                                                                                                 ttctTTTAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGCATc                                                                          <  2:989073/85‑1 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  2:914504/76‑1 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  1:914504/1‑76 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              >  2:720233/1‑78 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              <  1:720233/78‑1 (MQ=255)
                                                                                                      ttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCaa                                                                                                    >  2:547300/1‑54 (MQ=255)
                                                                                                      ttAATGTGCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCaa                                                                                                    <  1:547300/54‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACg                                                                                                                 <  1:873390/36‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACg                                                                                                                 >  2:873390/1‑36 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAAc                                                                                                          >  2:422373/1‑43 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAAc                                                                                                          <  1:422373/43‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  1:1218498/69‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  2:1218498/1‑69 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGc          >  2:363304/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           >  1:739032/1‑55 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           <  2:739032/55‑1 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:639926/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:1368132/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:553692/1‑139 (MQ=255)
                                                                                                                  gCGGATTATCTTCCGCGG‑CCCCCCTGACGTACAAACGTGCAACGAAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAAATTTTTCGCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGc   <  2:1083036/139‑1 (MQ=255)
                                                                                                                   cGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  1:488037/1‑61 (MQ=255)
                                                                                                                   cGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  2:488037/61‑1 (MQ=255)
                                                                                                                   cGGATTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTAGCACAGGAGTAACCACGGAACGCCGCCAGACTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGcc  >  1:751277/1‑139 (MQ=255)
                                                                                                                      aTTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACg                                                                                                  <  1:1209936/40‑1 (MQ=255)
                                                                                                                      aTTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  2:1209936/1‑40 (MQ=255)
                                                                                                                      aTTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         <  1:1504998/65‑1 (MQ=255)
                                                                                                                      aTTATCTTCCGC‑GCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         >  2:1504998/1‑65 (MQ=255)
                                                                                                                                    |                                                                                                                          
TGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGC‑G‑CCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCC  >  NC_000913/3709425‑3709677

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑