Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,438 A→T 5.7% F34Y (TTT→TAT) ‡ phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,4380AT 5.7% 514.3 / 14.9 141F34Y (TTT→TAT) ‡phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (64/69);  new base T (5/3);  total (69/72)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 145 total aligned reads displayed.

CCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACC  >  NC_000913/4315292‑4315582
                                                                                                                                                  |                                                                                                                                                
ccaACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAc                                                                                                                                                <  2:1010631/149‑1 (MQ=255)
  aaCTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgc                                                                                                                                              >  1:51039/1‑149 (MQ=255)
     tCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgc                                                                                                                                              <  4:862578/146‑1 (MQ=255)
     tCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtg                                                                                                                                           <  2:811567/149‑1 (MQ=255)
       cTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggt                                                                                                                                         <  3:340565/149‑1 (MQ=255)
        tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAGTGCAGGCCGATCATGCCGTCAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggtg                                                                                                                                        >  3:1235171/1‑149 (MQ=255)
        tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                        <  1:1248816/149‑3 (MQ=255)
        tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggtg                                                                                                                                        <  4:704063/149‑1 (MQ=255)
        tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgt                                                                                                                                            <  2:168588/145‑1 (MQ=255)
        tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgt                                                                                                                                            >  1:168588/1‑145 (MQ=255)
           aTTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                        >  3:1241317/1‑144 (MQ=255)
                gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACt                                                                                                                                >  3:54993/1‑149 (MQ=255)
                  cGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTcc                                                                                                                              <  2:51039/149‑1 (MQ=255)
                   gATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCg                                                                                                                             >  4:580824/1‑149 (MQ=255)
                         gTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCGTGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGaaa                                                                                                                                                 >  1:108443/1‑123 (MQ=255)
                            gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                    <  3:1072679/149‑1 (MQ=255)
                             aCATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAg                                                                                                                   <  1:430929/149‑1 (MQ=255)
                                 gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAc                                                                                                                                 >  4:854294/1‑131 (MQ=255)
                                 gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAc                                                                                                                                 <  3:854294/131‑1 (MQ=255)
                                 gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                  <  3:1264685/146‑1 (MQ=255)
                                 gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGATGGTCAAACTCCGCCTGTTTTAGc                                                                                                                  >  4:1264685/1‑146 (MQ=255)
                                       gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAc                                                                                                                                                >  1:1333243/1‑110 (MQ=255)
                                       gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGt                                                                                                                        >  4:38454/1‑134 (MQ=255)
                                       gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                    >  3:260037/1‑138 (MQ=255)
                                       gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATc                                                                                                         >  3:646099/1‑149 (MQ=255)
                                          gAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggtg                                                                                                                                        <  1:853115/115‑1 (MQ=255)
                                               tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGaaa                                                                                                                                                 <  1:43575/101‑1 (MQ=255)
                                               tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGaaa                                                                                                                                                 >  2:43575/1‑101 (MQ=255)
                                               tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgc                                                                                                                                              >  1:390478/1‑104 (MQ=255)
                                               tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgc                                                                                                                                              >  1:248048/1‑104 (MQ=255)
                                               tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgc                                                                                                                                              <  2:248048/104‑1 (MQ=255)
                                                gCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCaaa                                                                                                                                  <  4:763007/115‑1 (MQ=255)
                                                  aaaTGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAc                                                                                              >  2:1316954/1‑149 (MQ=255)
                                                   aaTGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                    <  3:738734/126‑1 (MQ=255)
                                                          gCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCgg                                                                                      >  1:876945/1‑149 (MQ=255)
                                                             gATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtg                                                                                                                                           >  3:935425/1‑93 (MQ=255)
                                                             gATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtg                                                                                                                                           <  4:935425/93‑1 (MQ=255)
                                                               tCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGt                                                                                                                        <  4:521645/110‑1 (MQ=255)
                                                               tCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGt                                                                                                                        >  3:521645/1‑110 (MQ=255)
                                                               tCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCg                                                                                 >  2:1136033/1‑149 (MQ=255)
                                                                         aaCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGt                                                                                         <  3:1424722/131‑1 (MQ=255)
                                                                         aaCTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGt                                                                                         >  4:1424722/1‑131 (MQ=255)
                                                                            ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                      <  4:1010251/115‑1 (MQ=255)
                                                                            ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                      >  3:1010251/1‑115 (MQ=255)
                                                                            ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                      <  2:1352503/115‑1 (MQ=255)
                                                                            ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCg                                                                    >  1:539849/1‑149 (MQ=255)
                                                                             tCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCgg                                                                   <  3:83783/149‑1 (MQ=255)
                                                                              cGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGt                                                                         <  2:145025/142‑1 (MQ=255)
                                                                                ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtggg                                                                                                                                         >  3:473133/1‑75 (MQ=255)
                                                                                 cATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                    >  3:951144/1‑96 (MQ=255)
                                                                                     aaGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCt                                                           <  1:1316954/149‑1 (MQ=255)
                                                                                      agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTc                                                          >  3:814235/1‑149 (MQ=255)
                                                                                      agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGCTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTc                                                          >  2:879144/1‑149 (MQ=255)
                                                                                       gcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTca                                                         >  3:814963/1‑149 (MQ=255)
                                                                                              cAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAg                                                  <  2:876945/149‑1 (MQ=255)
                                                                                                gATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGc                                                >  1:752344/1‑149 (MQ=255)
                                                                                                   ggTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                  >  4:520229/1‑112 (MQ=255)
                                                                                                   ggTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                  <  3:76951/112‑1 (MQ=255)
                                                                                                   ggTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                  >  4:76951/1‑112 (MQ=255)
                                                                                                       gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAcc                                                                                             <  1:1256749/97‑1 (MQ=255)
                                                                                                       gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtgt                                         <  2:521948/149‑1 (MQ=255)
                                                                                                       gcgcATGGTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAcc                                                                                             >  2:1256749/1‑97 (MQ=255)
                                                                                                         gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGttt                                       >  1:1007894/1‑149 (MQ=255)
                                                                                                          cATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgtg                                                                                    >  3:376876/1‑103 (MQ=255)
                                                                                                          cATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgtg                                                                                    <  4:376876/103‑1 (MQ=255)
                                                                                                          cATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                      <  2:1007894/149‑1 (MQ=255)
                                                                                                             gTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTcc                                   <  1:50414/149‑1 (MQ=255)
                                                                                                              tttGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                  >  3:1166230/1‑149 (MQ=255)
                                                                                                                   gTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGa                             <  4:1166230/149‑1 (MQ=255)
                                                                                                                    tCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGAt                            <  2:429193/149‑1 (MQ=255)
                                                                                                                       cgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                         <  4:646099/149‑1 (MQ=255)
                                                                                                                       cgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                         <  4:814963/149‑1 (MQ=255)
                                                                                                                        gcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACt                                                                             >  3:1450893/1‑96 (MQ=255)
                                                                                                                        gcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACt                                                                             <  4:1450893/96‑1 (MQ=255)
                                                                                                                                gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                <  2:1306617/149‑1 (MQ=255)
                                                                                                                                gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                <  2:281714/149‑1 (MQ=255)
                                                                                                                                  gTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACt                                                                             >  2:558235/1‑86 (MQ=255)
                                                                                                                                  gTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACt                                                                             <  1:558235/86‑1 (MQ=255)
                                                                                                                                  gTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                  >  1:1361719/1‑129 (MQ=255)
                                                                                                                                   ttAAAACCCACGCGAAACGCGTGGAGGTCAAACTCCGCCTGt                                                                                                                        >  1:224399/1‑42 (MQ=255)
                                                                                                                                     aaaaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTc                 >  4:1263067/1‑143 (MQ=255)
                                                                                                                                        caccACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGc                    <  3:1140113/135‑1 (MQ=255)
                                                                                                                                         cccACGCGAAACGCGTGGTGGTCAAACTCCTCCTGCTTTATCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTGAGCAGCgtgt       <  4:950858/149‑1 (MQ=255)
                                                                                                                                         accACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCgtgt       >  1:349196/2‑149 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAg                                                     >  3:130980/1‑102 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATg                                             <  2:688252/110‑1 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATg                                             >  1:688252/1‑110 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCa                      >  1:564448/1‑133 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCa                      <  2:564448/133‑1 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagca             <  4:822912/142‑1 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTc      >  2:83654/1‑149 (MQ=255)
                                                                                                                                          ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCCAAGCTCACAAg                                                     <  4:130980/102‑1 (MQ=255)
                                                                                                                                           cACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                    <  1:1103098/118‑1 (MQ=255)
                                                                                                                                           cACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                    >  2:1103098/1‑118 (MQ=255)
                                                                                                                                            aCGCGAAACGCGTGGTGGTCAAACTCCGCTTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGTCATGGTGTTTCTCCt                                  <  2:1353336/119‑1 (MQ=255)
                                                                                                                                            aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGTCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                  >  1:1353336/1‑119 (MQ=255)
                                                                                                                                            aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTGTCCTTc                                >  3:879893/1‑121 (MQ=255)
                                                                                                                                             cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAg                                                  <  3:443329/102‑1 (MQ=255)
                                                                                                                                              gcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCAt                                              >  4:1229508/1‑105 (MQ=255)
                                                                                                                                               cgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc  <  2:931073/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACC  >  NC_000913/4315292‑4315582

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