Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 312,189 C→T G26G (GGC→GGT ykgL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913312,1890CT98.6% 230.8 / NA 72G26G (GGC→GGTykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/38);  minor base G (1/0);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312068‑312323
                                                                                                                         |                                                                                                                                      
ccAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  1:1256784/139‑1 (MQ=255)
ccAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  1:642132/139‑1 (MQ=255)
 cAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  2:1634259/139‑1 (MQ=255)
   aCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  1:1030311/1‑136 (MQ=255)
   aCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:850755/137‑1 (MQ=255)
   aCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:850755/1‑137 (MQ=255)
   aCAGATTTACCACATTTATCATCTTAACAGTGTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  2:1030311/136‑1 (MQ=255)
        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               >  1:105489/1‑139 (MQ=255)
        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               <  2:572805/139‑1 (MQ=255)
        tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               <  2:5831/139‑1 (MQ=255)
           aCCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGa                                                                                                            >  1:1527479/1‑139 (MQ=255)
             cacaTTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          <  2:1527479/139‑1 (MQ=255)
                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaac                                                                                                                           >  1:1108938/1‑109 (MQ=255)
                          ttAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaac                                                                                                                           <  2:1108938/109‑1 (MQ=255)
                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               <  1:1258017/119‑1 (MQ=255)
                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               >  2:1258017/1‑119 (MQ=255)
                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  1:1428798/1‑139 (MQ=255)
                            aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  1:220876/1‑139 (MQ=255)
                              cAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                         >  2:261609/1‑139 (MQ=255)
                                 tttGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                                      <  2:220876/139‑1 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAATGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  2:1239719/105‑1 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  2:426822/1‑104 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  1:426822/104‑1 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  1:1239719/1‑105 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:1286882/1‑105 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:1286882/105‑1 (MQ=255)
                                   tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  1:776191/1‑139 (MQ=255)
                                      gTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          >  2:1429798/1‑114 (MQ=255)
                                      gTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          <  1:1429798/114‑1 (MQ=255)
                                           aaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAgg                                                                               >  2:566025/1‑136 (MQ=255)
                                           aaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAgg                                                                               <  1:566025/136‑1 (MQ=255)
                                               cAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaa                                                                                       >  1:1252496/1‑124 (MQ=255)
                                               cAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaa                                                                                       <  2:1252496/124‑1 (MQ=255)
                                                   cTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTc                                                                    >  1:1727854/1‑139 (MQ=255)
                                                   cTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTc                                                                    >  2:389541/1‑139 (MQ=255)
                                                    tGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGc                                                                                                                                   >  1:1221211/1‑75 (MQ=255)
                                                    tGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGc                                                                                                                                   <  2:1221211/75‑1 (MQ=255)
                                                     gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCg                                                                                              <  1:1385178/111‑1 (MQ=255)
                                                     gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCg                                                                                              >  2:1385178/1‑111 (MQ=255)
                                                     gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaa                                                                  <  1:389541/139‑1 (MQ=255)
                                                      cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:390786/1‑86 (MQ=255)
                                                      cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:390786/86‑1 (MQ=255)
                                                         cgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGa                                                              >  2:1491287/1‑139 (MQ=255)
                                                          gcAAAGGGTGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             <  2:1054859/139‑1 (MQ=255)
                                                                ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         <  2:69798/137‑1 (MQ=255)
                                                                ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         >  1:69798/1‑137 (MQ=255)
                                                                ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaa                                                       <  1:261609/139‑1 (MQ=255)
                                                                ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaa                                                       <  2:105489/139‑1 (MQ=255)
                                                                   ggagACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaagaa                                                    >  2:670364/1‑139 (MQ=255)
                                                                             tGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaa                                                                  <  2:30675/115‑1 (MQ=255)
                                                                             tGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaa                                                                  >  1:30675/1‑115 (MQ=255)
                                                                              ggCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAg                                         <  1:1221669/139‑1 (MQ=255)
                                                                              ggCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAg                                         <  1:1491287/139‑1 (MQ=255)
                                                                                cGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             >  2:702376/1‑117 (MQ=255)
                                                                                cGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             <  1:702376/117‑1 (MQ=255)
                                                                                         tCTAAAGAAGCGTTGATGCCAGCCTTCGGAGGGGAGGCGGTTCAACAAGTCAGCAGTAAGATATTCGAATATTCGGAGACCAACGTCAGGCAGTGTGCATCAAATGAACAGGAAGAATATATCCACAGGAAAATTAATg                              >  2:1221669/1‑139 (MQ=255)
                                                                                         tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGTATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATg                              >  2:1011325/1‑139 (MQ=255)
                                                                                         tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATg                              >  1:150530/1‑139 (MQ=255)
                                                                                              aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAg                                                                            >  1:697007/1‑88 (MQ=255)
                                                                                              aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAg                                                                            <  2:697007/88‑1 (MQ=255)
                                                                                                   ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  1:684307/1‑41 (MQ=255)
                                                                                                   ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  2:684307/41‑1 (MQ=255)
                                                                                                     ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAtt                                                                                                                <  2:603849/45‑1 (MQ=255)
                                                                                                     ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAtt                                                                                                                >  1:603849/1‑45 (MQ=255)
                                                                                                      tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  2:150530/139‑1 (MQ=255)
                                                                                                      tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  2:1428798/139‑1 (MQ=255)
                                                                                                           ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        >  2:292294/1‑95 (MQ=255)
                                                                                                           ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        <  1:292294/95‑1 (MQ=255)
                                                                                                           ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac            <  1:670364/139‑1 (MQ=255)
                                                                                                                     ggggTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAt  <  1:800545/139‑1 (MQ=255)
                                                                                                                      gggTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtata                   <  1:981089/121‑1 (MQ=255)
                                                                                                                      gggTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtata                   >  2:981089/1‑121 (MQ=255)
                                                                                                                         |                                                                                                                                      
CCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312068‑312323

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