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,1890CT99.0% 299.1 / NA 96G26G (GGC→GGTykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (50/45);  minor base A (1/0);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312063‑312323
                                                                                                                              |                                                                                                                                      
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGATGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                            <  1:1388984/139‑1 (MQ=255)
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                            >  2:1639302/1‑139 (MQ=255)
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                            <  1:1361820/139‑1 (MQ=255)
  gTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaaca                                                                                                                          >  1:897907/1‑139 (MQ=255)
  gTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAAAGGGTTGATGCCAGCATAAAGGGCAGATGCTGGTcaaca                                                                                                                          >  2:656123/1‑139 (MQ=255)
     ccAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  1:946601/139‑1 (MQ=255)
       aaCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCa                                                                                                                     >  1:1305800/1‑139 (MQ=255)
        aCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCaa                                                                                                                    >  1:888077/1‑139 (MQ=255)
             tttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               >  2:798792/1‑139 (MQ=255)
                     aTTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATg                                                                                                       >  2:1607405/1‑139 (MQ=255)
                      tttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATgg                                                                                                      >  1:362829/1‑139 (MQ=255)
                      tttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATgg                                                                                                      >  2:425180/1‑139 (MQ=255)
                       ttATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGa                                                                                                     <  2:845292/139‑1 (MQ=255)
                           catcTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAatat                                                                                                 <  1:798792/139‑1 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:1409150/112‑1 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:1409150/1‑112 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTaa                                                                                                              <  1:205365/120‑1 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTaa                                                                                                              >  2:205365/1‑120 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  2:1294612/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGga                                                                                            <  2:1205918/138‑1 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGga                                                                                            >  1:1205918/1‑138 (MQ=255)
                                     gTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGa                                                                                                            >  1:132186/1‑118 (MQ=255)
                                     gTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGa                                                                                                            <  2:132186/118‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:393247/1‑105 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:393247/105‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCg                                                                                              >  2:610022/1‑129 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCg                                                                                              <  1:610022/129‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  1:1089453/1‑139 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  2:1291639/1‑139 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  1:1118232/1‑139 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  1:276515/1‑139 (MQ=255)
                                           gTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                                 >  2:696799/1‑139 (MQ=255)
                                              aaaaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGc                                                                              >  1:62903/1‑139 (MQ=255)
                                                    cAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTgg                                                                        >  2:59544/1‑139 (MQ=255)
                                                          gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          <  2:119758/99‑1 (MQ=255)
                                                          gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          >  1:119758/1‑99 (MQ=255)
                                                          gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaa                                                                  <  1:553675/139‑1 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          <  2:188878/98‑1 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          >  1:188878/1‑98 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAAt                                                                                                        <  1:1148245/100‑1 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAAt                                                                                                        >  2:1148245/1‑100 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaaa                                                                 >  2:742330/1‑139 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGTTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAACCGTCCGGCAGTTGGCCTCAAATGAAc                                                            <  2:438149/139‑1 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  2:1118232/139‑1 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  2:1305800/139‑1 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  2:62903/139‑1 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  1:696799/139‑1 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGa                                                                                                                                      <  1:847316/62‑1 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGa                                                                                                                                      >  2:847316/1‑62 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         <  1:59544/139‑1 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         <  2:1089453/139‑1 (MQ=255)
                                                                    gggAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        <  2:362829/139‑1 (MQ=255)
                                                                     ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  2:1117989/75‑1 (MQ=255)
                                                                     ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  1:1117989/1‑75 (MQ=255)
                                                                            aCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatata                                                >  2:159393/1‑139 (MQ=255)
                                                                             cTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTc                                                                                               >  2:307809/1‑91 (MQ=255)
                                                                             cTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTc                                                                                               <  1:307809/91‑1 (MQ=255)
                                                                                aGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATccc                                            >  1:498767/1‑139 (MQ=255)
                                                                                    gCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAgg                                        >  1:1421328/1‑139 (MQ=255)
                                                                                         ttCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                         <  1:284051/85‑1 (MQ=255)
                                                                                         ttCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                         >  2:284051/1‑85 (MQ=255)
                                                                                          tCAGTCTAAAATAGGGGTGATGCCAGCATAAAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAAtt                                  >  2:270704/1‑139 (MQ=255)
                                                                                                 aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:888077/139‑1 (MQ=255)
                                                                                                 aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:1421328/139‑1 (MQ=255)
                                                                                                 aaaaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAAAAAGTTAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:217875/139‑1 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        >  1:1026503/1‑108 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        <  2:1026503/108‑1 (MQ=255)
                                                                                                        ggctGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       <  1:566789/85‑1 (MQ=255)
                                                                                                        ggctGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       >  2:566789/4‑88 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          <  2:707528/53‑1 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATa                                                                                                          >  1:707528/1‑53 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       <  2:190364/88‑1 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       >  2:898569/1‑88 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       <  1:898569/88‑1 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGc                                                                       >  1:190364/1‑88 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtat                    >  1:1327785/1‑139 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtat                    >  2:1217302/1‑139 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtat                    >  1:586961/1‑139 (MQ=255)
                                                                                                          ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaa                                                       >  2:970088/1‑102 (MQ=255)
                                                                                                          ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaa                                                       <  1:970088/102‑1 (MQ=255)
                                                                                                          ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaa                  >  2:1516347/1‑139 (MQ=255)
                                                                                                          ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaa                  >  1:1463004/1‑139 (MQ=255)
                                                                                                          ttGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaa                  >  1:411806/1‑139 (MQ=255)
                                                                                                           tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  1:1639302/139‑1 (MQ=255)
                                                                                                           tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  1:1607405/139‑1 (MQ=255)
                                                                                                            gATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                                 <  2:1236730/74‑1 (MQ=255)
                                                                                                            gATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCa                                                                                 >  1:1236730/1‑74 (MQ=255)
                                                                                                             aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           >  1:791335/1‑127 (MQ=255)
                                                                                                             aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:791335/127‑1 (MQ=255)
                                                                                                             aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAaca               <  2:897907/139‑1 (MQ=255)
                                                                                                                ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac            <  1:1217302/139‑1 (MQ=255)
                                                                                                                 cAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACAcc           <  1:159393/139‑1 (MQ=255)
                                                                                                                     aTACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACa       <  1:1291639/139‑1 (MQ=255)
                                                                                                                       aCAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           >  1:317198/1‑117 (MQ=255)
                                                                                                                       aCAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:317198/117‑1 (MQ=255)
                                                                                                                          ggggTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAt  >  1:819503/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                      
TAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312063‑312323

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