Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,449,050 T→A 11.9% intergenic (+32/+26) feaB → / ← tynA phenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,449,0500TA11.9% 207.4 / 17.9 67intergenic (+32/+26)feaB/tynAphenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base A (5/3);  total (38/29)
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.67e-01

GCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACATGTGT  >  NC_000913/1448911‑1449198
                                                                                                                                           |                                                                                                                                                    
gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtt                                                                                                                                             >  1:239950/1‑149 (MQ=255)
 cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTATCGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGCGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACTTGGTGCCGttt                                                                                                                                            <  4:53558/149‑1 (MQ=255)
  tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttt                                                                                                                                           <  1:115408/149‑1 (MQ=255)
  tAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttt                                                                                                                                           <  4:391536/149‑1 (MQ=255)
   aaCTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                                                                                                                                          <  3:55445/149‑1 (MQ=255)
     cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttttt                                                                                                                                        <  1:192016/149‑1 (MQ=255)
     cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttttta                                                                                                                                        <  3:287926/149‑2 (MQ=255)
      ttACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                                                                       >  1:95391/1‑149 (MQ=255)
                 tggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACCGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAAt                                                                                                                            >  2:132263/1‑149 (MQ=255)
                   ggATGCAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGATCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa                                                                                                                          >  2:584945/1‑149 (MQ=255)
                   ggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa                                                                                                                          >  2:315618/1‑149 (MQ=255)
                    gATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtt                                                                                                                                             <  4:176518/129‑1 (MQ=255)
                    gATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtt                                                                                                                                             >  3:176518/1‑129 (MQ=255)
                         aGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc                                                                                                                    >  3:221144/1‑149 (MQ=255)
                               cAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctttctt                                                                                                              <  4:215087/149‑1 (MQ=255)
                                   aaCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgc                                                                                                          >  4:103630/1‑149 (MQ=255)
                                        gCCGTGATTTTGGCCCCGACTGGCTGGGCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCa                                                                                                            >  3:519305/1‑142 (MQ=255)
                                        gCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCa                                                                                                            <  4:519305/142‑1 (MQ=255)
                                          cGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTa                                                                                                   >  3:187901/1‑149 (MQ=255)
                                                tttGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTg                                                                                             <  1:315618/149‑1 (MQ=255)
                                                     ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGtt                                                                                        >  2:316918/1‑149 (MQ=255)
                                                      cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGttt                                                                                       >  4:59413/1‑149 (MQ=255)
                                                         gACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                              <  4:365081/139‑1 (MQ=255)
                                                         gACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                                                                              >  3:365081/1‑139 (MQ=255)
                                                         gACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGt                                                                                    <  1:114057/149‑1 (MQ=255)
                                                             ggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                                                                                <  1:330805/149‑1 (MQ=255)
                                                             ggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaaa                                                                                <  3:280871/149‑1 (MQ=255)
                                                                tggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaa                                                                             <  2:296240/149‑1 (MQ=255)
                                                                 ggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                                                                            <  4:221144/149‑1 (MQ=255)
                                                                  gACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                                                                           <  4:498095/149‑1 (MQ=255)
                                                                   aCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGtt                                                                          <  2:396644/149‑1 (MQ=255)
                                                                      gTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                                                                       <  3:103630/149‑1 (MQ=255)
                                                                      gTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                                                                       >  2:274437/1‑149 (MQ=255)
                                                                             tgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc                                                                >  4:248928/1‑149 (MQ=255)
                                                                                 aCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca                                                            <  2:45367/149‑1 (MQ=255)
                                                                                      gTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa                                                       >  4:527757/1‑149 (MQ=255)
                                                                                       tCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                                                                            >  4:463590/1‑127 (MQ=255)
                                                                                       tCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaagaag                                                                            <  3:463590/127‑1 (MQ=255)
                                                                                       tCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTAt                                                      >  2:77713/1‑149 (MQ=255)
                                                                                          gtgtgtgtACCGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTa                                                   <  1:132263/149‑1 (MQ=255)
                                                                                           tgtgtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTAc                                                  <  3:527757/149‑1 (MQ=255)
                                                                                                 tACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCaga                                                          >  2:257998/1‑135 (MQ=255)
                                                                                                 tACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCaga                                                          <  1:257998/135‑1 (MQ=255)
                                                                                                  aCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATtcg                                           <  3:203994/149‑1 (MQ=255)
                                                                                                        ttAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGc                                     >  4:360461/1‑149 (MQ=255)
                                                                                                           aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTAt                                                      <  4:564814/129‑1 (MQ=255)
                                                                                                           aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTAt                                                      >  3:564814/1‑129 (MQ=255)
                                                                                                           aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATa                                  >  4:583758/1‑149 (MQ=255)
                                                                                                           aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATa                                  >  2:423751/1‑149 (MQ=255)
                                                                                                                gTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCgg                             <  4:191424/149‑1 (MQ=255)
                                                                                                                gTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCgg                             >  3:463746/1‑149 (MQ=255)
                                                                                                                 ttCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGc                            <  1:367229/149‑1 (MQ=255)
                                                                                                                       cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCActct                      >  4:322341/1‑149 (MQ=255)
                                                                                                                        acaaGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCtt                     >  1:330998/3‑149 (MQ=255)
                                                                                                                         tAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTc                    >  3:294502/1‑149 (MQ=255)
                                                                                                                           aGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCgg                  >  1:191643/1‑149 (MQ=255)
                                                                                                                            gTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgc                 >  3:466328/1‑149 (MQ=255)
                                                                                                                             tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcg                >  2:85567/1‑149 (MQ=255)
                                                                                                                             tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcg                >  1:567071/1‑149 (MQ=255)
                                                                                                                             tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcg                >  3:99219/1‑149 (MQ=255)
                                                                                                                              aaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGgcgc               >  3:222223/1‑149 (MQ=255)
                                                                                                                                aaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg             >  2:29577/1‑149 (MQ=255)
                                                                                                                                aaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg             <  3:59413/149‑1 (MQ=255)
                                                                                                                                aaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCgg             >  4:291283/1‑149 (MQ=255)
                                                                                                                                 aaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCggg            >  3:272718/1‑149 (MQ=255)
                                                                                                                                  aaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCAGGc           >  1:568531/1‑149 (MQ=255)
                                                                                                                                        aCCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtg     >  1:375955/1‑149 (MQ=255)
                                                                                                                                           tGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtgtgt  >  3:310508/1‑149 (MQ=255)
                                                                                                                                           tGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACAtgtgt  <  1:584945/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTATGTACCCATTCGGTCGGCATAATCGGCCACTCTTCGGCGCGGGCCACATGTGT  >  NC_000913/1448911‑1449198

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