Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,375,431 A→C V76G (GTA→GGA)  epmB ← EF‑P‑Lys34 lysylation protein; weak lysine 2,3‑aminomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,375,4310AC97.4% 293.3 / NA 76V76G (GTA→GGA) epmBEF‑P‑Lys34 lysylation protein; weak lysine 2,3‑aminomutase
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (38/36);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTACCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTT  >  NC_000913/4375283‑4375575
                                                                                                                                                    |                                                                                                                                                
gCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTACGTccc                                                                                                                                                >  1:191516/1‑151 (MQ=255)
     aaTTTACCGCGCCGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCCGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTACCTgcc                                                                                                                                            >  1:257418/1‑148 (MQ=255)
      aTTTACCGCGCATCCACCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGTGCGCGATGACAAACTCGTCTTGCGAGGTAAGTCCCTCAc                                                                                                                                            <  1:626334/149‑1 (MQ=255)
        ttACCGCGCAGCCACCTGTGACCAGCAAAATCGCCCGGTTGTGGTATTTATGCAACAAACGAGGCACTATGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCGTCTTGCGAGGTAAGTCCCTGACGCa                                                                                                                                         >  1:401842/1‑150 (MQ=255)
          aCCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCaaa                                                                                                                                       <  1:786612/150‑1 (MQ=255)
             gcgcAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAagag                                                                                                                                   <  1:901355/151‑1 (MQ=255)
                cAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGAt                                                                                                                                <  2:1052395/151‑1 (MQ=255)
                     gccTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTcc                                                                                                                           >  2:260399/1‑151 (MQ=255)
                      ccTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCg                                                                                                                          >  2:194416/1‑151 (MQ=255)
                        tttGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCgg                                                                                                                         >  2:477066/1‑150 (MQ=255)
                           gACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGAtt                                                                                                                      <  1:1071380/150‑1 (MQ=255)
                           gACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTg                                                                                                                     <  2:170988/151‑1 (MQ=255)
                            acaaGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCCTCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGAtt                                                                                                                      >  1:348780/4‑149 (MQ=255)
                            aCCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGc                                                                                                                    >  2:107486/1‑151 (MQ=255)
                            aCCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGc                                                                                                                    <  2:141529/151‑1 (MQ=255)
                               aGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCtt                                                                                                                 <  1:826320/151‑1 (MQ=255)
                                gCAAAAGCGCCCGGTTGTGGTCTTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCttt                                                                                                                >  1:995792/1‑151 (MQ=255)
                                   aaaGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTtctc                                                                                                             >  2:888460/1‑151 (MQ=255)
                                      gcgcCCGGTTTTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCa                                                                                                           <  1:64952/150‑1 (MQ=255)
                                      gcgcCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCAt                                                                                                          <  1:1047645/151‑1 (MQ=255)
                                        gcCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCACAAGAGGATCGTCCGGATTGCCTTTCTCCATGc                                                                                                        >  1:800987/1‑151 (MQ=255)
                                        gcCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGc                                                                                                        >  1:384762/1‑151 (MQ=255)
                                        gcCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGc                                                                                                        <  2:795281/151‑1 (MQ=255)
                                           gggTTGTGGTATTTATGTAACAAAAAAGGAAATACGATGTGATGTTATTGCAGTGGGTAGGTGGAGAATACGGGCGCGATGACAAATTCATCTTGCGAGGTAAGTACCTGACGCACACGAGGATCGTCAGGATTGCCTTTATTCATGcgat                                                                                                     >  1:359090/2‑151 (MQ=255)
                                            ggTTGTGGTATTTATGCAACAAAACAGGCCCTAAGCTGAGCTGCTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGGCAAACTCATCTTGTGAGGTAAGTCCCTGACGCAAAAGAGGACCGTCCGGATTGCCTTTCTCCATGc                                                                                                        >  1:261797/1‑147 (MQ=255)
                                               tgtgGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGCGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGCCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGcgatcg                                                                                                   >  1:917686/1‑149 (MQ=255)
                                               tgtgGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGcgatcga                                                                                                  <  2:924907/150‑1 (MQ=255)
                                                   gTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAAt                                                                                             >  1:546757/1‑151 (MQ=255)
                                                   gTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAAt                                                                                             <  2:720962/151‑1 (MQ=255)
                                                         aTGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCGGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAg                                                                                          <  1:204803/148‑1 (MQ=255)
                                                           gCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCgggg                                                                                     <  2:387335/151‑1 (MQ=255)
                                                           gCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCggg                                                                                      <  2:59506/150‑1 (MQ=255)
                                                                   ccAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGc                                                                              >  1:708115/1‑150 (MQ=255)
                                                                    cAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAg                                                                            <  2:1030315/151‑1 (MQ=255)
                                                                        cACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCa                                                                        <  2:871284/151‑1 (MQ=255)
                                                                         aCTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCaa                                                                       <  1:666150/151‑1 (MQ=255)
                                                                             cGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCaaaa                                                                     >  1:922537/1‑149 (MQ=255)
                                                                              gctgtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGc                                                                  >  1:127626/1‑151 (MQ=255)
                                                                              gctgtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGc                                                                  >  2:922989/1‑151 (MQ=255)
                                                                              gctgtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGc                                                                  >  2:502633/1‑151 (MQ=255)
                                                                              gctgtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGc                                                                  >  2:613312/1‑151 (MQ=255)
                                                                                tgtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCtt                                                                >  2:731715/1‑151 (MQ=255)
                                                                                 gtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCttt                                                               <  1:570695/151‑1 (MQ=255)
                                                                                 gtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCttt                                                               <  2:431464/151‑1 (MQ=255)
                                                                                 gtgctgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCtt                                                                <  2:1073902/150‑1 (MQ=255)
                                                                                     tgtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTgg                                                           >  2:10591/1‑151 (MQ=255)
                                                                                      gtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTgg                                                           <  2:302009/150‑1 (MQ=255)
                                                                                      gtTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGgc                                                          <  2:378893/151‑1 (MQ=255)
                                                                                        tcttcCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGgcgc                                                        >  1:271154/1‑151 (MQ=255)
                                                                                         cttcCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGgcg                                                         <  1:197853/149‑1 (MQ=255)
                                                                                         cttcCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCt                                                       >  2:1011342/1‑151 (MQ=255)
                                                                                              agagggTCGGTGGAGAATCCGGGGGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCAAACGCCGGGCAAAAAGCTTTTTGGCGCTgag                                                    >  1:22862/4‑147 (MQ=255)
                                                                                                  ggTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCgg                                               >  1:585607/1‑150 (MQ=255)
                                                                                                     cGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTa                                            >  2:1066288/1‑150 (MQ=255)
                                                                                                     cGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGCAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGAGCTGCGTCCGGCTa                                            >  2:434865/1‑150 (MQ=255)
                                                                                                          gagaATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGtt                                      <  1:584420/151‑1 (MQ=255)
                                                                                                            gaATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGtttt                                    >  2:952903/1‑151 (MQ=255)
                                                                                                                   gccgcgATGAAAAACTCCTCTTGCGAGGTAAGTCCCTGACGAAAAAGAGGATCGTCCGGATTGCCTCTCTCCATGCGATCGATAAATGCGCGGGGCACACGCAGGGCAAAAAGCCTTTTGGCGCTGCGCTCGGCCAACAGCTTTttcttc                              >  1:474680/3‑148 (MQ=255)
                                                                                                                    gcgcGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCcgc                            <  2:1090768/151‑1 (MQ=255)
                                                                                                                    gcgcGATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCcgc                            <  2:552800/151‑1 (MQ=255)
                                                                                                                       cgATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTc                         <  1:214587/151‑1 (MQ=255)
                                                                                                                        gATGACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCt                        >  2:241025/1‑151 (MQ=255)
                                                                                                                           gACAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCtat                      >  1:173340/1‑150 (MQ=255)
                                                                                                                           gACAAACGCCTCTTGCGAGGGAAGGCCCTGACGCAAAAGAGGATCGTCCGGAGTGCCTTTCTCCATGCGGTCGGTAAATGAGCGGGGCACAAGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCTGGCTAACAGTTTTTCCTCcgcg                           <  1:793446/145‑1 (MQ=255)
                                                                                                                            aCAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTatat                    <  2:346887/151‑1 (MQ=255)
                                                                                                                             cAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATAtt                   >  1:412587/1‑151 (MQ=255)
                                                                                                                             cAAACTCATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATAtt                   >  1:568740/1‑151 (MQ=255)
                                                                                                                                  tCATATTGCGAGGTAAGCCCATGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAACGGGGCACACGCCGCGCAAAAAGCCTTTTCGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTCTATTCac                >  1:18914/1‑148 (MQ=255)
                                                                                                                                   cATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCaaaaa             <  1:260295/151‑2 (MQ=255)
                                                                                                                                   cATCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCaaaa              <  1:929469/150‑1 (MQ=255)
                                                                                                                                   cATCTTGCGAGGTAAATCCCTGACGCAAAAGAGGATCGTCCGGAATGCCTTTCTCCATGCGATCGATAAATGAGCGGTGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCTGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAg             <  1:31054/151‑1 (MQ=255)
                                                                                                                                    aTCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCaaa               >  2:633902/1‑148 (MQ=255)
                                                                                                                                    aTCTTGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGATAACAGTTTTTCCTCCGCGTCTATATTcta                >  1:652645/1‑145 (MQ=255)
                                                                                                                                        tGCGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGa            <  1:628846/147‑1 (MQ=255)
                                                                                                                                          cGAGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGc         >  2:368758/1‑148 (MQ=255)
                                                                                                                                            aGGTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAg    <  2:504414/151‑1 (MQ=255)
                                                                                                                                              gTAAGTCCCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGtt  >  2:425889/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCGGCAATTTACCGCGCAGCCGCCTTTGACCAGCAAAAGCGCCCGGTTGTGGTATTTATGCAACAAACCAGGCACTACGCTGTGCTGTTCTTCCAGTGGGTCGGTGGAGAATCCGGGCGCGATGACAAACTCATCTTGCGAGGTAAGTACCTGACGCAAAAGAGGATCGTCCGGATTGCCTTTCTCCATGCGATCGATAAATGAGCGGGGCACACGCAGGGCAAAAAGCTTTTTGGCGCTGCGTCCGGCTAACAGTTTTTCCTCCGCGTCTATATTCAAAAGACGCAGAAGTT  >  NC_000913/4375283‑4375575

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