Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,922 G→A 5.3% Q75* (CAA→TAA) ‡ higA ← antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,9220GA 5.3% 246.1 / 3.0 75Q75* (CAA→TAA) ‡higAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (34/37);  new base A (2/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

TTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTC  >  NC_000913/3233775‑3234070
                                                                                                                                                   |                                                                                                                                                    
ttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCGTATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGa                                                                                                                                                     <  2:510305/149‑1 (MQ=255)
ttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGa                                                                                                                                                     <  3:312929/149‑1 (MQ=255)
     gCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCAt                                                                                                                                                <  2:437159/149‑1 (MQ=255)
        tgtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcc                                                                                                                                             <  3:164946/149‑3 (MQ=255)
           tCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt                                                                                                                                          <  2:266535/149‑1 (MQ=255)
              aGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTaaa                                                                                                                                       <  2:118326/149‑1 (MQ=255)
                  ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                                                                   >  3:210720/1‑149 (MQ=255)
                  ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                                                                   <  3:103661/149‑1 (MQ=255)
                      tttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATCTTGATCCGTAATGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCaa                                                                                                                               >  2:426241/1‑149 (MQ=255)
                       ttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCaaa                                                                                                                              <  4:561448/149‑1 (MQ=255)
                           actTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCAGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCg                                                                                                                          >  2:505505/2‑149 (MQ=255)
                            ctTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCgg                                                                                                                         <  2:197470/149‑1 (MQ=255)
                            ctTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCgg                                                                                                                         <  1:81572/149‑1 (MQ=255)
                             ttCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg                                                                                                                        >  4:190105/1‑149 (MQ=255)
                             ttCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCAt                                                                                                                                                <  3:239877/125‑1 (MQ=255)
                             ttCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCAt                                                                                                                                                >  4:239877/1‑125 (MQ=255)
                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                                     <  3:392949/149‑1 (MQ=255)
                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                                     >  3:422335/1‑149 (MQ=255)
                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                                     <  1:302689/149‑1 (MQ=255)
                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                                     <  1:158349/149‑1 (MQ=255)
                                ccGCNCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                                                                     >  3:459075/1‑149 (MQ=255)
                                  gCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACTGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg                                                                                                                   >  2:360390/1‑149 (MQ=255)
                                   cTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa                                                                                                                  <  1:87869/149‑1 (MQ=255)
                                   cTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGa                                                                                                                  <  2:57711/149‑1 (MQ=255)
                                    tCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAt                                                                                                                 <  4:472999/149‑1 (MQ=255)
                                      aaaaCGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttc                                                                                                               >  2:514789/1‑149 (MQ=255)
                                      aaaaCGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttc                                                                                                               >  1:333021/1‑149 (MQ=255)
                                         aCGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttcttc                                                                                                            <  4:178524/149‑1 (MQ=255)
                                               gACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgc                                                                                                      >  3:216010/1‑149 (MQ=255)
                                               gACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgc                                                                                                      >  3:247908/1‑149 (MQ=255)
                                                acacCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg                                                                                                     >  2:184144/1‑149 (MQ=255)
                                                  acCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGt                                                                                                   >  3:29267/1‑149 (MQ=255)
                                                  acCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGt                                                                                                   >  3:17601/1‑149 (MQ=255)
                                                    cATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTa                                                                                                 >  1:585439/1‑149 (MQ=255)
                                                         aTTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTg                                                                                            >  3:84911/1‑149 (MQ=255)
                                                           ttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc                                                                                          <  1:314383/149‑1 (MQ=255)
                                                                 ccAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCacacac                                                                                    >  2:318495/1‑149 (MQ=255)
                                                                  cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACAcc                                                                                   >  3:565205/1‑149 (MQ=255)
                                                                  cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACAcc                                                                                   >  4:414953/1‑149 (MQ=255)
                                                                      ttCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGAt                                                                               <  2:353492/149‑1 (MQ=255)
                                                                          ggCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATTCCGCCAGGCATGGCTTGAGCCATGGTATTAAATTCCGCAAATTCGGGCGCTGATTCTTCACACGCGGTTATTTTGGCACACACCAGATCCAg                                                                           <  2:133379/149‑1 (MQ=255)
                                                                                tCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAggg                                                                     >  3:114927/1‑149 (MQ=255)
                                                                                tCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAggg                                                                     <  2:286630/149‑1 (MQ=255)
                                                                                   gAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGtt                                                                  >  2:533829/1‑149 (MQ=255)
                                                                                   gAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGtt                                                                  >  4:168484/1‑149 (MQ=255)
                                                                                         gTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAgg                                                            >  2:341844/1‑149 (MQ=255)
                                                                                          ttAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGa                                                           >  3:300883/1‑149 (MQ=255)
                                                                                               ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGtt                                                      <  1:342133/149‑1 (MQ=255)
                                                                                               ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGtt                                                      <  4:247908/149‑1 (MQ=255)
                                                                                                cATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTg                                                     >  1:353033/1‑149 (MQ=255)
                                                                                                 atatTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGa                                                    >  2:376998/1‑149 (MQ=255)
                                                                                                 atatTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGa                                                    <  3:190105/149‑1 (MQ=255)
                                                                                                    ttGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGagca                                                 <  4:17601/149‑1 (MQ=255)
                                                                                                         ccATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAga                                            <  3:105289/149‑1 (MQ=255)
                                                                                                          cATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAgat                                           >  1:201220/1‑149 (MQ=255)
                                                                                                           aTAAGGGAACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATCCCGCAAATTCGGGCGCTGAATCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTATCAGGCTCGTTGAACAGCAgatg                                          >  1:33459/1‑149 (MQ=255)
                                                                                                            tAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAgatga                                         <  1:376998/149‑1 (MQ=255)
                                                                                                             aaGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAgatgat                                        <  1:514789/149‑1 (MQ=255)
                                                                                                              aGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATc                                       <  3:168484/149‑1 (MQ=255)
                                                                                                                 gTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCTGGATCGTTGAGCAGCAGATGATCTAc                                    <  2:585439/149‑1 (MQ=255)
                                                                                                                 gTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTAc                                    >  2:480956/1‑149 (MQ=255)
                                                                                                                     gAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGt                                >  1:417834/1‑149 (MQ=255)
                                                                                                                        tCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTcc                             <  1:30480/149‑1 (MQ=255)
                                                                                                                         cACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCa                            <  4:29267/149‑1 (MQ=255)
                                                                                                                          aCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAg                           >  1:51207/1‑149 (MQ=255)
                                                                                                                          aCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAg                           >  3:115955/1‑149 (MQ=255)
                                                                                                                              cTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGCTTTCAGGTTCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGcc                       >  1:494419/1‑149 (MQ=255)
                                                                                                                               tataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCt                      <  1:505505/149‑1 (MQ=255)
                                                                                                                                   ccGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgt                  <  3:141621/149‑1 (MQ=255)
                                                                                                                                   ccGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGgt                  <  4:88251/149‑1 (MQ=255)
                                                                                                                                        tgccATGGCTCAAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTAtt             >  3:214396/4‑149 (MQ=255)
                                                                                                                                          gCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGt           <  1:360390/149‑1 (MQ=255)
                                                                                                                                            aTGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTc         <  4:114927/149‑1 (MQ=255)
                                                                                                                                               tcTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCAAACACCAGATCCATCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTTTATTGTTCCTc      <  2:33459/148‑1 (MQ=255)
                                                                                                                                                 ttGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAACAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGt    >  3:367575/1‑149 (MQ=255)
                                                                                                                                                   gAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAATTCCAGCGCCTGGGTGTATTGTTCCTCGTTc  <  4:216010/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAGCAGCAGATGATCTACCAGTTCCAGCGCCTGGGTGTATTGTTCCTCGTTC  >  NC_000913/3233775‑3234070

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