Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 717,368 A→C 6.7% F300V (TTT→GTT)  potE ← putrescine/proton symporter: putrescine/ornithine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913717,3680AC 6.7% 189.9 / 6.0 71F300V (TTT→GTT) potEputrescine/proton symporter: putrescine/ornithine antiporter
Reads supporting (aligned to +/- strand):  ref base A (34/31);  new base C (3/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AACTGACTGTTTAGCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCAT  >  NC_000913/717220‑717512
                                                                                                                                                    |                                                                                                                                                
aaCTGACTGTTTAGCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTaaa                                                                                                                                                  <  6:186433/149‑1 (MQ=255)
     actgTTTAGCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCt                                                                                                                                             >  5:114440/1‑149 (MQ=255)
     actgTTTAGCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCt                                                                                                                                             <  7:138382/149‑1 (MQ=255)
             gCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACAGGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGCTTTAAACACCTGGGCAATg                                                                                                                                     >  2:342391/1‑149 (MQ=255)
             gCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATg                                                                                                                                     <  5:159838/149‑1 (MQ=255)
             gCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATg                                                                                                                                     <  4:223084/149‑1 (MQ=255)
                aCGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGt                                                                                                                                   >  1:109548/1‑148 (MQ=255)
                  gggCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGaa                                                                                                                                >  5:61715/1‑149 (MQ=255)
                   ggCTAATGGGCATCAGCGACAAACCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACATTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTAAAACACCTGGGCAATGGTCAAc                                                                                                                               <  5:107027/149‑1 (MQ=255)
                    gCTAATGGTCATCAGTGCGAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACt                                                                                                                              <  5:76162/149‑1 (MQ=255)
                    gCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACt                                                                                                                              >  1:155/1‑149 (MQ=255)
                           gTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTgccagc                                                                                                                        >  3:206121/1‑148 (MQ=255)
                              atcaGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAg                                                                                                                    >  3:176591/1‑149 (MQ=255)
                               tcaGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGt                                                                                                                   <  2:109548/149‑1 (MQ=255)
                                    gCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCg                                                                                                               >  3:60768/1‑148 (MQ=255)
                                     ccAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGGAATGGTGAACTGCCAGCCAAGTAGCGaa                                                                                                             >  3:13328/1‑149 (MQ=255)
                                       aaTCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAAcc                                                                                                           >  3:289033/1‑149 (MQ=255)
                                            aCTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCgcagc                                                                                                      <  4:13328/149‑1 (MQ=255)
                                              tctGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGcagcag                                                                                                    >  2:70215/1‑149 (MQ=255)
                                                tGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTgg                                                                                                                                           <  5:37764/108‑1 (MQ=255)
                                                tGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTgg                                                                                                                                           >  6:37764/1‑108 (MQ=255)
                                                  aataatCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGAc                                                                                                 <  4:20811/148‑1 (MQ=255)
                                                    taatCACAATGGTCACCATTCCCTGCACAGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACATGGGCAATTGTGAAATGCCAGCCAAGTAGCGAACCGCAGCAGGACAt                                                                                               >  6:314856/1‑148 (MQ=255)
                                                             tGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAgc                                                                                     >  2:311397/1‑149 (MQ=255)
                                                             tGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAgc                                                                                     <  4:60768/149‑1 (MQ=255)
                                                               gTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAACCACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAgcgc                                                                                   >  8:1780/1‑149 (MQ=255)
                                                               gTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAACCACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAgcgc                                                                                   >  8:6542/1‑149 (MQ=255)
                                                                  aaCATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAACCACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACAATCAGCGCCa                                                                                 >  1:270595/1‑148 (MQ=255)
                                                                  aaCATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCa                                                                                 >  5:160577/1‑148 (MQ=255)
                                                                      ttCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAAt                                                                             >  3:28060/1‑148 (MQ=255)
                                                                         ccTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACt                                                                         >  6:40034/1‑149 (MQ=255)
                                                                          cTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACtt                                                                        >  3:262749/1‑149 (MQ=255)
                                                                             cACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTa                                                                      <  1:70215/148‑1 (MQ=255)
                                                                              aCCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTAcc                                                                    >  1:36472/1‑149 (MQ=255)
                                                                              aCCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTAcc                                                                    >  8:84651/1‑149 (MQ=255)
                                                                                cGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAACCACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCa                                                                  <  7:1780/149‑1 (MQ=255)
                                                                                cGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAACCACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCa                                                                  <  7:6542/149‑1 (MQ=255)
                                                                                  gTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGTAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCACCAGGACATCACCATCAGCGCCATAATGACTTTACCCACt                                                                >  6:21372/1‑149 (MQ=255)
                                                                                     cATCCACTTTGGGTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTcc                                                             >  6:30813/1‑149 (MQ=255)
                                                                                         cACTTTGGTTACACGTGAGAAAATTTAAGGGAAGTAGCCATCATCAGATGAAGATTTAATCACCTGGGCAAGGGTGAACTGGCAGCCAAGTAGCGATCCGCTGCAGGTCATCACCATCAGCGCCATAATGACTCTACCCACTTCCCGc                                                          <  4:279122/148‑1 (MQ=255)
                                                                                         cACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGc                                                          <  2:270595/148‑1 (MQ=255)
                                                                                            tttGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGATCCGCAGCAGGACATCACCATCAGCGCCATAAAGACTTTACCCACTTCCGGCGTg                                                       >  1:125757/1‑148 (MQ=255)
                                                                                              tGGTTAAACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGaa                                                     <  2:347755/148‑1 (MQ=255)
                                                                                                gTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACAt                                                  <  5:40034/149‑1 (MQ=255)
                                                                                                    cacGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGATCCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTgc                                              <  5:314856/149‑1 (MQ=255)
                                                                                                    cacGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTg                                               <  6:160577/148‑1 (MQ=255)
                                                                                                      cGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTgcgc                                            <  4:289033/149‑1 (MQ=255)
                                                                                                          agaAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAg                                        >  7:232900/1‑149 (MQ=255)
                                                                                                            aaaaTTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGc                                      >  7:265449/1‑149 (MQ=255)
                                                                                                               aTTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAg                                   <  2:36472/149‑1 (MQ=255)
                                                                                                                     gggAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCaaa                             <  5:21372/149‑1 (MQ=255)
                                                                                                                       gAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATg                           <  8:108760/149‑1 (MQ=255)
                                                                                                                             gCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCg                     >  3:185849/1‑149 (MQ=255)
                                                                                                                               cTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGtt                   >  2:133032/1‑149 (MQ=255)
                                                                                                                               cTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGtt                   <  1:133032/149‑1 (MQ=255)
                                                                                                                               cTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGtt                   <  4:176591/149‑1 (MQ=255)
                                                                                                                                 tcaacaGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTg                  >  1:210435/1‑148 (MQ=255)
                                                                                                                                   atcaGATGAAGAATTAAACACCTGGGCAATGGTGCACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCATGAACATATGCGCGAAGGCCAGACCAAATGGTGCCGTTGAAt               >  5:297809/1‑149 (MQ=255)
                                                                                                                                    tcaGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAAtt              <  8:265449/149‑1 (MQ=255)
                                                                                                                                      aGATGAAGATTTAAACACCTGGTCAATTGTTAACTGCCATCCAAGTAGAGAACCGAAGCATTACATCACAATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATTGTGCCGTTGAAttt             <  2:155/148‑1 (MQ=255)
                                                                                                                                      aGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAAttt             <  6:114440/148‑1 (MQ=255)
                                                                                                                                      aGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTg            >  7:118770/1‑149 (MQ=255)
                                                                                                                                       gATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGc           <  1:311397/149‑1 (MQ=255)
                                                                                                                                       gATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGGGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGc           <  4:262749/149‑1 (MQ=255)
                                                                                                                                         tGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCa         <  4:206121/149‑1 (MQ=255)
                                                                                                                                         tGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCa         <  1:342391/149‑1 (MQ=255)
                                                                                                                                          gAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAACGCCAGACCAAATGGTGCCGTTGAATTTGCCAg        >  4:298135/1‑149 (MQ=255)
                                                                                                                                          gAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACAATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAGATGGTGCCGTTGAATTTGCCAg        >  8:340246/1‑149 (MQ=255)
                                                                                                                                           aaGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGc       >  5:225109/1‑149 (MQ=255)
                                                                                                                                             gATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGAcc                                >  4:121625/1‑122 (MQ=255)
                                                                                                                                             gATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGAcc                                <  3:121625/122‑1 (MQ=255)
                                                                                                                                                 tAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCat  <  8:200549/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AACTGACTGTTTAGCGACGGGCTAATGGTCATCAGTGCCAATCCACTCTGAATAATCACAATGGTCAACATTCCCTGCACCGGTGCATCCACTTTGGTTACACGGGAGAAAATTTTAGGGAAGTAGCCTTCATCAGATGAAGATTTAAACACCTGGGCAATGGTGAACTGCCAGCCAAGTAGCGAACCGCAGCAGGACATCACCATCAGCGCCATAATGACTTTACCCACTTCCGGCGTGAACATCTGCGCGAAGGCCAGACCAAATGGTGCCGTTGAATTTGCCAGCTCCAT  >  NC_000913/717220‑717512

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