Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,793,026 G→T 8.7% A178E (GCG→GAG)  btuD ← vitamin B12 ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,793,0260GT 8.7% 850.7 / 52.4 242A178E (GCG→GAG) btuDvitamin B12 ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (102/119);  new base T (14/7);  total (116/126)
Fisher's exact test for biased strand distribution p-value = 1.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 248 total aligned reads displayed.

GCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCC  >  NC_000913/1792877‑1793174
                                                                                                                                                     |                                                                                                                                                    
    cAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCt                                                                                                                                                  <  1:4273147/150‑1 (MQ=255)
           ggcggcGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAAtt                                                                                                                                           <  1:924275/150‑1 (MQ=255)
            gcggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGTGCATGACGCAATGTGTGGTTGAGATCGAGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtat                                                                                                                                          >  1:2538299/1‑148 (MQ=255)
            gcggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAttt                                                                                                                                          >  2:2380660/1‑150 (MQ=255)
            gcggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAttt                                                                                                                                          >  1:2118445/1‑150 (MQ=255)
             cggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtttt                                                                                                                                         <  1:2428902/150‑1 (MQ=255)
             cggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtttt                                                                                                                                         <  1:1589341/150‑1 (MQ=255)
             cggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtttt                                                                                                                                         <  1:240107/150‑1 (MQ=255)
             cggcgTGAGCACCTCTTCCCAGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAAtttt                                                                                                                                         <  1:2735813/150‑1 (MQ=255)
              ggcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTg                                                                                                                                        <  1:2517770/150‑1 (MQ=255)
               gcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGt                                                                                                                                       >  2:1592943/1‑150 (MQ=255)
               gcgTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGt                                                                                                                                       >  1:4194935/1‑150 (MQ=255)
                  tGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTa                                                                                                                                    >  2:2116631/1‑150 (MQ=255)
                   gAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTaa                                                                                                                                   <  2:2875969/150‑1 (MQ=255)
                     gCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACg                                                                                                                                 >  2:2622899/1‑150 (MQ=255)
                      cACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGc                                                                                                                                >  1:3960299/1‑150 (MQ=255)
                         ctctTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACt                                                                                                                             >  1:1388892/1‑150 (MQ=255)
                            ttCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttg                                                                                                                          <  1:3856644/150‑1 (MQ=255)
                             tCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttgt                                                                                                                         >  1:3251360/1‑150 (MQ=255)
                              cccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttgtt                                                                                                                        <  2:2219513/150‑1 (MQ=255)
                              cccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTAGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttgtt                                                                                                                        <  2:3554741/150‑1 (MQ=255)
                               ccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCTTGACTGCTAATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATATTGTCTAACGCACTttgttg                                                                                                                       >  1:3199727/1‑150 (MQ=255)
                               ccTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTttgttg                                                                                                                       <  2:1919865/150‑1 (MQ=255)
                                cTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTgc                                                                                                                      <  1:2921006/150‑1 (MQ=255)
                                cTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTgc                                                                                                                      >  2:2364921/1‑150 (MQ=255)
                                 tGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGTCGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTgcg                                                                                                                     <  2:2281487/150‑1 (MQ=255)
                                    gTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCaa                                                                                                                  >  2:1029368/1‑150 (MQ=255)
                                    gTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCaa                                                                                                                  >  1:3412523/1‑150 (MQ=255)
                                    gTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCACaa                                                                                                                  >  2:2960230/1‑150 (MQ=255)
                                      ccACTGGCCAGCATTTTTCCAACTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACa                                                                                                                <  2:2705270/150‑1 (MQ=255)
                                          tGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCaa                                                                                                            <  1:3528199/150‑1 (MQ=255)
                                             ccAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGAc                                                                                                         >  1:31310/1‑150 (MQ=255)
                                               aGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTg                                                                                                       >  2:417109/1‑150 (MQ=255)
                                               aGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTg                                                                                                       >  2:2460225/1‑150 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             >  2:1701958/1‑150 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             >  2:2928416/1‑150 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             >  1:1852087/1‑150 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             >  1:3561403/1‑150 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             <  2:4071897/150‑1 (MQ=255)
                                                         cACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCt                                                                                             >  1:3165988/1‑150 (MQ=255)
                                                          aCCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCtc                                                                                            >  1:3695595/1‑150 (MQ=255)
                                                             ttttAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCtcatc                                                                                         >  2:69175/1‑150 (MQ=255)
                                                              tttAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCtcatca                                                                                        <  1:1560579/150‑1 (MQ=255)
                                                                tAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGCCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCaag                                                                                      <  2:2109465/150‑1 (MQ=255)
                                                                 aGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCaaga                                                                                     <  1:1941475/150‑1 (MQ=255)
                                                                 aGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCaaga                                                                                     <  1:4063155/150‑1 (MQ=255)
                                                                 aGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATATTGTCTAACGCACTTTGTTGCGCACCATCAAGACTGTTCATCGGCTCATCaaga                                                                                     >  1:2271834/1‑150 (MQ=255)
                                                                    aaCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagc                                                                                  >  2:77133/1‑150 (MQ=255)
                                                                    aaCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagc                                                                                  <  2:4154791/150‑1 (MQ=255)
                                                                    aaCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagc                                                                                  >  1:807391/1‑150 (MQ=255)
                                                                    aaCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTGGCTGACCCAGCGCGCTCAGAATTTTGTCTAACGCCCTTTGGTGCGCAACATCAAGACTGTTCAGCGGCTCATCAAGAagc                                                                                  >  2:2126846/1‑150 (MQ=255)
                                                                     aCCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagca                                                                                 >  1:3031999/1‑150 (MQ=255)
                                                                      ccACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAagcag                                                                                <  2:3051192/150‑1 (MQ=255)
                                                                       cacgcacgATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACAAAGCGCGCTCAGAATTTTGTCTAACGCACATTGTTGCGCAACATCAACACTGTGCATATGCTCATCAAGAagcacc                                                                               >  2:1946404/1‑148 (MQ=255)
                                                                            ccGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTg                                                                          <  1:2093281/150‑1 (MQ=255)
                                                                                  cgcATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGc                                                                    <  2:1768137/150‑1 (MQ=255)
                                                                                   gcATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCg                                                                   >  2:4089763/1‑150 (MQ=255)
                                                                                    cATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCgg                                                                  >  1:1089563/1‑150 (MQ=255)
                                                                                     aTGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCggg                                                                 >  2:2040670/1‑150 (MQ=255)
                                                                                         cGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTGCTTGCTGACACAGCGCGATCAGAATTTGGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTg                                                             <  2:553601/150‑1 (MQ=255)
                                                                                         cGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTg                                                             <  1:1201289/150‑1 (MQ=255)
                                                                                           cAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGc                                                           <  2:1787965/150‑1 (MQ=255)
                                                                                            aaTGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCt                                                          >  2:1770713/1‑150 (MQ=255)
                                                                                             aTGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCtt                                                         <  1:659294/150‑1 (MQ=255)
                                                                                              tgtgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtg                                                        <  1:2172597/150‑1 (MQ=255)
                                                                                               gtgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtgt                                                       <  2:3687252/150‑1 (MQ=255)
                                                                                               gtgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCCCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtgt                                                       <  1:3481818/150‑1 (MQ=255)
                                                                                                tgtgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCTCTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTtgtg                                                      <  2:224656/150‑1 (MQ=255)
                                                                                                  tgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTgtggt                                                    <  2:3236508/150‑1 (MQ=255)
                                                                                                  tgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGGTCTGAATTTTGTCTAACGCTGTTTGTTGCGCTACTTGATGAGTGTTCTTCGGGTGTTGAAGAAGCTGCTATTGGCGTGGGGGATTGGGTTgtggt                                                    >  2:437416/1‑150 (MQ=255)
                                                                                                  tgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTgtggt                                                    >  2:2320689/1‑150 (MQ=255)
                                                                                                  tgGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTgtggt                                                    >  1:3344712/1‑150 (MQ=255)
                                                                                                    gTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGgtgt                                                  >  1:2175807/1‑150 (MQ=255)
                                                                                                     ttGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGtgtg                                                 >  2:773095/1‑150 (MQ=255)
                                                                                                       gagaTCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGAt                                               >  2:4036590/1‑150 (MQ=255)
                                                                                                         gaTCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGAttt                                             >  2:748707/1‑150 (MQ=255)
                                                                                                          aTCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTg                                            <  2:2773943/150‑1 (MQ=255)
                                                                                                          aTCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTg                                            >  1:2763775/1‑150 (MQ=255)
                                                                                                          aTCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTg                                            <  2:1388892/150‑1 (MQ=255)
                                                                                                           tCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGc                                           >  2:386931/1‑150 (MQ=255)
                                                                                                           tCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGc                                           <  1:3899259/150‑1 (MQ=255)
                                                                                                           tCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGc                                           <  2:3318456/150‑1 (MQ=255)
                                                                                                            cGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCa                                          <  2:826750/150‑1 (MQ=255)
                                                                                                                aCTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAacac                                      <  1:4121493/150‑1 (MQ=255)
                                                                                                                     tcatcaCAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCg                                 >  1:2354272/1‑150 (MQ=255)
                                                                                                                        tcaCAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCACCACCACCgcag                              >  1:742529/1‑150 (MQ=255)
                                                                                                                         cacaATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTCCTCAACAACAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATAGGCTTGTGGTGTGATTTGCAACACCACCgcagc                             >  2:4227395/1‑150 (MQ=255)
                                                                                                                              tCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGAc                        >  2:2808895/1‑150 (MQ=255)
                                                                                                                              tCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGAc                        <  1:4036590/150‑1 (MQ=255)
                                                                                                                                 ccAGTCCTTGCTGACACAGCTCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTa                     >  1:438878/1‑150 (MQ=255)
                                                                                                                                      ccTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGt                <  2:3431089/150‑1 (MQ=255)
                                                                                                                                       cTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGtt               >  2:1575115/1‑150 (MQ=255)
                                                                                                                                        ttGCTGACACAGCGCGCTCAGACTTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTg              <  1:1289406/150‑1 (MQ=255)
                                                                                                                                        ttGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTg              <  1:26604/150‑1 (MQ=255)
                                                                                                                                        ttGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTg              <  1:69175/150‑1 (MQ=255)
                                                                                                                                          gCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGcc            >  2:3756013/1‑150 (MQ=255)
                                                                                                                                              acacAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTc        >  1:1059088/1‑150 (MQ=255)
                                                                                                                                               cacaGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCa       <  1:1365080/150‑1 (MQ=255)
                                                                                                                                                acaGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCAc      >  1:2634421/1‑150 (MQ=255)
                                                                                                                                                   gcgcgcTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCAccgc   <  1:1756846/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACGCCCGATGCGCATGACGCAATGTGTGGTTGAGATCGTGACTGCTCATCACAATCGCCAGTCCTTGCTGACACAGCGCGCTCAGAATTTTGTCTAACGCACTTTGTTGCGCAACATCAAGACTGTTCATCGGCTCATCAAGAAGCAGCAATTGGCCTGCGGGATTGGCTTGTGGTGTGATTTGCAACACCACCGCAGCAAGACGTACGCGTTGCCATTCACCGCC  >  NC_000913/1792877‑1793174

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