Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC95.9% 141.7 / ‑5.6 49Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/24);  minor base A (0/2);  total (23/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCG  >  NC_000913/1358614‑1359081
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
aGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGGTAATCATCGGCAAAAAGAGACAATcc                                                                                                                                                                                                                           >  1:395122/1‑251 (MQ=255)
  aCGGTGCCAACGCCTTCTCCTTTAGGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGTGAAGAGGGCGGAGAGATAGATGTTCGCCAACAAGACGGTAATCATCGGCAAAAATAGACAAGCcag                                                                                                                                                                                                                         <  2:631273/251‑1 (MQ=255)
   cGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGCAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCCCCCACGGCGGCCCTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACCAGACGTTCATCAACGGGCAAAAGAGACAATCcagc                                                                                                                                                                                                                        >  2:33909/1‑251 (MQ=255)
      tGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACCCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCagcagc                                                                                                                                                                                                                     >  1:418170/1‑251 (MQ=255)
       gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCg                                                                                                                                                                                                                    <  2:147717/251‑1 (MQ=255)
            cGCCTTCTCCTTTATGCAGTCCCTGAACTAACAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGAAGAAGGTTACCACCCACACCCACGGCGGCACTTTCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACCAGACGTTAATCATCGGCCAAAAGAGACAATCCAGCAGCGATGAc                                                                                                                                                                                                               >  1:201288/1‑251 (MQ=255)
             gCCTTCTCCTTTATGCAGGCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTGCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGAcc                                                                                                                                                                                                              <  2:418170/251‑1 (MQ=255)
                 tctcCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAg                                                                                                                                                                                                            <  1:250691/249‑1 (MQ=255)
                 tctcCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAg                                                                                                                                                                                                            <  1:293980/249‑1 (MQ=255)
                    ccTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTCACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATGTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGcc                                                                                                                                                                                                          <  2:278349/248‑1 (MQ=255)
                             aGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGCAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCCACCATAACTcc                                                                                                                                                                                              >  2:30197/1‑251 (MQ=255)
                              gTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCg                                                                                                                                                                                             <  1:598552/251‑1 (MQ=255)
                                    gAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGc                                                                                                                                                                                       <  2:546611/251‑1 (MQ=255)
                                                  aTAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCCCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAAATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGCCCATAAATCCGACGTGCGGGTTAATCGCCtt                                                                                                                                                                         >  1:264555/1‑251 (MQ=255)
                                                    aaaCACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACtttt                                                                                                                                                                       <  2:670309/251‑1 (MQ=255)
                                                       caccacCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGAAGGTTAACACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGGGATAGATTTTCGCCAACAAGACGTTAATCAACGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGAATTTTgcg                                                                                                                                                                    >  1:595087/1‑251 (MQ=255)
                                                         ccaccaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCACCCCCACCCACAGCGGCACTTCCGGGGAGAAGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATAGGCAAAAAGAGACAATCCCGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTTcgcg                                                                                                                                                                  >  1:442968/1‑251 (MQ=255)
                                                                     agaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGcc                                                                                                                                                      <  1:46790/251‑1 (MQ=255)
                                                                       aTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCCCCCACAACCCCGGCGGCACTTCCGGGAGGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCCGCCGCGGTGACCAACCGAACATCAATCCGACGTGCGGGTTAATTGACTTTTGCGCGGTGGGATAGGCCgg                                                                                                                                                    >  2:87334/1‑250 (MQ=255)
                                                                       aTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAgg                                                                                                                                                              <  1:246128/241‑1 (MQ=255)
                                                                       aTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTTCCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGGAAAAAGAGACAAACCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGCATTTTGCGCGGAgg                                                                                                                                                              >  2:246128/1‑241 (MQ=255)
                                                                                   aCAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTTCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCg                                                                                                                                         <  1:318412/250‑1 (MQ=255)
                                                                                     aaaCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGGCGCAGGGTACCACCCACCCCCACGGGGGCCCTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCCACAAGAAGGTAAACATCGGGAAAAAGCGACAATCCAGCAGCGATCACCCGCCGACCATAAATCCGACGCGCGGGCTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGa                                                                                                                                       >  1:670162/1‑250 (MQ=255)
                                                                                     aaaCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACCCCCCCGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCCACCATAAATCCGACGCGCGGGTTAATCGACTTTTGCGCGGAGGTATAGGCCGAACCGGCCTCCGGa                                                                                                                                       >  1:666040/1‑250 (MQ=255)
                                                                                          gggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTTCCAACCCCACCCACGGCAGCACTTCCCGGAAGAGGGCGGGGAGATAGATTTTCGCCAACAAGAAGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTTCGCGTAGGTATAGGCCGAACCGGCCTCCGGGAACttg                                                                                                                                 >  1:85690/1‑249 (MQ=255)
                                                                                          gggTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTgg                                                                                                                                 <  2:148115/251‑1 (MQ=255)
                                                                                               ttGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAACCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGCACCGGCCTCCGGAAACTGGCGAAc                                                                                                                            >  1:135011/1‑251 (MQ=255)
                                                                                                  aaGTTAGCGACCAGGTTTACGCTCTTCAGACTCGCTGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGTAAGAGGGCGGAGAGATAGATGTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGGCCATAAATCCGACGTGCGGTTTAATCGACTTTTCCGCGTAGGTATAGGCCGAACCGGCGTCCGGAAATTGGCGAACCAg                                                                                                                         <  2:358885/251‑1 (MQ=255)
                                                                                                           aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAAACCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGGATAGGCCGGACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGt                                                                                                                 >  1:254157/1‑250 (MQ=255)
                                                                                                                  tGACGCTCTTCAGATTCGCGGCGGATAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGCGGAGAGATAGATTTTTGCCAAAAAAACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAACACCGACGTGCGGGTTAATCGAATTTTGCGCGGAGGTAGAGGCCGGACCGGCCTCCGGGAACTGGCGAACCAGTTTGCCGTAGCTggt                                                                                                          >  1:628982/1‑248 (MQ=255)
                                                                                                                   gACGCTCTTCAGATGCGCGGCGGTTAAAATTGCGACGAAGGTTACCACCCACTCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATGTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCATCAGCGATTACCAGCCGACCATAAAGCCGACGTGCGGGTTAATCGACTTTTGCGCGTTGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTGTGCCGTAGCTGATat                                                                                                        <  2:429398/251‑2 (MQ=255)
                                                                                                                         ctTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGt                                                                                                    <  1:33909/249‑1 (MQ=255)
                                                                                                                            cAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACa                                                                                               <  1:213599/251‑1 (MQ=255)
                                                                                                                                                        aaGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAg                                                                                                               <  1:442717/207‑1 (MQ=255)
                                                                                                                                                        aaGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAg                                                                                                               >  2:442717/1‑207 (MQ=255)
                                                                                                                                                            ttACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGcc                                                               >  1:263858/1‑251 (MQ=255)
                                                                                                                                                                cacccacacccaCGGGGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGGGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCTAACCGGCCTCCGAAAACTGGCGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGaa                                                           <  1:276404/251‑1 (MQ=255)
                                                                                                                                                                                   ccTTCCGGGAAGAGGGCGGAGAGAAAGATTTCCGCCAACAAGACGTTAATCATCGGAAAAAAGAGACAATCCAGAAGCGAGGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTGGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAAACTGGCAAACCAGTTTGCCGTGGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGcc                                        <  1:235503/250‑1 (MQ=255)
                                                                                                                                                                                        cGGGAAGGGGGCGGAGAGATAGATTTGCGCCAACAAGACGTTAATCACCGGCAAAAAGAGACAACCCAGCTGCGATGACCAGCCGACCATAAACCCGACGGGCGGGTTAATCGCCTTTTGCGCTTAGGTATAGGCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgaga                                   <  2:392218/251‑1 (MQ=255)
                                                                                                                                                                                         gggAAGAGGGGGGAGAGATAGATTTCCGCCAACAAGACGTTAACCAACGGAAAAAAGAGAAAACCCAGCGGCGAGGACCAGCCGACCATAAATCCGACGTGGGGGGTAAGCGCCTTTTGCGCGAAGTTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTGGCCGTGGCTGATAGGGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCTCTAATGCCGGAGAc                                  <  2:666040/251‑1 (MQ=255)
                                                                                                                                                                                          ggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACa                                 <  1:30197/251‑1 (MQ=255)
                                                                                                                                                                                            aaGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAACCATCGGCAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCCGTCGTGCGGGTTAATCGACTTTTGCTCGTAGGTAGAGGCCGAACGGGCCCCCGGAAACTGGCGACCCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                               <  1:111280/251‑1 (MQ=255)
                                                                                                                                                                                                 ggCGGATAGATCGATTTTCGACAACAAGACGGTACTGAACGGACGAAAGAGACAATCCCGCAGCGATGACCCGCCGACCGTAAGTCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACAGGCCTCCGGAAACTGGCGAAACAGTTTGCCGTAGCTGATAGCGGAAAACAGCACACCCGCCAGCGCCAGCAACTAGGACGCCGGAACGCGACCGTCGCTAATGCCGGGGACAATaacc                           >  2:81584/1‑247 (MQ=255)
                                                                                                                                                                                                    ggAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGTCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCAGGAGACAATACCaaa                         >  1:14823/1‑249 (MQ=255)
                                                                                                                                                                                                                  ccGCCAACAAGACGTAAATCACGGGAAAAAAGAGACAATCCAGCATCGATGACCAGCCGACCAAAAACCCGCCGTGCGGGTAAACCGACTTTTGCGCGTAGGTAGAGGCCGAGCCGGCCTCCGGAACCTGGCGAACAGGTTGGCCGTAGCTGAGAGCGGAAAACAGCCCACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGCCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTc         <  2:670162/250‑1 (MQ=255)
                                                                                                                                                                                                                   cGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTc         >  1:522516/1‑250 (MQ=255)
                                                                                                                                                                                                                        aaaaaccgttaacaaccggaaaaaaGAGAAAATCCAGCAGCGAGGACCAGCCGACCATAAACCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAACCTGGCGAACCAGTTTGCCGTAGCGGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGc   <  2:395122/231‑1 (MQ=255)
                                                                                                                                                                                                                         caaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCg  >  1:95563/1‑251 (MQ=255)
                                                                                                                                                                                                                         caaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGAGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCg  >  1:327444/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                              
AGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCG  >  NC_000913/1358614‑1359081

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