Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,425,968 T→G 8.4% intergenic (‑52/+38) hisQ ← / ← hisJ lysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisQ/histidine ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,425,9680TG 8.4% 154.0 / 3.1 69intergenic (‑52/+38)hisQ/hisJlysine/arginine/ornithine ABC transporter/histidine ABC transporter, membrane subunit HisQ/histidine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (24/36);  new base G (6/2);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

ATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTC  >  NC_000913/2425831‑2426110
                                                                                                                                         |                                                                                                                                                
atTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTATGGTCTCCCTCCGGTT‑‑Acca                                                                                                                                           >  2:495713/1‑142 (MQ=255)
    cAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAg                                                                                                                                       <  1:1075138/143‑1 (MQ=255)
     aGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTTCGGGT‑‑AACCGGAgg                                                                                                                                      <  1:207193/143‑1 (MQ=255)
       cGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCGCTTGTCTCCCTCCGGGT‑‑ACCCGAAGGGa                                                                                                                                    >  2:81230/1‑143 (MQ=255)
         cgCACTACAGCGCTGATAGCCCGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAACACCCATCCACCATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGc                                                                                                                                  <  2:447702/142‑1 (MQ=255)
         aGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGCGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCCCCCCTTCTCGCCCTCCGGCT‑‑ACCCCgagcgagc                                                                                                                                  >  2:446181/1‑143 (MQ=255)
              tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGct                                                                                                                             >  1:1009677/1‑143 (MQ=255)
              tACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGCCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCCCCGGGT‑‑ACCCGAACGgggcgggcc                                                                                                                             >  1:417438/1‑142 (MQ=255)
                cAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGcttt                                                                                                                           <  1:746685/143‑1 (MQ=255)
                 agagCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAACCACGTTAAACATTCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGctttc                                                                                                                          <  1:242109/143‑1 (MQ=255)
                   agCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACGGGCTCTGGGCTCCCTCCGGTT‑‑ACCCGGAGGGAGCGGGctttctt                                                                                                                        >  1:726107/1‑143 (MQ=255)
                   agCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑GCCCGAAGGGAGCGGGctttctt                                                                                                                        >  1:65996/1‑143 (MQ=255)
                      tGATAGCCAGCTCCAGCGTGGCGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGATT‑‑ACCCGAAGGGGGCGGGcttccttggc                                                                                                                     >  1:509368/1‑140 (MQ=255)
                       gATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCAc                                                                                                                    <  1:717282/143‑1 (MQ=255)
                            ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTGGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCTCCgggg                                                                                                               >  1:696467/1‑140 (MQ=255)
                            ccAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCAGACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGGGCGGGGTTTCTTTCACCGGTg                                                                                                               >  1:382374/1‑143 (MQ=255)
                             cAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATgg                                                                                                              <  2:581096/143‑1 (MQ=255)
                                cTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAATGGAGCGGGCTTTCTTTCACCGATGGCGa                                                                                                           >  2:253776/1‑143 (MQ=255)
                                     gCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCATCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGc                                                                                                      <  1:464855/143‑1 (MQ=255)
                                      cGTGACGAGCGCACCCTGTAAAATGACTCACCTACAGGAATACAACAGGCTGCCTGTCCTGTCGTGTGTGGTGGATGACCG‑‑TTTTGTTTCTCTTTGGTTTGACCCTCAGCGTGCGTGTTTTCTTTCACAGATGGCGATTGGca                                                                                                     >  2:271833/1‑142 (MQ=255)
                                        ttaCGAGCACTCCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTTTCCTGTGCGGTGTGGTAAATTAGCGCCTCTTGTCTCCCTCCGGTT‑‑AGCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccac                                                                                                   <  2:993100/141‑1 (MQ=255)
                                        tGACGAGCGCACCCTGTAAAATAACACCTGGAAACCCATACGACATGCTGCCTGTCCGGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccac                                                                                                   <  2:897053/143‑1 (MQ=255)
                                        tGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGccac                                                                                                   <  2:575335/143‑1 (MQ=255)
                                              gcgcACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCAAATTATTACCCTCCGGGT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATaaa                                                                                             >  1:493326/1‑143 (MQ=255)
                                              gcgcACCATGTATAAGAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATaaa                                                                                             <  1:532243/143‑1 (MQ=255)
                                                  aCCCTGTCAAATAACGCCTGAAAGCCCATACACCATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATc                                                                                         <  2:1048600/143‑1 (MQ=255)
                                                      tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGGCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGACTAGCCACCATAAACATCaaaa                                                                                     >  2:761179/1‑143 (MQ=255)
                                                      tGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCaaaa                                                                                     >  1:766405/1‑143 (MQ=255)
                                                         aaaaTAACGCCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCg                                                                                  <  1:254745/143‑1 (MQ=255)
                                                         aaaaGAACACCTGAAAACCCATACAACATGCTGTCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCg                                                                                  <  1:28507/143‑1 (MQ=255)
                                                           caTAACACCTGAAAACCCATACAACATGCTGCCTGTCCCGTCGTGGGTGGTGAATTACCGCCACTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  1:397096/142‑1 (MQ=255)
                                                           agtaccacCGGAAAACCCATACAACCGGCTGCCGGTCCTGTCGTGTGTGGTGAATGACCGCCTCTTGTCTCCCTCCGGTG‑‑ACCCGAAGGGAGCGGGCTTTCTTCGACCGATTGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  2:996941/138‑1 (MQ=255)
                                                           aataatATCTGAAAAACCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCTGTT‑‑ACCCGAAGGGCGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  2:89387/143‑1 (MQ=255)
                                                           aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  1:913018/143‑1 (MQ=255)
                                                           aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  1:644001/143‑1 (MQ=255)
                                                           aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGaa                                                                                <  1:727535/143‑1 (MQ=255)
                                                           aaTAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCAATAGCGATTAGCCACCATAAACATCCAAATCGaa                                                                                >  1:34285/1‑143 (MQ=255)
                                                                cacCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACt                                                                           <  2:1063748/143‑1 (MQ=255)
                                                                   cTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACtttt                                                                        >  1:179125/1‑143 (MQ=255)
                                                                     gcaaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTc                                                                      <  2:154337/141‑1 (MQ=255)
                                                                      aaaaCCCATACAGCGTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑CCCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCg                                                                     <  2:426162/143‑1 (MQ=255)
                                                                      aaaaCCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCATCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATCACCATCACACTCAAAGTACTTTTTCg                                                                     >  2:450297/1‑143 (MQ=255)
                                                                                caccaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCGGTCGTGGGCATCCATTCTTTCTCCGATTGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTtctc                                                           <  1:495713/140‑1 (MQ=255)
                                                                                 aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGGTTAGCCACCATAAACGTCCACATCGAAGTACTTTTGCGCCAATTTCTCg                                                          >  2:85150/1‑143 (MQ=255)
                                                                                 aacaTGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCg                                                          <  2:446815/143‑1 (MQ=255)
                                                                                               cctttcTTGCGTGGTGAATTACCGCCTCATGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgc                                            <  1:386594/143‑1 (MQ=255)
                                                                                                 tgtcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgc                                          >  1:517056/1‑143 (MQ=255)
                                                                                                 ggtcGTGTGGGGTGAATGACCGCCTCGTGTGTCCCTCCGGTT‑‑ACCCGTAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgc                                          <  2:823233/142‑1 (MQ=255)
                                                                                                  gtcGGGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                         <  1:864931/143‑1 (MQ=255)
                                                                                                  gccGTGTGGCGTGAAGTACCGCCTCTTGTCTCCCTCCGGTC‑‑ACCCGACGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAgcgcg                                         <  1:281058/141‑1 (MQ=255)
                                                                                                   tcGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACAGTCAgcgcgc                                        >  2:233729/1‑143 (MQ=255)
                                                                                                      tgtgtgGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTACTTTCTCGTAGGTACCGTAAGCGCGCAtt                                     >  2:7824/1‑143 (MQ=255)
                                                                                                           ggtgAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGc                                >  2:459816/1‑143 (MQ=255)
                                                                                                           ggtgAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGc                                <  1:610259/143‑1 (MQ=255)
                                                                                                               caTTACCCCCTCCTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCTGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                            <  1:309426/142‑1 (MQ=255)
                                                                                                               aaTTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCTATTAGTCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAg                            <  1:886711/143‑1 (MQ=255)
                                                                                                                aTTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAgg                           >  1:752401/1‑143 (MQ=255)
                                                                                                                  gaCCGCCTCGTGTGTCGCTCCGGTG‑‑ACCCGAAGGGAGGAGGCTTTCTTTCACCGATGGCGATTAGCCACCGTAAACCTCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCt                         <  2:417438/142‑1 (MQ=255)
                                                                                                                     cGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGCTTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCCAAGGCTTTg                      >  2:681407/1‑143 (MQ=255)
                                                                                                                      gCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGt                     >  1:624586/1‑143 (MQ=255)
                                                                                                                         tctTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCa                  <  2:1009677/143‑1 (MQ=255)
                                                                                                                            ccgctcccttcGGGT‑‑AACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTAGATTTTCGCTAATTTCTCGTAAGTACCGTCGGCGCGCATTTCGGCAAAGGCTTTGTTCAggg               >  2:641606/4‑141 (MQ=255)
                                                                                                                              tctcCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCGTTGTTCAGTGCt             >  1:916027/1‑143 (MQ=255)
                                                                                                                                tccctccGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc           <  2:179125/143‑1 (MQ=255)
                                                                                                                                tccatgcGGTT‑‑ACCCGAAGGGAGCGGGAGTTGTTACACCGAAGGGGAGGAGCCACAAGAAACGTCAAAATCGAAGTGCTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTc           <  1:271833/137‑1 (MQ=255)
                                                                                                                                 cccttcGGGT‑‑AACCGGAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTGGTTCAGTgtgtcg          >  2:770940/1‑138 (MQ=255)
                                                                                                                                  ccttcGGGT‑‑AACCGGAGGGACCGGGCTTTCTTTCACCGATGTCGATTAGCCACCATAAAGCTCAAAATCGACGTACTTTTTCGCTAATTTCTCGTCAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgc         >  2:235316/1‑143 (MQ=255)
                                                                                                                                   ctccGGGT‑‑ACCCGAAGGGAGCGGGCTTTCTTGCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTGTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTcgcg        <  2:645182/143‑1 (MQ=255)
                                                                                                                                     ccGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTGTTTCGCTAATTTCTCGTACCTACAGTTCGCGCCCATTTCGCCAAAGGCTTTGTTCCACGCTTCGCGCa      >  2:220682/1‑143 (MQ=255)
                                                                                                                                         tt‑‑ACCCGAAGGGAGCGGGCTATCTTTCACCGATGGCGATTAGCCACCATAACCATCAAAATCGAAGAACTTTTTAGCTAATTTCTCGTCAGTACCGACAGCGCGCATATCAGCAAGGCCTTTGTTCAGTGCTTCACGCAgata  >  2:713490/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                                
ATTACAGCGAGCACTACAGAGCTGATAGCCAGCTCCAGCGTGACGAGCGCACCCTGTAAAATAACACCTGAAAACCCATACAACATGCTGCCTGTCCTGTCGTGTGTGGTGAATTACCGCCTCTTGTCTCCCTCCGGTT‑‑ACCCGAAGGGAGCGGGCTTTCTTTCACCGATGGCGATTAGCCACCATAAACATCAAAATCGAAGTACTTTTTCGCTAATTTCTCGTAAGTACCGTCAGCGCGCATTTCGGCAAAGGCTTTGTTCAGTGCTTCGCGCAGTTC  >  NC_000913/2425831‑2426110

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: