Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,555,741 C→A 5.6% A886E (GCG→GAG)  malT → mal regulon transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,555,7410CA 5.6% 293.0 / 3.6 86A886E (GCG→GAG) malTmal regulon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (42/39);  new base A (2/3);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTGAATGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAAC  >  NC_000913/3555597‑3555876
                                                                                                                                                |                                                                                                                                       
aaCGTGAATGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCgg                                                                                                                                       >  7:277073/1‑147 (MQ=255)
     gAATGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTacaa                                                                                                                                  >  2:299379/1‑147 (MQ=255)
       aTGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAacac                                                                                                                               >  5:277391/1‑148 (MQ=255)
        tGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTacaacaag                                                                                                                              >  8:135949/1‑146 (MQ=255)
        tGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACg                                                                                                                              >  8:181094/1‑148 (MQ=255)
         ggCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTacaac                                                                                                                                 >  7:220442/1‑144 (MQ=255)
         ggCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTacaac                                                                                                                                 <  8:220442/144‑1 (MQ=255)
          gCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGcc                                                                                                                            <  1:203650/148‑1 (MQ=255)
          gCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGc                                                                                                                             <  1:44514/147‑1 (MQ=255)
          gCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGc                                                                                                                             <  6:277391/147‑1 (MQ=255)
            aGGTACTGGGGCTGATCTACTCTGGTTACAGAAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGAATATCCGCAATCTGTATAAGAAACTCGGCGTGTCCCATCGCCACTATCAGGTACAACACGcc                                                                                                                            >  3:286674/1‑146 (MQ=255)
               tACTGGGGCTGATCTACTCTGGTTACAGCAATGAGGAAATTGCCGGAGAACTGGAAGTCGCGGCAACAACGATCAAAACGCATATCCGCAATCTGTATCAGAAACTGGGCGTGGCCCATCGCCTGGATGCGGTACAACACGCCCAGCa                                                                                                                       <  6:149576/148‑1 (MQ=255)
                  tGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAAtt                                                                                                                    >  3:42853/1‑148 (MQ=255)
                   ggggCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGccc                                                                                                                           >  6:155189/1‑140 (MQ=255)
                   ggggCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGccc                                                                                                                           <  5:155189/140‑1 (MQ=255)
                   ggggCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGAGGTACAACACGCCCAGCAATTg                                                                                                                   >  7:273692/1‑148 (MQ=255)
                      gCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCAGCATCAAAACGCATATCCGCAATCGGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCt                                                                                                                 >  4:105595/1‑147 (MQ=255)
                       cTGATATACTCTGGTTACAGAAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCAGATCCGCAATCTATAGAAGAAACTCGGCGTGGCCCAGCTCCAGGATGCGGTACACCACGCCAAGCAAGTgagg                                                                                                                >  6:164724/1‑144 (MQ=255)
                          aTCTACTCTTGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGTTACAACACGCCCAGCAATTGCTg                                                                                                                >  2:49877/1‑144 (MQ=255)
                          aTCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTg                                                                                                                >  8:204154/1‑144 (MQ=255)
                          aTCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAg                                                                                                             <  8:277073/147‑1 (MQ=255)
                          aTCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAg                                                                                                             <  1:299379/147‑1 (MQ=255)
                           tCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgat                                                                                                           >  8:121122/1‑148 (MQ=255)
                              aCTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        >  5:273021/1‑148 (MQ=255)
                               ctctGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        <  7:135949/147‑1 (MQ=255)
                               ctctGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          <  4:42853/145‑1 (MQ=255)
                               ctctGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          <  7:181094/145‑1 (MQ=255)
                               ctctGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGatgatg                                                                                                       <  3:235829/148‑1 (MQ=255)
                                   ggTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        <  4:77371/143‑1 (MQ=255)
                                   ggTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        >  3:77371/1‑143 (MQ=255)
                                   ggTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATgggg                                                                                                    >  4:273449/1‑147 (MQ=255)
                                    gTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATggg                                                                                                     <  6:273021/145‑1 (MQ=255)
                                        cAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGc                                                                                              >  1:202727/1‑148 (MQ=255)
                                          gCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCgt                                                                                            >  5:71398/1‑148 (MQ=255)
                                             aTGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTa                                                                                         <  1:160802/148‑1 (MQ=255)
                                                     aTTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        >  6:277679/1‑125 (MQ=255)
                                                     aTTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatgat                                                                                                        <  5:277679/125‑1 (MQ=255)
                                                     aTTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGc                                                                                 >  8:192142/1‑148 (MQ=255)
                                                         ccGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTAGCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCgg                                                                              <  5:138006/147‑1 (MQ=255)
                                                           ggAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATa                                                                           <  1:19100/148‑1 (MQ=255)
                                                                 cTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACgcgc                                                                     <  7:121122/148‑1 (MQ=255)
                                                                  tGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACgcgc                                                                     >  6:232939/1‑147 (MQ=255)
                                                                           ctgCAACCACCATCAAAACGCATATCCGCAATCTGTATCATAAACTCGGCGTGGCCCATCGCCAGGATGAGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGc                                                           <  7:70784/146‑1 (MQ=255)
                                                                           cgGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGc                                                           <  5:154630/148‑1 (MQ=255)
                                                                               aaCCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTAc                                                       <  3:105595/148‑1 (MQ=255)
                                                                                  caccaTCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGACCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTAGAtc                                                    <  4:286674/148‑1 (MQ=255)
                                                                                   accaTCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACAtct                                                   <  2:202727/148‑1 (MQ=255)
                                                                                     caTCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          >  4:55928/1‑91 (MQ=255)
                                                                                     caTCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          <  3:55928/91‑1 (MQ=255)
                                                                                      aTCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCATGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTg                                                >  2:211346/1‑148 (MQ=255)
                                                                                       tCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTAc                                                       <  5:89853/140‑1 (MQ=255)
                                                                                       tCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTAc                                                       >  6:89853/1‑140 (MQ=255)
                                                                                        cAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGca                                              >  3:16263/1‑148 (MQ=255)
                                                                                              gCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGAGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTc                                        <  8:273692/148‑1 (MQ=255)
                                                                                               cATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCa                                       <  7:204154/148‑1 (MQ=255)
                                                                                                   tCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGc                                   <  4:16263/148‑1 (MQ=255)
                                                                                                    ccGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAAt                                     >  7:97992/1‑145 (MQ=255)
                                                                                                    ccGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAAt                                     <  8:97992/145‑1 (MQ=255)
                                                                                                     cGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCt                                  <  6:258082/147‑1 (MQ=255)
                                                                                                      gCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          >  3:229722/1‑74 (MQ=255)
                                                                                                      gCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          <  4:229722/74‑1 (MQ=255)
                                                                                                      gCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          >  7:200666/1‑74 (MQ=255)
                                                                                                      gCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAgatg                                                                                                          <  8:200666/74‑1 (MQ=255)
                                                                                                      gCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCa                                >  7:332876/1‑148 (MQ=255)
                                                                                                       cAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAg                                                                                                             >  6:66049/1‑70 (MQ=255)
                                                                                                       cAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAg                                                                                                             <  5:66049/70‑1 (MQ=255)
                                                                                                       cAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCAc                               >  2:281133/1‑148 (MQ=255)
                                                                                                             gTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGt                         >  2:308485/1‑148 (MQ=255)
                                                                                                              tATCAGAAACTCGTAGTGGCCAATCGCCAGGATGAGGTACAACACGCCCAGCAAGTGCTGAAGATTATTGGGTACTGCGTGTAAGTTTAGCCGGATAACGCGCCATATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGtt                        <  5:164724/148‑1 (MQ=255)
                                                                                                               aTCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGtga                       >  4:179059/1‑146 (MQ=255)
                                                                                                                 cAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACgcgc                                                                     >  1:59243/1‑100 (MQ=255)
                                                                                                                 cAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACgcgc                                                                     <  2:59243/100‑1 (MQ=255)
                                                                                                                  aGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGACAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGc                    >  1:230036/1‑148 (MQ=255)
                                                                                                                      aCTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATc                >  7:255816/1‑148 (MQ=255)
                                                                                                                       cTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCt               >  5:174224/1‑148 (MQ=255)
                                                                                                                         cGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGt             >  6:204735/1‑148 (MQ=255)
                                                                                                                           gCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGttt           <  5:204735/148‑1 (MQ=255)
                                                                                                                             gTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGTGGTACGGCATGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCt         <  1:211346/148‑1 (MQ=255)
                                                                                                                                  ccATCGCCAGGATGCGTTACAACACGCCAATCAATTTATCAAGATGATGGCCTACGGCGTGTAATTTTAGCCGGATAACTCGCCAGATCAGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCaa    <  1:49877/148‑1 (MQ=255)
                                                                                                                                    aTCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACAtctc                                                  <  5:63613/100‑1 (MQ=255)
                                                                                                                                    aTCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACAtctc                                                  >  6:63613/1‑100 (MQ=255)
                                                                                                                                    aTCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAAc  >  5:99764/1‑148 (MQ=255)
                                                                                                                                     tCGCGAGGATGCGGTACAACACGCCGAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTa                       <  6:315927/126‑1 (MQ=255)
                                                                                                                                      cGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACa                                                      >  2:176185/1‑94 (MQ=255)
                                                                                                                                      cGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACa                                                      <  1:176185/94‑1 (MQ=255)
                                                                                                                                      cGCCAGGATGCGGTACAACACGCCCAGCAATTGCAGAAGATGATGCGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAAAGCTCACACGCCTTAc                      >  5:315927/1‑126 (MQ=255)
                                                                                                                                                |                                                                                                                                       
AACGTGAATGGCAGGTACTGGGGCTGATCTACTCTGGTTACAGCAATGAGCAAATTGCCGGAGAACTGGAAGTCGCGGCAACCACCATCAAAACGCATATCCGCAATCTGTATCAGAAACTCGGCGTGGCCCATCGCCAGGATGCGGTACAACACGCCCAGCAATTGCTGAAGATGATGGGGTACGGCGTGTAAGTTTAGCCGGATAACGCGCCAGATCCGGCTTACATCTCTGCATCATTCAATGCTCACCCGCGTTACGCCATCTGTTTCTATCAAAC  >  NC_000913/3555597‑3555876

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