Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 336,866 G→C 5.9% G30A (GGC→GCC)  yahF → putative NAD(P)‑binding succinyl‑CoA synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913336,8660GC 5.9% 290.4 / 8.6 85G30A (GGC→GCC) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base C (3/2);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

ACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAAC  >  NC_000913/336721‑337008
                                                                                                                                                 |                                                                                                                                              
aCATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGt                                                                                                                                             <  3:12402/149‑1 (MQ=255)
 cATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTc                                                                                                                                            <  3:147012/149‑1 (MQ=255)
     cTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgt                                                                                                                                        <  2:12563/149‑2 (MQ=255)
      tGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCa                                                                                                                                       >  4:178907/1‑149 (MQ=255)
       gTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      <  1:136702/149‑1 (MQ=255)
       gTATGGCATGGCCGATGGATGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      >  4:45884/1‑149 (MQ=255)
         atggcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGc                                                                                                                                    >  4:95283/1‑149 (MQ=255)
         atggcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGc                                                                                                                                    <  3:126426/149‑1 (MQ=255)
             catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAttt                                                                                                                                >  4:300687/1‑149 (MQ=255)
                ggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTtgtg                                                                                                                             <  2:214934/149‑1 (MQ=255)
                 gcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtgg                                                                                                                            >  2:668722/1‑149 (MQ=255)
                 gcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtgg                                                                                                                            <  1:278055/149‑1 (MQ=255)
                   cGATGGTTGTGCGCTTAGCCAAACCTACGTAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCCAATAAACTCGACGGCGTGGAGCAGGCATTtgtgcga                                                                                                                          <  4:354630/149‑4 (MQ=255)
                    gATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcga                                                                                                                         <  2:390069/149‑1 (MQ=255)
                      tGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatg                                                                                                                       >  3:209595/1‑149 (MQ=255)
                       ggTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgatgg                                                                                                                      >  1:510009/1‑149 (MQ=255)
                            tgCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGAcc                                                                                                                 <  3:45884/149‑1 (MQ=255)
                                cTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGaaa                                                                                                             <  3:178907/149‑1 (MQ=255)
                                 tGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAAt                                                                                                            >  3:92228/1‑149 (MQ=255)
                                   aGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                          <  4:607775/149‑1 (MQ=255)
                                    gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGaa                                                                                                         <  1:70397/149‑1 (MQ=255)
                                    gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGaa                                                                                                         <  4:104756/149‑1 (MQ=255)
                                    gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGaa                                                                                                         <  4:700060/149‑1 (MQ=255)
                                       aaaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaa                                                                                                      >  1:173478/1‑149 (MQ=255)
                                       aaaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgttta                                                                                                                                     >  3:415211/1‑114 (MQ=255)
                                       aaaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgttta                                                                                                                                     <  4:415211/118‑5 (MQ=255)
                                         aCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                    <  4:209595/149‑1 (MQ=255)
                                         aCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                    >  2:768191/1‑149 (MQ=255)
                                          ccTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGTAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                   >  3:526221/1‑149 (MQ=255)
                                          ccTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAgg                                                                                                   <  3:255846/149‑1 (MQ=255)
                                            tACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      >  2:493341/1‑112 (MQ=255)
                                            tACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      <  3:697543/112‑1 (MQ=255)
                                            tACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      <  1:493341/112‑1 (MQ=255)
                                            tACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAg                                                                                                                                      >  4:697543/1‑112 (MQ=255)
                                             aCGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGt                                                                                                >  1:545675/1‑149 (MQ=255)
                                               ggAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATATAGGCGTGc                                                                                              <  4:66442/149‑1 (MQ=255)
                                                 aGGGTATTATGTAAGTTAAAATAGTCATTAAACCGAATACCTATTGTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGAAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGACGTGCTg                                                                                            >  4:529123/1‑149 (MQ=255)
                                                   ggAATTATGTCAGTTAAAATAGTCAGTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTgaa                                                                                          <  1:342232/149‑1 (MQ=255)
                                                       ttATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                                      <  2:705496/149‑1 (MQ=255)
                                                        tATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAtt                                                                                     <  1:635573/149‑1 (MQ=255)
                                                          tGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                   >  4:530269/1‑149 (MQ=255)
                                                            tCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                          <  4:511892/124‑1 (MQ=255)
                                                            tCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGa                                                                                                          >  3:511892/1‑124 (MQ=255)
                                                             cAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAACCTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGa                                                                                <  3:426786/149‑1 (MQ=255)
                                                               gTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGAct                                                                              >  4:348066/1‑149 (MQ=255)
                                                               gTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGAct                                                                              >  3:667741/1‑149 (MQ=255)
                                                                       aGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGc                                                                      <  3:665639/149‑1 (MQ=255)
                                                                          cATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCgg                                                                   <  1:381113/149‑1 (MQ=255)
                                                                           aTTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGa                                                                  <  4:108187/149‑1 (MQ=255)
                                                                            ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGACGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAg                                                                 >  4:349489/1‑149 (MQ=255)
                                                                               aaCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTg                                                              <  4:526221/149‑1 (MQ=255)
                                                                               aaCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCAGGAGCTg                                                              <  2:415207/149‑1 (MQ=255)
                                                                                aCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTgg                                                             <  4:659409/149‑1 (MQ=255)
                                                                                  cGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAg                                                           >  4:702990/1‑149 (MQ=255)
                                                                                         cTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAGTTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGa                                                    >  4:495867/1‑149 (MQ=255)
                                                                                         cTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGa                                                    >  1:120696/1‑149 (MQ=255)
                                                                                            ttttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGaaaa                                                 >  2:396377/1‑149 (MQ=255)
                                                                                            ttttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGaaaa                                                 >  2:714058/1‑149 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCg                                           >  2:4317/1‑149 (MQ=255)
                                                                                                         tcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACTGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtgat                                    >  2:626245/1‑149 (MQ=255)
                                                                                                            cTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTgatgat                                 <  1:396377/149‑1 (MQ=255)
                                                                                                             tGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGAtt                                <  4:123176/149‑1 (MQ=255)
                                                                                                             tGAAGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGTTGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAATTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAATACGGCTACCTGATGAtt                                >  4:266366/1‑149 (MQ=255)
                                                                                                               aTGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGt                              <  4:745853/149‑1 (MQ=255)
                                                                                                                 gTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATCGtca                            <  2:702598/149‑1 (MQ=255)
                                                                                                                  tCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGtcat                           >  3:423738/1‑149 (MQ=255)
                                                                                                                   cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGtcatc                          <  3:702990/149‑1 (MQ=255)
                                                                                                                    tATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGtcatca                         >  3:449199/1‑149 (MQ=255)
                                                                                                                     aTCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCaa                        >  2:66254/1‑149 (MQ=255)
                                                                                                                      cctcCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAAt                       <  1:668722/148‑1 (MQ=255)
                                                                                                                       ctcCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGa                                                                         <  1:214153/98‑1 (MQ=255)
                                                                                                                       ctcCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGa                                                                         >  2:214153/1‑98 (MQ=255)
                                                                                                                       ctcCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGATCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATg                      <  1:89565/149‑1 (MQ=255)
                                                                                                                         ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGt                    >  2:697266/1‑149 (MQ=255)
                                                                                                                         ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGt                    >  3:33156/1‑149 (MQ=255)
                                                                                                                         ccACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGt                    >  4:576034/1‑149 (MQ=255)
                                                                                                                           aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTaa                  >  1:650430/1‑149 (MQ=255)
                                                                                                                                   aaaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTg          >  1:427716/1‑149 (MQ=255)
                                                                                                                                   aaaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTg          >  2:43841/1‑149 (MQ=255)
                                                                                                                                      tAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCgg       <  3:427098/149‑1 (MQ=255)
                                                                                                                                       ctgcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGa      >  4:360334/4‑149 (MQ=255)
                                                                                                                                        tgcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGAc     >  3:175501/3‑149 (MQ=255)
                                                                                                                                          cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACaa   <  1:43841/149‑1 (MQ=255)
                                                                                                                                          cTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACaa   <  4:396312/149‑1 (MQ=255)
                                                                                                                                           tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAAc  >  1:333430/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
ACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAAC  >  NC_000913/336721‑337008

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: