Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,303,529 C→T H28Y (CAC→TAC)  yhbW → putative luciferase‑like monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,303,5290CT96.4% 254.7 / 0.1 84H28Y (CAC→TAC) yhbWputative luciferase‑like monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (41/40);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

TCTTTAATGCAGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATA  >  NC_000913/3303402‑3303658
                                                                                                                               |                                                                                                                                 
tctTTAATGCAGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGAt                                                                                                                        >  1:155175/1‑139 (MQ=255)
 ctTTAATGCAGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCAATCTCTCGAtc                                                                                                                       >  2:768846/1‑139 (MQ=255)
        gCAGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCg                                                                                                                <  1:1674573/139‑1 (MQ=255)
          aGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTc                                                                                                              <  1:482846/139‑1 (MQ=255)
            tATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTg                                                                                                            >  1:55272/1‑139 (MQ=255)
            tATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTg                                                                                                            >  1:549001/1‑139 (MQ=255)
              ttAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGc                                                                                                          >  2:365287/1‑139 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCg                                                                                                                    >  2:182335/1‑121 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCg                                                                                                                    <  1:182335/121‑1 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                >  1:603371/1‑125 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                <  2:603371/125‑1 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGt                                                                                                               >  2:1553292/1‑126 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGt                                                                                                               <  1:1553292/126‑1 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAACATTCTCCTACTCTCTCGATCTCGCCCGTc                                                                                                              <  1:779884/127‑1 (MQ=255)
                      tAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAACATTCTCCTACTCTCTCGATCTCGCCCGTc                                                                                                              >  2:779884/1‑127 (MQ=255)
                        aTCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgcgc                                                                                                >  2:584355/1‑139 (MQ=255)
                            gTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCt                                                                                            <  2:55272/139‑1 (MQ=255)
                              aaCAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGAtctc                                                                                                                     >  2:1556035/1‑112 (MQ=255)
                              aaCAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGAtctc                                                                                                                     <  1:1556035/112‑1 (MQ=255)
                                       aGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCg                                                                                                                    >  1:849939/1‑104 (MQ=255)
                                       aGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCg                                                                                                                    <  2:849939/104‑1 (MQ=255)
                                           gTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggc                                                                             <  2:155175/139‑1 (MQ=255)
                                              aTGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          >  1:297454/1‑139 (MQ=255)
                                              aTGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          <  2:297454/139‑1 (MQ=255)
                                                    gATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCGGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          >  2:1395460/1‑133 (MQ=255)
                                                    gATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          <  1:1395460/133‑1 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                >  1:1715830/1‑93 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                <  2:1715830/93‑1 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     >  2:1174532/1‑136 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     <  1:1174532/136‑1 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAacac                                                                  >  1:618779/1‑139 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAacac                                                                  >  2:94162/1‑139 (MQ=255)
                                                      tAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAacac                                                                  >  2:1629175/1‑139 (MQ=255)
                                                          aCCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCAATCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACcaca                                                              <  1:768846/139‑1 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                >  1:1204738/1‑84 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                <  2:1204738/84‑1 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                >  2:430132/1‑84 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCg                                                                                                                <  1:430132/84‑1 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     >  1:570222/1‑127 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     <  2:570222/127‑1 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATg                                                         >  2:977024/1‑139 (MQ=255)
                                                               tGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATg                                                         >  1:156466/1‑139 (MQ=255)
                                                                gCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        <  1:94162/139‑1 (MQ=255)
                                                                gCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        <  1:1629175/139‑1 (MQ=255)
                                                                 cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCCCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggc                                                                             <  1:1190785/117‑1 (MQ=255)
                                                                 cGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCCCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggc                                                                             >  2:1190785/1‑117 (MQ=255)
                                                                                   tGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGtgctgc                                     <  1:977024/139‑1 (MQ=255)
                                                                                    ggCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCa                                                                        >  1:255868/1‑103 (MQ=255)
                                                                                    ggCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCa                                                                        <  2:255868/103‑1 (MQ=255)
                                                                                      cccccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgc                                                                                                  >  2:826335/1‑75 (MQ=255)
                                                                                      cccccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAgc                                                                                                  <  1:826335/75‑1 (MQ=255)
                                                                                      cccccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     >  1:1313466/1‑104 (MQ=255)
                                                                                      cccccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     <  2:1313466/104‑1 (MQ=255)
                                                                                        cccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGt                                <  2:156466/139‑1 (MQ=255)
                                                                                        cccATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGt                                <  1:584355/139‑1 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTc                                                                                                              >  2:517998/1‑58 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTc                                                                                                              <  1:517998/58‑1 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGACAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAat                                                            >  2:468157/1‑108 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGt                                <  2:881829/136‑1 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGt                                >  1:881829/1‑136 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCgg                             >  2:512752/1‑139 (MQ=255)
                                                                                           aTTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCGGGCAGAACACCACAat                                                            <  1:468157/108‑1 (MQ=255)
                                                                                               ccGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        >  2:1282052/1‑108 (MQ=255)
                                                                                               ccGAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        <  1:1282052/108‑1 (MQ=255)
                                                                                                 gAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGc                                                                                                                   >  2:1346560/1‑47 (MQ=255)
                                                                                                 gAAGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGc                                                                                                                   <  1:1346560/47‑1 (MQ=255)
                                                                                                  aaGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatcat                                                                                       >  1:1204925/1‑74 (MQ=255)
                                                                                                  aaGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatcat                                                                                       <  2:1204925/74‑1 (MQ=255)
                                                                                                  aaGGTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGAt                      <  2:618779/139‑1 (MQ=255)
                                                                                                    ggTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGAt                                                                                                                        >  2:413578/1‑39 (MQ=38)
                                                                                                    ggTTCTTCAGCGCGAGAAGCATTCTCCTACTCTCTCGAt                                                                                                                        <  1:413578/39‑1 (MQ=255)
                                                                                                         ttcAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatcat                                                                                       <  2:431733/67‑1 (MQ=255)
                                                                                                         ttcAGCGCGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTatcat                                                                                       >  1:431733/1‑67 (MQ=255)
                                                                                                             gcgcGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          <  2:70481/76‑1 (MQ=255)
                                                                                                             gcgcGAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTActggctgg                                                                          >  1:70481/1‑76 (MQ=255)
                                                                                                              cgcgAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     >  2:111873/1‑80 (MQ=255)
                                                                                                              cgcgAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGaa                                                                     <  1:111873/80‑1 (MQ=255)
                                                                                                               gcgAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACg                                 >  2:629834/1‑115 (MQ=255)
                                                                                                               gcgAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACg                                 <  1:629834/115‑1 (MQ=255)
                                                                                                               gcgAGAAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTggcg         <  1:365287/139‑1 (MQ=255)
                                                                                                                  agaAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGgcggcg      >  1:704028/1‑139 (MQ=255)
                                                                                                                   gaAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        <  1:1082848/88‑1 (MQ=255)
                                                                                                                   gaAGCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGa                                                        >  2:1082848/1‑88 (MQ=255)
                                                                                                                      gCATTCTCCTACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATa  <  2:704028/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
TCTTTAATGCAGTATTAAAGATTAATCGGTAACAAAGTGAGCTGTTATGACTGATAAAACCATTGCGTTTTCGCTACTCGATCTGGCCCCCATTCCCGAAGGTTCTTCAGCGCGAGAAGCATTCTCCCACTCTCTCGATCTCGCCCGTCTGGCTGAAAAGCGCGGCTATCATCGCTACTGGCTGGCAGAACACCACAATATGACTGGCATTGCCAGTGCTGCCACGTCGGTATTGATCGGCTATCTGGCGGCGAATA  >  NC_000913/3303402‑3303658

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