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

TGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCC  >  NC_000913/3709425‑3709677
                                                                                                                                    |                                                                                                                          
tgtCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                    <  2:770449/139‑3 (MQ=255)
       cATCGTGTTGCTGTAGAACCACATTTCGCTCGGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTAc                                                                                                              <  2:368128/139‑1 (MQ=255)
                   gTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  1:1024893/1‑139 (MQ=255)
                           acaTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTgaga                                                                                          >  1:1083036/1‑139 (MQ=255)
                            caTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         >  1:378148/1‑139 (MQ=255)
                             aTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   >  2:1231936/1‑108 (MQ=255)
                             aTTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   <  1:1231936/111‑4 (MQ=255)
                                    tctGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAGAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCCGTGAGAGTTTGGTTg                                                                                 <  2:1024893/139‑1 (MQ=255)
                                         aTAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTa                                                                            <  2:751455/139‑1 (MQ=255)
                                           ctgaagcagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  2:868696/8‑139 (MQ=255)
                                                agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                    >  2:466111/1‑89 (MQ=255)
                                                agagGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                    <  1:466111/91‑3 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgcg                                                                                                                             <  2:346701/81‑1 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgcg                                                                                                                             >  1:346701/1‑81 (MQ=255)
                                                  agGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAg                                                                   >  1:1502155/1‑139 (MQ=255)
                                                   ggTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGt                                                                  <  1:1151908/139‑1 (MQ=255)
                                                       aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAAcc                                                              >  1:1254559/1‑139 (MQ=255)
                                                       aGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAAcc                                                              >  2:762465/1‑139 (MQ=255)
                                                         aaCTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCAc                                                            >  1:400033/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           <  1:783797/103‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           >  2:783797/1‑103 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       >  1:247903/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       <  1:868696/139‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAc                                                       >  2:740026/1‑139 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTg                                                                                                                    <  1:5416/77‑1 (MQ=255)
                                                              aaTCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTg                                                                                                                    >  2:5416/1‑77 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                  <  2:1254559/139‑1 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAAcgccgc                                                  >  1:559240/1‑139 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   >  1:1106449/1‑70 (MQ=255)
                                                                   tAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   <  2:1106449/73‑4 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         <  2:949715/95‑1 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   >  1:259938/1‑65 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccctg                                                                                                                   <  2:259938/68‑4 (MQ=255)
                                                                        gCTTCAGAACCGCGAAAATGTTCTGTTCATTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGagag                                                                                         >  1:949715/1‑95 (MQ=255)
                                                                         cTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGcccc                                                                                                                        <  1:344168/62‑1 (MQ=255)
                                                                         cTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGcccc                                                                                                                        >  2:344168/1‑62 (MQ=255)
                                                                           tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccc                                                                                                                         <  1:333423/59‑1 (MQ=255)
                                                                           tCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccc                                                                                                                         >  2:333423/1‑59 (MQ=255)
                                                                                aCCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     <  2:378148/139‑1 (MQ=255)
                                                                                 ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgcg                                                                                                                             <  1:706478/50‑1 (MQ=255)
                                                                                 ccGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCcgcg                                                                                                                             >  2:706478/1‑50 (MQ=255)
                                                                                   gcgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                   <  2:640583/90‑1 (MQ=255)
                                                                                   gcgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGt                                                                                   >  1:640583/1‑90 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                    <  1:820454/72‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                    >  2:820454/1‑72 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                  <  2:1266046/74‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  1:1266046/1‑74 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         <  1:565379/99‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         >  2:565379/1‑99 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  2:98826/1‑98 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          <  1:98826/98‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     <  1:1195238/135‑1 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTcc                                     >  2:1195238/1‑135 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTggg                                 >  1:1326564/1‑139 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                    <  1:1510755/55‑3 (MQ=255)
                                                                                    cgAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccct                                                                                                                    >  2:1510755/1‑53 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              >  1:471276/1‑89 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              <  2:471276/89‑1 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccc                                                                                                                     <  2:464859/49‑2 (MQ=255)
                                                                                         aTGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGccccccc                                                                                                                     >  1:464859/1‑48 (MQ=255)
                                                                                                gttctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                     <  2:1502155/139‑1 (MQ=255)
                                                                                                gttctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCg                     <  2:247903/139‑1 (MQ=255)
                                                                                                 ttctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATc                                                                          >  1:989073/1‑85 (MQ=255)
                                                                                                 ttctTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGCATc                                                                          <  2:989073/85‑1 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  2:914504/76‑1 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  1:914504/1‑76 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              >  2:720233/1‑78 (MQ=255)
                                                                                                    ttttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCgg                                                                              <  1:720233/78‑1 (MQ=255)
                                                                                                      ttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                    >  2:547300/1‑54 (MQ=255)
                                                                                                      ttAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCaa                                                                                                    <  1:547300/54‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACg                                                                                                                 <  1:873390/36‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACg                                                                                                                 >  2:873390/1‑36 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAAc                                                                                                          >  2:422373/1‑43 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAAc                                                                                                          <  1:422373/43‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  1:1218498/69‑1 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  2:1218498/1‑69 (MQ=255)
                                                                                                           gtgCGCTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGc          >  2:363304/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           >  1:739032/1‑55 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTgag                                                                                           <  2:739032/55‑1 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:639926/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:1368132/1‑139 (MQ=255)
                                                                                                              cgcTGCGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCa       >  1:553692/1‑139 (MQ=255)
                                                                                                                  gCGGATTATCTTCCGCGGCCCCCCTGACGTACAAACGTGCAACGAAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAAATTTTTCGCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGc   <  2:1083036/139‑1 (MQ=255)
                                                                                                                   cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                >  1:488037/1‑61 (MQ=255)
                                                                                                                   cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGc                                                                                <  2:488037/61‑1 (MQ=255)
                                                                                                                   cGGATTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTAGCACAGGAGTAACCACGGAACGCCGCCAGACTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGcc  >  1:751277/1‑139 (MQ=255)
                                                                                                                      aTTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                  <  1:1209936/40‑1 (MQ=255)
                                                                                                                      aTTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACg                                                                                                  >  2:1209936/1‑40 (MQ=255)
                                                                                                                      aTTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         <  1:1504998/65‑1 (MQ=255)
                                                                                                                      aTTATCTTCCGCGCCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATca                                                                         >  2:1504998/1‑65 (MQ=255)
                                                                                                                                    |                                                                                                                          
TGTCCGCCATCGTGTTGCTGTAGAACCACATTTCGCTCTGCATAGCGTAGAGGTCAGAACTGAATCCTAACTGCTTCAGAACCGCGAAAATGTTCTGTTCTTTTAATGTGCGCTGCGGATTATCTTCCGCGCCCCCCTGACGTACAAACATGCAACGCAGTGAGAGTTTGGTTGCGGTATCACAGGAGTAACCACGGAACGCCGCCAGATTTTTCTCCTGGGCCAGTTTCGGCGTGGTATTACGCTCATAGCC  >  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: