Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,505,884 A→G L46P (CTG→CCG)  yhfS ← FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,505,8840AG99.1% 350.2 / ‑3.8 107L46P (CTG→CCG) yhfSFNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (54/51);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 107 total aligned reads displayed.

CCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTGCCTTC  >  NC_000913/3505738‑3506030
                                                                                                                                                   |                                                                                                                                                  
ccAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTc                                                                                                                                                 <  2:570678/150‑1 (MQ=255)
    aagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGCAACGCGCGGTTGATTCGGTCCTGGc                                                                                                                                            >  2:368937/1‑151 (MQ=255)
      gaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGt                                                                                                                                          <  2:383774/151‑1 (MQ=255)
       aCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTTGGTCCTGGCGtt                                                                                                                                         <  1:379289/151‑1 (MQ=255)
       aCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGtt                                                                                                                                         >  2:1205596/1‑151 (MQ=255)
       aCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGtt                                                                                                                                         <  2:1261241/151‑1 (MQ=255)
       aCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGtt                                                                                                                                         <  1:203937/151‑1 (MQ=255)
            gCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAg                                                                                                                                      <  1:439557/149‑1 (MQ=255)
            gCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGcc                                                                                                                                    >  1:313867/1‑151 (MQ=255)
             ccccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCt                                                                                                                                   <  3:277300/151‑1 (MQ=255)
                cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTaaa                                                                                                                                <  1:238898/151‑1 (MQ=255)
                cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTAAAGCCTaaa                                                                                                                                <  1:433407/151‑1 (MQ=255)
                 ggTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAAt                                                                                                                               <  2:13952/151‑1 (MQ=255)
                     tGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCGGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAATTCAcc                                                                                                                           >  3:177103/1‑151 (MQ=255)
                      gAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCg                                                                                                                          <  2:1241574/151‑1 (MQ=255)
                          aaaGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAg                                                                                                                      <  1:66947/151‑1 (MQ=255)
                            aGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGt                                                                                                                     >  2:248539/1‑150 (MQ=255)
                              ccgccAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTa                                                                                                                    >  4:339084/1‑149 (MQ=255)
                               cgccAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTa                                                                                                                    >  2:1020760/1‑148 (MQ=255)
                               cgccAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTa                                                                                                                    >  3:536281/1‑148 (MQ=255)
                                  cAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGaaa                                                                                                               >  2:942796/1‑150 (MQ=255)
                                   aGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACt                                                                                                             >  2:644249/1‑151 (MQ=255)
                                     cccGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTCCCACCCGCTGGGTAACGCGCGGTTGATTCGGTCTGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCg                                                                                                           <  2:968143/151‑1 (MQ=255)
                                     cccGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCg                                                                                                           >  2:878659/1‑151 (MQ=255)
                                        ggCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTg                                                                                                        >  2:1343146/1‑151 (MQ=255)
                                             gAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCggg                                                                                                   >  2:897820/1‑151 (MQ=255)
                                              aaccgcgCCAGCCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTCCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCg                                                                                                     <  2:125632/145‑1 (MQ=255)
                                              aaccgcccctccccccGCGCCTTGCACCAGCGCCGCAGCTTGTTCGTAAATTCCATCACCTAACACCTGTCCCACACGCGGGTTAACGCGCGGTTGTTTCGGTCTTGGCTTTAACCGTAAATCACCGCCAGTAAGAAACTCGCTGCCgg                                                                                                    <  2:182452/137‑1 (MQ=255)
                                              aaTCGCGCCAGTCCCCGCCCCCTGCCCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTCCCACACGCTGGGTACCGCGGGGTTGATTCGGTCCTGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCggg                                                                                                   <  2:302382/150‑1 (MQ=255)
                                               aTCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACATGTTCCACACGCTGCGTAACGCGCGGTTGATTCGGTCCTGACGTTAAACCTAAATCACCGCCAGTAAGAAACTCGCTGCCggg                                                                                                   >  3:848401/1‑149 (MQ=255)
                                               aTCGCGCAAGCCCCCGCGCCCTGCCCCAGCGCCGCAGCCTGTGCGTAAAATGCATCAGCCAGCACCTGTTCCCCACGCGGGGTAACGCGCGTTGGATTCGGTCCTGGGGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGa                                                                                                 <  1:19881/151‑1 (MQ=255)
                                                 cgcgCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaa                                                                                                >  3:410826/1‑150 (MQ=255)
                                                 cgcgCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaa                                                                                                >  2:1380238/1‑150 (MQ=255)
                                                     ccAGCCCCCGCCCCCTGCCCCAGCGCCGCGGCCTGGGCGTGAAATGCATCAGCCAGCACCTGTTCCACCCGCTGGTTAACGCGGGGTTGATTGGGTCCGGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaaa                                                                                               <  1:220379/147‑1 (MQ=255)
                                                         cccccGCCCCCGGCCCCGGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCCCTTTTTCCCCACGCTGGGTAACGCGCGGTTGATTGGTTCCTGGCGTAAAGCTAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACa                                                                                       <  3:615825/150‑1 (MQ=255)
                                                          ccccGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACa                                                                                       >  2:708097/1‑150 (MQ=255)
                                                             cgcgCCCTGCCCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATg                                                                                   <  1:59212/151‑1 (MQ=255)
                                                                gcCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCt                                                                                 >  1:282406/1‑150 (MQ=255)
                                                                gcCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCt                                                                                 >  2:423037/1‑150 (MQ=255)
                                                                 ccctcaaCCAGCCCCGAAGCCTGTGCGTGAAGTCCATATGCCACCCCCTGTCCCACTGGCTGGGTAGCGCGCGTTTGTTTCGTCCCTGGGGATAAGCCAAAATCGCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTAt                                                                               <  2:410796/145‑1 (MQ=255)
                                                                 cccTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTAt                                                                               >  2:911401/1‑151 (MQ=255)
                                                                     gCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccac                                                                           >  3:1074524/1‑151 (MQ=255)
                                                                       aCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacc                                                                          >  2:826140/1‑150 (MQ=255)
                                                                       aCCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacc                                                                          >  2:770906/1‑150 (MQ=255)
                                                                        ccAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacc                                                                          >  1:355945/1‑149 (MQ=255)
                                                                          aGCCCCGCGGCCTGTGCGTGAAATCCTTCGCCCAGCCCCTTTCCCACACGCGGGTAACCGCGCGTTTGTTTGGGTCCGGGCGTAAAGCAAAATTCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATccacca                                                                         <  3:1132767/148‑1 (MQ=255)
                                                                            cgccgcAGCTTGTGCGTGAAATGCCTCAGCCAGCCCCTGTTCCACACGCTGGGTACCGCGCGGTTGATTCGGTCTTGGCGTAAAGCCTAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgcgc                                                                     <  4:368670/150‑1 (MQ=255)
                                                                              ccgcGGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTTCCACACGCGGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgcg                                                                      <  2:824454/147‑1 (MQ=255)
                                                                              ccgcAGCCTGCGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCaa                                                                   >  3:1011612/1‑150 (MQ=255)
                                                                               cgcAGCCTGTGCGTGAAATGCATCAGCCAGCCCCTGTCCCCCACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAACCCTAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                                                                 <  2:382040/151‑1 (MQ=255)
                                                                               cgcAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                                                                 >  3:758370/1‑151 (MQ=255)
                                                                               cgcAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                                                                 >  2:760899/1‑151 (MQ=255)
                                                                               ccccgccggggCGTGAAGTGCCCCGCCCAGCTCTTTTTCCACCGGATGGCCAACGCGCGGTTTATTCGGTCCTGGGGTGACGCCAAAACCTCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCaa                                                                   <  3:868841/139‑1 (MQ=255)
                                                                                 cAGCCTGGGCGTGAAATGCCCCACCCAGCCCCTGTCCCCCCCGCGGGGTACCGCGCGGTTGATCGGGTCCGGGCGTTAAGCCAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                                                                 <  3:418331/149‑1 (MQ=255)
                                                                                   gCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCt                                                             >  1:188528/1‑151 (MQ=255)
                                                                                    cgcgcgggtaataTGCATCAGCAACCACCTGTTCCACACGCTGGGTAACGGGCGGGTTGTTGTGTCCCGGGGTTAAGCCTAAAACGCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAAc                                                                 <  4:265401/134‑1 (MQ=255)
                                                                                    ccTGGGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCt                                                             <  3:533416/150‑1 (MQ=255)
                                                                                     cTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCttt                                                           >  3:285946/1‑151 (MQ=255)
                                                                                       gtgCGTAAAATGCATCAGCCAGCCCCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCttt                                                           <  1:133701/149‑1 (MQ=255)
                                                                                       gtgCGTAAAATCAATCACCAAGCACCTGTACCACCCGCTGGGTAACGCGCGGTTGATTCGGCCCGGGCGTAACGCTAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtttt                                                          <  3:124288/150‑1 (MQ=255)
                                                                                        tgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgc                                                        >  2:724720/1‑151 (MQ=255)
                                                                                        tgCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgc                                                        >  2:587694/1‑151 (MQ=255)
                                                                                         gCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgct                                                       >  3:707784/1‑151 (MQ=255)
                                                                                         gCGTAAAATGCATCAGCCAGCCCCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgct                                                       <  2:536247/151‑1 (MQ=255)
                                                                                          cGTGAAATGCATCAGCCAGCCCCTGTCCCACACGCTGGGAAACGCGCGGTGGATTGGGTCCTGGCTTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtttt                                                          <  3:354427/147‑1 (MQ=255)
                                                                                           gTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgctgc                                                     >  1:247970/1‑151 (MQ=255)
                                                                                              aaaTGCATCAGCCAGCACCTTTCCCACCCGCGGGTTAACGCGGGTTTGATCGGGTCCGGGCTTAAAGCTAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgct                                                       <  3:1128902/146‑1 (MQ=255)
                                                                                              aaaTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGcc                                                  >  3:29423/1‑151 (MQ=255)
                                                                                                aTGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTc                                                >  3:1088726/1‑151 (MQ=255)
                                                                                                   cATCAGCCAGCCCCTGTCCCCCACGCTGGGTACCGCGGGGTTGATGCGGTCCGGGCGTAAACCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCa                                               <  3:248564/149‑1 (MQ=255)
                                                                                                   cATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaat                                             >  3:624876/1‑151 (MQ=255)
                                                                                                    aTCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaat                                             >  2:494227/1‑150 (MQ=255)
                                                                                                     tCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaata                                            >  1:306367/1‑150 (MQ=255)
                                                                                                       agccagcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATc                                         >  3:1234580/1‑151 (MQ=255)
                                                                                                            gcACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTc                                      >  3:1210494/1‑149 (MQ=255)
                                                                                                               ccTGTCCCCCCCGCTGGGTACGGCGCGGTTGATTCGGTCCTGGGGTGAAGCATAAACCCCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGc                                  <  3:651010/150‑1 (MQ=255)
                                                                                                                   ttCCACCCGCGGGGTACCGCGCGGTTGATTCGGTCTTGGCGTAAAGCAAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCtt                                <  3:691344/148‑1 (MQ=255)
                                                                                                                   ttCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGc                             >  1:182287/1‑151 (MQ=255)
                                                                                                                    tCCACACGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCACCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGc                             >  2:1240346/1‑150 (MQ=255)
                                                                                                                     ccccacccggggAAACGCGGGTTTTATTCGGGCCTGGGGTTAACCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCa                            <  3:368969/141‑1 (MQ=255)
                                                                                                                       acacGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCag                           >  2:324958/1‑149 (MQ=255)
                                                                                                                        cacGCTGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAgg                        >  1:111945/1‑151 (MQ=255)
                                                                                                                         acGCTGGGTACCGGGGGTTTGTTCCGGTCCTGGGGTTAAGCTAAAAACACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCaga                          <  2:544992/148‑1 (MQ=255)
                                                                                                                           gCTGGGAAACGCGCGGTTTATTGGGCCCGGGCGTTAAGCCTAAACCCCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCagag                         <  2:177077/147‑1 (MQ=255)
                                                                                                                             tGGGTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACg                   >  1:196616/1‑151 (MQ=255)
                                                                                                                              gggTAACGCGGGGTTGATTCGGTCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGaaa                     <  2:49419/148‑1 (MQ=255)
                                                                                                                              gggTAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGt                  >  2:819947/1‑151 (MQ=255)
                                                                                                                                 tAACGCGCGGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCtt               >  4:28536/1‑151 (MQ=255)
                                                                                                                                   accgcgggTTTGATTCGGTCCTGGCTTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCt                <  1:85178/147‑1 (MQ=255)
                                                                                                                                   accgcgcgGTTGATTCGGTCCTGGCGTTAAGCCTAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCtt               <  3:43765/148‑1 (MQ=255)
                                                                                                                                        ggggTTGATTGGGTCCGGGCGTTAAGCCAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCAt            <  3:1371816/146‑1 (MQ=255)
                                                                                                                                        gcgGTTGATTCGGTCCTGGCGTTAACCCTAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAgg         <  2:1234521/151‑1 (MQ=255)
                                                                                                                                        gcgGTTGATTCGGTCCTGGCGTAAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAg          <  2:1011569/150‑1 (MQ=255)
                                                                                                                                         cgGTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGt        >  2:456488/1‑151 (MQ=255)
                                                                                                                                          ggTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTg       >  3:508091/1‑151 (MQ=255)
                                                                                                                                          ggTTGATTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTg       >  2:1364259/1‑151 (MQ=255)
                                                                                                                                          ggTTGATTCGGTCCGGGCGTAAAGCCAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCa             <  3:897856/145‑1 (MQ=255)
                                                                                                                                          ggTTGATCCGGTCGTGGCGTTAAGCCTAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAgg         <  3:1205657/149‑1 (MQ=255)
                                                                                                                                              gATTCGGTCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTGcct    <  1:257472/150‑1 (MQ=255)
                                                                                                                                               aTTCGGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTGccttc  >  4:281124/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGCAGAGGAAACGTCTTCATAGGTGCCTTC  >  NC_000913/3505738‑3506030

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