Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 7,711 A→T A83A (GCT→GCA yaaJ ← putative transporter YaaJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009137,7110AT97.6% 122.7 / NA 43A83A (GCT→GCAyaaJputative transporter YaaJ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (18/23);  minor base G (0/1);  total (19/24)
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.95e-01

CTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGC  >  NC_000913/7473‑7957
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         
cctGAATTCCGCTGAAAATTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTTCCCGCGCAAcc                                                                                                                                                                                                                                              >  1:382952/3‑250 (MQ=255)
   gaattCCGCCGAAAATTATTCCGTCGGCGATGATCAAAAAGCCGGCGAACAGAACCCCCATCCAGCTCATCCCCAGCCCGCGCGCCATATCCCAAGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCTAGCTCAACGCCg                                                                                                                                                                                                                                           <  1:166098/247‑1 (MQ=255)
                aaTTATGCCGTAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGGTGCCGc                                                                                                                                                                                                                              >  1:67150/1‑251 (MQ=255)
                 aTTATGCCGTAGGCAATGAGAAAAAAGCCGGCGAACAGACCGCCCAGCCAGCGCATCCCCAACCCGCGCGCCATATACCAGGCCGGTCCGCCACGAAACTGCCCAGTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGACGCCCTGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCt                                                                                                                                                                                                                             <  1:46497/251‑1 (MQ=255)
                              cGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGcgcc                                                                                                                                                                                                                <  1:154775/251‑1 (MQ=255)
                               gATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGACCAACGCCGGCCAGGGTGCCCCCACCCACGCGcgccc                                                                                                                                                                                                               >  1:284770/1‑250 (MQ=255)
                                                          ccccctCCAGTGCACCCCCAGCCCGCGCGCCATATCCCACTCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAACCGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATg                                                                                                                                                                                   <  1:565219/246‑1 (MQ=255)
                                                               ccccgccaacccccgccccGCGCCCCATATCCCATCACGGTCCGCCAAGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGGCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAGTGCCAGCGCAACGCCGGCCAGGGTGCCGCGACCCACGCGCGTCGCAAGACTGGTACACAATGACTGAAATGAGGtt                                                                                                                                                                             <  2:544714/235‑1 (MQ=255)
                                                                ccAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCGCGCAACGCCCGCCCGGGTGCCGCTACCCACGCGCGCCGCCATGCTGGGACACAATGaagatcggaagcgcc                                                                                                                                                                             >  1:480280/1‑237 (MQ=255)
                                                                    caccacccccccccGCGCCCCATATCCCATGCGGGTCCGCAACAAAACTGCCCATTGGCACCACGTGCTTTATAAAGCTGTGGGAGGGACCATCCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACAGCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAAcc                                                                                                                                                                        <  2:629151/240‑1 (MQ=255)
                                                                     gcATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCACGCAACGCCGGCCCGGGTGCCGCTAACCACGCGCGCCGCCAGACTGGGACACAAAGACTGAAATGAGGTTaaaacc                                                                                                                                                                       >  1:20508/1‑248 (MQ=255)
                                                                            cAGCCCGCGCGCCATATCCCATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGTTACACAATGACTGAAATGAGGTTAAACCGCCTGGCt                                                                                                                                                                <  1:571211/251‑1 (MQ=255)
                                                                              ccccgcgccccaaaacacATGCCGGTCCACCACAAAGCTGCCCGTTGGCATCAGTTGCTTTATAAGGCGGTGCGGGGGAACATCCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCACCCCACATCCAGAAGACGGCACCAGGTCAACCGGCGGTAATTGCCAGCGTAACGCCGGCCAGGTGGCCGCTCCCCACGCGCGCCGAAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCtgt                                                                                                                                                              <  2:107214/234‑1 (MQ=255)
                                                                                       ccATATCCCATGCCGGCCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGAAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATg                                                                                                                                                     <  1:603899/251‑1 (MQ=255)
                                                                                        cATAGACCATGCCGGGCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATGCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGAc                                                                                                                                                                                                         <  1:86826/200‑1 (MQ=255)
                                                                                           ataCCATGCCGTCCCCCCACGAACCTGCCCATGGGCATCACGTTCTTAATAAACCGGTGCGAGGTACCATCCGCCAAACGAGGTCGCCTTGCCGATAAACGCGGCAACCCACACCCAGAATACGGCTCCAGGTCCACCGGCGGTAATGGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCCCGCGCGCCGTAAGACTTGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTAt                                                                                                                                                 <  1:340879/251‑1 (MQ=255)
                                                                                               cATGCCGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAAATTCCAGCGCCAACCCGGACAGGGTGCCGCCACCCCCACGGCCCGCAAGGCTGGGCAACACTGAATGAAATGAGGGTAAACCCCCTGGCTGTGGATGAATGCTCttttt                                                                                                                                             >  2:415830/1‑251 (MQ=255)
                                                                                                    cGGTCCGCCACGAACCTGCCCATTGGCATCCCTTCCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGAc                                                                                                                                        <  1:144357/251‑1 (MQ=255)
                                                                                                    cGGTCCGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTa                                                                                                                                                                            <  1:668087/215‑1 (MQ=255)
                                                                                                         cGCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGAAATTGCCAGCGCAAAGCCGGCCAGGATGCCGCTACCCACGCGCGCCGCCAGACTGGTACACCATGACTGGAATGAGGGTAAACCGCCTGGCTTTGGGTGTATGCTATTTTTAAggcctttg                                                                                                                                   >  1:480929/1‑244 (MQ=255)
                                                                                                          gCCACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCCTCGCAACGCCGGCCAGGGTGCCGCCACCCCCCCGCGCCGCATGACTGGGACACAATGGCTGAACAGAGGTTAAACCGCACGGCTTGGGATGAATGATCTTTTTAAAACtttttc                                                                                                                                  >  2:526969/1‑249 (MQ=255)
                                                                                                            cACGAAACTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCCGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGGACACCATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGGTGAATGCTATTTTTAAGACTTTTGccc                                                                                                                                >  1:364508/1‑250 (MQ=255)
                                                                                                               caaaaccgcccattggcatcccttacTTTATAAAGCTGTGCGAGGGAACATGCGGCAAATGGGGTCGCCAGGCCGAAAAACGGGGCAACCCACACCCAGAAGACGGCTCCTGGTCCACGGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAGGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCaaa                                                                                                                              <  2:675250/226‑1 (MQ=255)
                                                                                                                  aCTGCCCATTGGCATCACGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGGAATTGCCAGCGCAAAGCCGGCCAGGGTGCCGCTACCCCCGCGCGCCGCAAGGCTGGTACACAATGACTGAAATGAGGTTAAACCGCCGGGGTGTGGGAGAATGCTATTTTTAAACCTTTTGCCAAACTgg                                                                                                                          >  1:227743/1‑251 (MQ=255)
                                                                                                                                   cGTTCTTTATAAAGCTGTGCGAGGGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAAATGCCCGCCCAACGCCGGCCAGGTTGCCGCTACCCCCGCTCGCAGCAAGACTGGACAACAATGCTTGAAATTGGGTTTAACCGCCTCGATGTGGGTGAATGCTCATTTTAACACTTTTTCCAAACTGGCggggctagaggaaaa                                                                                                          >  2:7721/1‑237 (MQ=255)
                                                                                                                                          ttATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTTATAGCCAGCCCCACCCCGGCCAGGGTGCCGCCACCCACCCGCGCCGCCAGAACGGGACACAATAGCTGCAATTAGGGTTAAACCCCCGGCTGGGGAGGAAGGCGATTTTTAAAAAATTTGACCAACTGGCgggggggaaggaaacacacac                                                                                                    >  2:120823/1‑231 (MQ=255)
                                                                                                                                                 gctgtgCGAGGGAACATGCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCCCCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTg                                                                                            <  1:40028/249‑1 (MQ=255)
                                                                                                                                                            gaACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGGTGCCGCTACCCCCGCGCGCCGCCAGACTGGTACACAATGACTGAAATGAGGGTAAACCGCCCGGCtgtg                                                                                                                                                             >  1:588181/1‑175 (MQ=255)
                                                                                                                                                             cgccgatctcAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCCCATCCATAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCATGTTGCCGCTACCCACGCGCGCCGCAAGACTGTTACACAATGACTGAAATGAGTTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAAc                                                                                <  1:25512/242‑1 (MQ=255)
                                                                                                                                                                  gaggaaaccGAGGTCGCCTGTCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCAACCGGCGGTAATTGCTAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACGGGTACACAAGGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAACCTGCACAAATCCGGTGCGAAAAGTGAACcaaca                                                                           <  1:514633/243‑1 (MQ=255)
                                                                                                                                                                    ggCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTTCCAGCGCAACGCCGGCCAGGATGCCCCTACCCCCCCGGGCCGCAAGACTGGTAAACAATGAATGAAAAGGGGTTAAAACGCCTGGATGTGGATGAATGCTATTTTTAAGGCTTTTTCCAAACTGGCGGATGGAACGAAAATGGAAAAATTCGGGGCGAAAAgcggaaacacaa                                                                          >  1:335089/1‑239 (MQ=255)
                                                                                                                                                                         aCGAGGCCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCCCCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGATGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTgcg                                                                    <  1:390048/251‑1 (MQ=255)
                                                                                                                                                                             ggTCGCCATGCCGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGGTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTTCCAAACTGGCGGGTGTAGCGAAAATGCACAAATCCGGTGCGAAAAAGGAAACAACAAACTggcccgg                                                                >  1:192887/1‑245 (MQ=255)
                                                                                                                                                                                cgcaaTGCCGATAAACGCGGCAACCCACACCCAGAAAACGGCTCCAGGCCCACCGGCGGTAATTGCCAGCCCACCGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAaga                                                             <  1:634138/247‑1 (MQ=255)
                                                                                                                                                                                cGCCATGCCGATCACCGCGTAAACCCCCAGCCATAAGACGTCTCCAGGTCCCCCGGCTGTAATTGCCAGCGCAACGCCGGCCAGGTTTCCGCTACCCACGCGCGCCGCAAGACTGGTACACAAGGTCTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACGGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAaga                                                             <  1:630165/251‑1 (MQ=255)
                                                                                                                                                                                  ccaggCCGATAAACGCGGCACCCCACATCCAGAAGACGGATCCAGGTCCACCGCCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCTCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGATCCAACAACCTGCGCCGAAGAGc                                                           <  1:540304/247‑1 (MQ=255)
                                                                                                                                                                                        ccGATAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGGTGCCGCCACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGGTAAAACGCCTGGCTGTGGGTGAATGCTTTTTTTAAGACTTTTGCCAAAATGGAGGGTGTAGCGAAACAGCACAAATCCCGTGCGACAAGTGAACCAACAACCCGGGCCGAAGAGCAggc                                                       >  1:116513/1‑249 (MQ=255)
                                                                                                                                                                                            tAAACGCGGCAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGGCAGGGTGGCGCCCCCCCCCCGCGCCGCAAGACTGGTACACCATAACTGAAATGAGGTTAAACCACCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAAATGGCGGATGGTGCTAAACTGCACACATCCGGTGAGAAAAAGGAACCAAAAAACTGCGACGCAGAGAAGGAAAATc                                                  >  1:662714/1‑251 (MQ=255)
                                                                                                                                                                                                     cAACCCACATCCAGAAGACGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGAt                                         <  1:187101/251‑1 (MQ=255)
                                                                                                                                                                                                                       cGGCTCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGGGGGTAAACCGCCTGGCTGTGGGTGAATGCTATTTTTAAGAATTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGGGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGGAAATCATTACCGATACCCAAAGGACGCTGtt                        >  1:522018/1‑250 (MQ=255)
                                                                                                                                                                                                                           tCCAGGTCCACCGGCGGTAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCaacaac                                                                         <  1:387766/197‑1 (MQ=255)
                                                                                                                                                                                                                                gtcacccGGGGGTAATTGCCGGCGCACGGCCGGCCAGGTGGCAGCTCCCCCCCCGCGCTGCAAGACGGGTACACAATGAGGTAAATGAGTGTAACCCGCCGGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGTGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGgag              <  2:283527/246‑1 (MQ=255)
                                                                                                                                                                                                                                            tAATTGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGGGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATAAAGGGGGAAAAATCTggg  >  2:663552/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         
CTTGAACTCCGCTGAAAATTATGCCATAGGCGATGAGCAAAAAGACGGCGAACAGAACGCCCATCCAGCGCATCCCCAGCCCGCGCGCCATATACCATGCCGGTCCGCCACGAAACTGCCCATTGACGTCACGTTCTTTATAAAGTTGTGCCAGAGAACATTCGGCAAACGAGGTCGCCATGCCGATAAACGCGGCAACCCACATCCAAAAGACGGCTCCAGGTCCACCGGCGGTAATAGCCAGCGCAACGCCGGCCAGGTTGCCGCTACCCACGCGCGCCGCAAGACTGGTACACAATGACTGAAATGAGGTTAAACCGCCTGGCTGTGGATGAATGCTATTTTTAAGACTTTTGCCAAACTGGCGGATGTAGCGAAACTGCACAAATCCGGTGCGAAAAGTGAACCAACAACCTGCGCCGAAGAGCAGGTAAATCATTACCGATCCCCAAAGGACGCTGTTAATGAAGGAGAAAAAATCTGGC  >  NC_000913/7473‑7957

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