Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,439 A→T 5.9% F34I (TTT→ATT) ‡ phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,4390AT 5.9% 296.1 / 7.5 85F34I (TTT→ATT) ‡phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (36/44);  new base T (2/3);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

ACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGAC  >  NC_000913/4315291‑4315585
                                                                                                                                                    |                                                                                                                                                  
accaACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGaaa                                                                                                                                                    <  1:780405/149‑1 (MQ=255)
      tCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtg                                                                                                                                              >  3:705057/1‑149 (MQ=255)
       ccTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgg                                                                                                                                             <  4:646275/149‑1 (MQ=255)
         tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGTCGTTAAAACCCACGCGAAACGCgtggtg                                                                                                                                           >  1:143771/1‑149 (MQ=255)
         tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                           <  3:780296/149‑3 (MQ=255)
         tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgtggtg                                                                                                                                           >  1:662666/1‑149 (MQ=255)
            aTTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTc                                                                                                                                        <  1:462634/149‑1 (MQ=255)
                cGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCACAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAc                                                                                                                                    <  1:131403/149‑1 (MQ=255)
                   cGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTcc                                                                                                                                 >  3:52451/1‑149 (MQ=255)
                    gATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCg                                                                                                                                <  1:333888/149‑1 (MQ=255)
                       ccAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCCTGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCt                                                                                                                             <  3:499618/149‑1 (MQ=255)
                       ccAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCt                                                                                                                             <  2:659412/149‑1 (MQ=255)
                        cAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTg                                                                                                                            >  3:524746/1‑149 (MQ=255)
                         aGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGt                                                                                                                           >  2:106756/1‑149 (MQ=255)
                          gTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                           <  1:667999/132‑3 (MQ=255)
                          gTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                           >  2:667999/1‑130 (MQ=255)
                           ttGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGttt                                                                                                                         >  1:16843/1‑149 (MQ=255)
                            tGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtttt                                                                                                                        <  1:477232/149‑1 (MQ=255)
                               cATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                     <  2:143771/149‑1 (MQ=255)
                               cATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                     >  4:423702/1‑149 (MQ=255)
                                 tgatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTc                                                                                                                   <  3:637744/149‑1 (MQ=255)
                                 tgatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTc                                                                                                                   <  4:705263/149‑1 (MQ=255)
                                     gCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgg                                                                                                                                             >  3:554967/1‑119 (MQ=255)
                                     gCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaa                                                                                                               <  4:789589/149‑1 (MQ=255)
                                     gCAGATGAAACTGCAAATGCAGGCCGACCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaa                                                                                                               >  1:222146/1‑149 (MQ=255)
                                       aGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCGATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                     <  2:648465/141‑1 (MQ=255)
                                       aGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGc                                                                                                                     >  1:648465/1‑141 (MQ=255)
                                       aGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAAt                                                                                                             >  1:201748/1‑149 (MQ=255)
                                       aGATGAAACTGCAAATGCAGGCCGACCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAAt                                                                                                             <  2:222146/149‑1 (MQ=255)
                                        gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATc                                                                                                            >  1:629552/1‑149 (MQ=255)
                                         aTGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCa                                                                                                           <  4:128898/149‑1 (MQ=255)
                                         aTGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCa                                                                                                           >  2:400798/1‑149 (MQ=255)
                                            aaaCTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgcg                                                                                                        >  2:553644/1‑149 (MQ=255)
                                                   aaaTGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAc                                                                                                 >  2:552514/1‑149 (MQ=255)
                                                         aGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGc                                                                                                                               >  3:492778/1‑113 (MQ=255)
                                                         aGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGc                                                                                                                               <  4:492778/113‑1 (MQ=255)
                                                             cGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGgtg                                                                                       >  1:608418/1‑149 (MQ=255)
                                                                            cTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGcc                                                                        >  2:639887/1‑149 (MQ=255)
                                                                              tCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCgg                                                                      <  3:423702/149‑1 (MQ=255)
                                                                                  cATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                  >  2:406708/1‑149 (MQ=255)
                                                                                  cATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                  <  2:201748/149‑1 (MQ=255)
                                                                                       agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTc                                                             >  1:278019/1‑149 (MQ=255)
                                                                                       agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTc                                                             <  2:278019/149‑1 (MQ=255)
                                                                                           gcgcCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACaa                                                         >  2:254393/1‑149 (MQ=255)
                                                                                              ccAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCa                                                      >  2:327686/1‑149 (MQ=255)
                                                                                              ccAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCa                                                      >  1:244781/1‑149 (MQ=255)
                                                                                              ccAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCa                                                      >  2:400894/1‑149 (MQ=255)
                                                                                               cAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAg                                                     <  2:32719/149‑1 (MQ=255)
                                                                                                       aGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtg                                             <  1:406708/149‑1 (MQ=255)
                                                                                                        gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtgt                                            <  4:705057/149‑1 (MQ=255)
                                                                                                        gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtgt                                            >  1:103476/1‑149 (MQ=255)
                                                                                                        gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtgt                                            >  3:693709/1‑149 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                     <  2:60855/106‑1 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                     >  1:60855/1‑106 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgt                                              >  1:132484/1‑145 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgt                                              <  2:132484/145‑1 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGttt                                          >  2:57142/1‑149 (MQ=255)
                                                                                                          gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGttt                                          <  2:244781/149‑1 (MQ=255)
                                                                                                           cATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                         >  4:635851/1‑149 (MQ=255)
                                                                                                             tGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                       >  3:397632/1‑149 (MQ=255)
                                                                                                             tGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCCGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                       <  4:449937/149‑1 (MQ=255)
                                                                                                              gTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTcc                                      <  1:677821/149‑1 (MQ=255)
                                                                                                                 tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                   <  2:608418/149‑1 (MQ=255)
                                                                                                                 tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                   <  1:57142/149‑1 (MQ=255)
                                                                                                                 tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                   <  4:258019/149‑1 (MQ=255)
                                                                                                                 tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                   <  1:254393/149‑1 (MQ=255)
                                                                                                                    gTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGa                                <  4:524746/149‑1 (MQ=255)
                                                                                                                     tCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGAt                               <  2:440489/149‑1 (MQ=255)
                                                                                                                      cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATg                              <  4:396793/149‑1 (MQ=255)
                                                                                                                       gcgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATgg                             >  1:674163/1‑149 (MQ=255)
                                                                                                                           aGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCa                         <  3:117115/149‑1 (MQ=255)
                                                                                                                             aTTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGc                       <  1:552514/149‑1 (MQ=255)
                                                                                                                             aTTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGc                       <  1:400798/149‑1 (MQ=255)
                                                                                                                               tGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGt                     <  4:52451/149‑1 (MQ=255)
                                                                                                                                 gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                   <  2:358928/149‑1 (MQ=255)
                                                                                                                                 gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                   >  2:651808/1‑149 (MQ=255)
                                                                                                                                 gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                   >  2:179723/1‑149 (MQ=255)
                                                                                                                                 gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTca                   <  2:691746/149‑1 (MQ=255)
                                                                                                                                   gTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagc                 >  4:334440/1‑149 (MQ=255)
                                                                                                                                     tAAAACCCACGCGAAACGCGTGGTGGTCACACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagcag               <  4:253001/149‑1 (MQ=255)
                                                                                                                                      aaaaCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCagcagc              <  2:16843/149‑1 (MQ=255)
                                                                                                                                           ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTc         <  2:629552/149‑1 (MQ=255)
                                                                                                                                             aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTc                    <  4:258063/136‑1 (MQ=255)
                                                                                                                                             aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTc                    >  3:258063/1‑136 (MQ=255)
                                                                                                                                             aCGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcga       <  1:106756/149‑1 (MQ=255)
                                                                                                                                              cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgac      <  3:334440/149‑1 (MQ=255)
                                                                                                                                                  aaaCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgaccgac  <  1:639887/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGAC  >  NC_000913/4315291‑4315585

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