Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,505,857 A→G V55A (GTG→GCG)  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,8570AG96.2% 226.1 / ‑3.0 82V55A (GTG→GCG) yhfSFNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (35/41);  total (36/46)
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.91e-01

GTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGC  >  NC_000913/3505710‑3506003
                                                                                                                                                   |                                                                                                                                                   
gtcgtcGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCgcac                                                                                                                                                  >  2:435573/1‑151 (MQ=255)
  cgtcgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGc                                                                                                                                                >  1:253264/1‑151 (MQ=255)
  cgtcgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGc                                                                                                                                                >  2:1363727/1‑151 (MQ=255)
     cgGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTgg                                                                                                                                             <  2:1358741/151‑1 (MQ=255)
       ggTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGt                                                                                                                                           <  1:1331444/151‑1 (MQ=255)
          aaaCAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCAGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTaa                                                                                                                                        <  2:1539595/150‑1 (MQ=255)
            aCAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAAcg                                                                                                                                      >  2:409916/1‑150 (MQ=255)
              aGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAAcgcgc                                                                                                                                   <  2:68105/151‑1 (MQ=255)
                gcgcGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCgg                                                                                                                                 <  2:438980/151‑1 (MQ=255)
                gcgcGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCgg                                                                                                                                 <  2:1563645/151‑1 (MQ=255)
                 cgcgTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGt                                                                                                                                >  2:9972/1‑151 (MQ=255)
                  gcgTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGtt                                                                                                                               <  2:336049/151‑1 (MQ=255)
                     tCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGAt                                                                                                                            >  2:40681/1‑151 (MQ=255)
                         gCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAg                                                                                                                        <  1:556107/151‑1 (MQ=255)
                             cAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCt                                                                                                                    <  2:224227/151‑1 (MQ=255)
                              agaagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTg                                                                                                                   <  1:1159826/151‑1 (MQ=255)
                                aagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGc                                                                                                                 >  1:1311843/1‑151 (MQ=255)
                                aagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGc                                                                                                                 >  2:267441/1‑151 (MQ=255)
                                aagaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGc                                                                                                                 <  1:1181853/151‑1 (MQ=255)
                                 agaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCg                                                                                                                >  2:1282571/1‑151 (MQ=255)
                                 agaCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCg                                                                                                                >  2:118886/1‑151 (MQ=255)
                                   aCGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGtt                                                                                                              >  1:1510945/1‑151 (MQ=255)
                                       tGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTaa                                                                                                            <  1:1071576/149‑1 (MQ=255)
                                       tGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGc                                                                                                          >  1:969178/1‑151 (MQ=255)
                                          cccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTa                                                                                                       <  2:420610/151‑1 (MQ=255)
                                           ccGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTa                                                                                                       >  1:183547/1‑150 (MQ=255)
                                            cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTaaa                                                                                                     <  1:394853/151‑1 (MQ=255)
                                            cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTaaa                                                                                                     <  1:587868/151‑1 (MQ=255)
                                            cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTaaa                                                                                                     <  2:80651/151‑1 (MQ=255)
                                            cGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTaaa                                                                                                     <  1:759595/151‑1 (MQ=255)
                                              gTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAAt                                                                                                    >  1:1504580/1‑150 (MQ=255)
                                                  gAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCAcc                                                                                                <  2:968434/150‑1 (MQ=255)
                                                     cAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGcc                                                                                             >  1:1305771/1‑150 (MQ=255)
                                                      aaaGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAg                                                                                           <  2:1620703/151‑1 (MQ=255)
                                                              cAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAAc                                                                                   >  1:133836/1‑151 (MQ=255)
                                                               aGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACt                                                                                  >  2:80624/1‑151 (MQ=255)
                                                                 ccgggCGCGAATCGCGCCAGCCCCCGCCCCCTGCCCCAGCGCCGCGGCTTGTGCGTAAAATGCATCAGCCAGCCCCTGTCCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTc                                                                                 <  1:1676170/147‑1 (MQ=255)
                                                                   cGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCt                                                                              >  2:1134997/1‑151 (MQ=255)
                                                                    ggCGCGACCCGCGCCAGCCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTTTTCCGCCCGCTGGGTACCGCGCGGTTGTTTCAGTCCTGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCt                                                                              <  2:247731/150‑1 (MQ=255)
                                                                    ggCGCGAACCGCGCCGGTCCCCGCGCCCTCCCCCAGCGCCGCGGCCGGGGCGTAAAATGCATCAGCCAGCACCTGTTCCGCACGTTGGGTAAGGCGGGGTTGTTCCAGCCCGGGCGTTAAGCCTAAACCACCGCCAGTAAGAAACTCGCt                                                                              <  2:1483511/150‑1 (MQ=255)
                                                                       gggAATCGCCCCAGCCCCCGCGCCCTGCCCCGGCCCCGCGGCCTGTGGGTGAAATGCATCACCCACCCCCTGTCCCCCACGCTGGGAACCGGGCGGTTGATTCAGCCCTGGCGTTAAGCCAAAACCACCGCCAGTAAGAAACTCGCTGcc                                                                           <  1:526184/148‑1 (MQ=255)
                                                                           accgcgCCAGCCCCCGCCCCCTGCCCCGGCGCCGCAGCTTGTGCGTGAAATGAACCAGCCAGCCCCTGTCCCGCCCGCTGGGTACCGCGCGGTTGATTCAGTCCTGGCTTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCggg                                                                        <  1:737107/147‑1 (MQ=255)
                                                                               cgcCAGCCCCCGCGCCCCGCACCGGCGCCGCGGCCTGTGCGTAAAATGCATCACCCAGCACCTGTTCCGCCCGCGGGGTAACGCGCGGTTGATACAGCCCTGGCGTAAAGCCAAAAACACCGCCAGTAAGAAACTCGCTGCCGGGGAAATg                                                                  <  2:1450024/151‑1 (MQ=255)
                                                                               cgcCAGCCCCCGCCCCCTCCCCCAGCGCCGCAGCCTGTGCGTAAAATGCATCAGCCAGCACCTGTCCCGCACGCTGGGTACCGCGGGGTTGATTCAGTCCGGGCGTTAAGCAAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGaa                                                                     <  2:510510/148‑1 (MQ=255)
                                                                                 cgattCCCCGCCCCCAGCCCCTGGGCCGCACGCTGGGGGTGAAATAAAAAAGCCGGCCCGTGTCCCGCCCGGTGGGTAACGGGGGGTGAATTCAGCCCGGGCTTGATACCTAATTCACCGCCAGTAAGAAACTCGCTGCCgggg                                                                       <  1:98175/140‑1 (MQ=255)
                                                                                 ccAGCCCCCGCCCCCTGCCCCGGCGCCGCAGCTTGGGCGTGAAAGGCACCAGCCAGCACCTGTTCCGCACGGTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTAAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATgg                                                                 <  1:45317/150‑1 (MQ=255)
                                                                                       cccGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaa                                                           >  1:1224536/1‑150 (MQ=255)
                                                                                       cccGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACaa                                                           >  1:101463/1‑150 (MQ=255)
                                                                                        ccGCGCCCGGCCCCGGCGCCGCGGCCTGTGCGTGAAATGCAACAGCCAGCCCCTGTCCCACACGCTGGGTAACGCGGGGTTGATTCAGCCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAAt                                                         <  1:1242343/151‑1 (MQ=255)
                                                                                         cgcgCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATg                                                        >  1:733135/1‑151 (MQ=255)
                                                                                           cgcCCTGCCCCAGCGCCGCAGCTTGTGCGTAAAATGCATCAGCCAGCCCCTGTCCCGCCCGCGGGGTAACGCGCGGTGGATTCAGTCCGGGCGTTAAGCCTAAATCTCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCt                                                      <  1:1557569/151‑1 (MQ=255)
                                                                                                tGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATcc                                                  >  2:1662427/1‑150 (MQ=255)
                                                                                                     cgccacggcggCCTGGGGGTGAAACGCATCCGCCACCCCTGGTACCGCACGCTGGGTAACGGGGGGTGGATTCGGCCCTGGCTTTAAGCCTAAAACACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgc                                            <  2:1491721/141‑1 (MQ=255)
                                                                                                        cgccgcAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgcgc                                          >  2:490126/1‑150 (MQ=255)
                                                                                                          ccgcGGCCGGTGGGTAAAATGCATCCGCCACCCCCTGTCCCGCCCGCTGGGTAACGCGGGGTTGATTCAGTCCGGGCGTAAACCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAgc                                            <  2:402525/146‑1 (MQ=255)
                                                                                                          ccgcAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCaa                                        >  1:1707889/1‑150 (MQ=255)
                                                                                                             cAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTg                                    >  2:450872/1‑151 (MQ=255)
                                                                                                                ccTGTGCGTAAAATGCATCAGCCAGCACCTGTCCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGc                                   <  2:904464/149‑1 (MQ=255)
                                                                                                                  tgggCGTGAATAGCTTCCGCCAGCACCTGTGCCGCACGTGGGTAAACGCGGGGTTTTTTCAGTCTTGGGGTAAAGCCAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTgtt                                  <  1:257449/145‑3 (MQ=255)
                                                                                                                   gtgCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtg                              >  2:1461928/1‑151 (MQ=255)
                                                                                                                   gggCGTGAAATGCACCAGCCAGCCCGTGTCCCGCACGCTGGGTACCGCGGGGTTGGTTCAGTCCGGGCGTTAAGCTTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCtt                                 <  1:1293476/146‑1 (MQ=255)
                                                                                                                   gggCGTAAAATGCACCGGCCAGCCCCTGTCCCGCACGCGGGGTACCGCGGGGTTGTTTCAGCCCGGGCGTAAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACt                                     <  1:1675119/142‑1 (MQ=255)
                                                                                                                    tgCGTGAAATGCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgc                             >  2:59842/1‑151 (MQ=255)
                                                                                                                          aaaTGCTTCAGCAGCCTCCTGTCCCGCCCGCGGGGTACCGGGCGGGTGATTCGGTCCTGGCGTAAAACCAAAACCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgctg                           <  1:1151121/147‑1 (MQ=255)
                                                                                                                          aaaTGCATCAGCCAGCACCTGTCCCACACGGGGGGTAACGCGGGGTTGATTCAGTCCGGGCGTAAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTtgctgc                          <  1:118878/148‑1 (MQ=255)
                                                                                                                              gCATCAGCCAGCACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCa                    >  1:316882/1‑150 (MQ=255)
                                                                                                                               caccagccggcCCCTTTCCCGCACGCTGGGTATCGCGCGGTTATTTCAGCCCGGGCGTAAGGCCAAAAACACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaat                  <  1:1363635/148‑1 (MQ=255)
                                                                                                                               cATCAGCAACCCCCTGTTCCGCACGCGGGGTAACGCGCGGTTGATTCAGTCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaa                   <  2:1224631/150‑1 (MQ=255)
                                                                                                                               cATCACCCAGCGCCTGTCCCCCACCCTGGGTAACGCGCGGTTTATACAGCCCTGGCGTTAAGCCCAAATCCCCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCa                    <  2:1716749/149‑1 (MQ=255)
                                                                                                                                   agccagcACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATc              >  1:615683/1‑151 (MQ=255)
                                                                                                                                     ccagcACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGt            >  1:1349004/1‑151 (MQ=255)
                                                                                                                                     ccagcACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGt            >  1:1202403/1‑151 (MQ=255)
                                                                                                                                     ccacccccTGTTCCTCCCGCTGGCTACCGCGCGTTTGATTCAGTCCGGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCaataa                <  1:424523/141‑1 (MQ=255)
                                                                                                                                      cagcACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTc           >  2:1380375/1‑151 (MQ=255)
                                                                                                                                         cACCTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAgg        >  2:1373579/1‑151 (MQ=255)
                                                                                                                                         cACCTGTCCCGCCCGTGGGGTACCGCGCGGTTGATTCAGCCCGGGCGTAAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAgg        <  2:746422/151‑1 (MQ=255)
                                                                                                                                          aCCTGTCCCGCACGCGGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCAAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAg         <  1:367675/149‑1 (MQ=255)
                                                                                                                                           ccTGTCCCGCACGCTGGGTAACGCGGGGTTGATTCAGCCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTc           <  2:521609/146‑1 (MQ=255)
                                                                                                                                           ccTGTCCCGCACGCTGGGTAACGCGGGGTTGATCCAGCCCGGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTc           <  1:9969/146‑1 (MQ=255)
                                                                                                                                            cTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCtt     >  2:462711/1‑151 (MQ=255)
                                                                                                                                            cTGTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCtt     >  2:1533888/1‑151 (MQ=255)
                                                                                                                                              gTTCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGAGACAAATGCTATCCACCAGCGCAAACTGCTTTTCTGCGCCTCAATAATCGTCAGGCTTTg   >  1:864717/1‑150 (MQ=255)
                                                                                                                                               ttCCGCACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGc  >  2:1497969/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTCGTCGGGTAAACAGGCGCGTCATGCACCAGAAGACGCTGCCCCGGTTTGAGCAAAGCCGCCAGCCCGGCGCGAATCGCGCCAGTCCCCGCGCCCTGCACCAGCGCCGCAGCCTGTGCGTGAAATGCATCAGCCAGCACCTGTTCCACACGCTGGGTAACGCGCGGTTGATTCAGTCCTGGCGTTAAGCCTAAATCACCGCCAGTAAGAAACTCGCTGCCGGGGAAATGGCGACAAATGCTATCCACCAGCGCAAACTGCTTTTGCTGCGCCTCAATAATCGTCAGGCTTTGC  >  NC_000913/3505710‑3506003

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