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,1890CT97.9% 291.5 / ‑2.4 94G26G (GGC→GGTykgLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (45/47);  total (46/48)
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.63e-01

TAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312063‑312323
                                                                                                                              |                                                                                                                                      
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTcaa                                                                                                                            <  2:1176765/139‑1 (MQ=255)
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAg                                                                                                                                    >  1:571848/1‑130 (MQ=255)
tAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAg                                                                                                                                    <  2:571848/130‑1 (MQ=255)
       aaCAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCa                                                                                                                     >  1:728678/1‑139 (MQ=255)
         cAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAAc                                                                                                                   >  1:598536/1‑139 (MQ=255)
              ttACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTaa                                                                                                              >  2:678996/1‑139 (MQ=255)
                aCCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGa                                                                                                            >  2:485777/1‑139 (MQ=255)
                        tATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGaa                                                                                                    >  1:633745/1‑139 (MQ=255)
                         atcatcTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTg                                                                                                                                 <  2:532859/109‑1 (MQ=255)
                         atcatcTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTg                                                                                                                                 >  1:532859/1‑109 (MQ=255)
                           catcTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAatat                                                                                                 <  1:1101669/139‑1 (MQ=255)
                              cTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCg                                                                                              >  2:1000912/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTACAAGAGGGATGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  2:1027071/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  2:1019204/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  1:753011/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  2:389537/1‑139 (MQ=255)
                                 aaCAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGACGTGTTCAGTCTAAAAGTGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTTAACAAGTCAACATTAAGATAATGGAATATTCGgag                                                                                           >  1:751747/1‑139 (MQ=255)
                                   cAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGgagag                                                                                         <  1:741040/139‑1 (MQ=255)
                                      tttGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                                      <  2:237313/139‑1 (MQ=255)
                                      tttGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaaa                                                                                      >  2:953789/1‑139 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  2:57804/1‑105 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  1:57804/105‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAt                                                                                                                 >  2:392886/1‑110 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACAt                                                                                                                 <  1:392886/110‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    <  2:753011/139‑1 (MQ=255)
                                        tGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACg                                                                                    >  2:1061452/1‑139 (MQ=255)
                                              aaaaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  2:155633/99‑1 (MQ=255)
                                              aaaaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  1:155633/1‑99 (MQ=255)
                                              aaaaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGc                                                                              >  1:124600/1‑139 (MQ=255)
                                                aaTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAg                                                                            <  1:1061452/139‑1 (MQ=255)
                                                 aTGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGt                                                                           <  2:598536/139‑1 (MQ=255)
                                                      ggCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  1:868152/1‑90 (MQ=255)
                                                      ggCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGcc                                                                      >  1:513368/1‑139 (MQ=255)
                                                      ggCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGcc                                                                      >  1:995783/1‑139 (MQ=255)
                                                      ggCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGCGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  2:868152/90‑1 (MQ=255)
                                                         tGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaa                                                                                       <  2:499345/119‑1 (MQ=255)
                                                         tGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGaa                                                                                       >  1:499345/1‑119 (MQ=255)
                                                          gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCaa                                                                  <  1:485777/139‑1 (MQ=255)
                                                          gcgcCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAAGGTTGATGCCAGCATACAGGGTTGAAGCTGGTCAACCAGTCAACCTTAAGATAATGGAATATTCGGAGAGAAACGCCAGGCAGTTGGCCTCaa                                                                  <  2:751747/139‑1 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAATGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      <  2:930171/86‑1 (MQ=255)
                                                           cgccgcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTc                                                                                                                      >  1:930171/1‑86 (MQ=255)
                                                               gcAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             <  1:678996/139‑1 (MQ=255)
                                                                cAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTAAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  1:752152/139‑1 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         <  2:124600/139‑1 (MQ=255)
                                                                   aGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAgg                                                         <  2:633745/139‑1 (MQ=255)
                                                                    gggAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGga                                                        <  2:1116086/139‑1 (MQ=255)
                                                                     ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  2:1052341/1‑75 (MQ=255)
                                                                     ggaggaGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  1:1052341/75‑1 (MQ=255)
                                                                        ggagACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGgaagaa                                                    >  1:1015345/1‑139 (MQ=255)
                                                                          agaCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAata                                                  >  1:622169/1‑139 (MQ=255)
                                                                           gaCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatat                                                 >  1:560587/1‑139 (MQ=255)
                                                                           gaCTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAatat                                                 <  1:378095/139‑1 (MQ=255)
                                                                               gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGc                                                                                                                                   <  1:1105698/53‑1 (MQ=255)
                                                                               gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGc                                                                                                                                   >  2:1105698/1‑53 (MQ=255)
                                                                               gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTg                                                                                                                                 >  2:239139/1‑55 (MQ=255)
                                                                               gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTg                                                                                                                                 <  1:239139/55‑1 (MQ=255)
                                                                               gAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGGCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAAAATATATcc                                             >  2:375996/1‑139 (MQ=255)
                                                                                aGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATccc                                            >  2:420586/1‑139 (MQ=255)
                                                                                 gTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATg                                                               >  1:1186151/1‑119 (MQ=255)
                                                                                 gTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATg                                                               <  2:1186151/119‑1 (MQ=255)
                                                                                  tGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCa                                          >  1:1211414/1‑139 (MQ=255)
                                                                                       tgtTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGaaa                                     >  2:359056/1‑139 (MQ=255)
                                                                                        gtTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             <  2:902387/114‑1 (MQ=255)
                                                                                        gtTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGaa                                                             >  1:902387/1‑114 (MQ=255)
                                                                                          tCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAAAAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAAAGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATACCAAGGaaaaat                                  >  2:237270/1‑137 (MQ=255)
                                                                                              tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGAt                                                                                                           >  2:1146472/1‑62 (MQ=255)
                                                                                              tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGAt                                                                                                           <  1:1146472/62‑1 (MQ=255)
                                                                                              tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATg                              <  1:47218/139‑1 (MQ=255)
                                                                                              tCTAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATg                              <  1:237270/139‑1 (MQ=255)
                                                                                                tAAAAGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGtt                            <  2:1211414/139‑1 (MQ=255)
                                                                                                 aaaaGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  1:389537/139‑1 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               >  2:962126/1‑53 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTa                                                                                                               <  1:962126/53‑1 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            <  2:198821/104‑1 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAAc                                                            >  1:198821/1‑104 (MQ=255)
                                                                                                   aaGAGGGTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGgt                         >  1:1305344/1‑139 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       <  2:952312/40‑1 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGt                                                                                                                       >  1:952312/1‑40 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGa                                                                                                     <  1:994329/58‑1 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGa                                                                                                     >  2:994329/1‑58 (MQ=255)
                                                                                                        ggTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtat                    >  1:1023285/1‑139 (MQ=255)
                                                                                                         gTTGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTtata                   <  2:1015345/139‑1 (MQ=255)
                                                                                                           tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  2:622169/139‑1 (MQ=255)
                                                                                                           tGATGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATaaa                 <  2:1305344/139‑1 (MQ=255)
                                                                                                             aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:422473/127‑1 (MQ=255)
                                                                                                             aTGCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           >  1:422473/1‑127 (MQ=255)
                                                                                                               gCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           <  2:13504/125‑1 (MQ=255)
                                                                                                               gCCAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTg                           >  1:13504/1‑125 (MQ=255)
                                                                                                                ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac            <  2:560587/139‑1 (MQ=255)
                                                                                                                ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac            <  2:1023285/139‑1 (MQ=255)
                                                                                                                ccAGCATACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAacacac            <  1:1027071/139‑1 (MQ=255)
                                                                                                                     aTACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACa       <  1:953789/139‑1 (MQ=255)
                                                                                                                     aTACAGGGGTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACa       <  2:995783/139‑1 (MQ=255)
                                                                                                                          ggggTGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAt  <  1:359056/139‑1 (MQ=255)
                                                                                                                              |                                                                                                                                      
TAGTACCAACAGATTTACCACATTTATCATCTTAACAGTTTGAGTGAAAATGCAGGCTGCGCCGCAAAGGGAGGAGACTGAGTGGCGTGTTCAGTCTAAAAGAGGGTTGATGCCAGCATACAGGGGCGAAGCTGGTCAACAAGTCAACATTAAGATAATGGAATATTCGGAGAGAAACGTCAGGCAGTTGGCCTCAAATGAACAGGAAGAATATATCCCCAGGAAAATTAATGTTGGTGTTATAAACACACCAACATTAAT  >  NC_000913/312063‑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: