Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 336,866 G→C 7.0% 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 7.0% 244.2 / 8.5 71G30A (GGC→GCC) yahFputative NAD(P)‑binding succinyl‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base C (2/3);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCA  >  NC_000913/336721‑337013
                                                                                                                                                 |                                                                                                                                                   
aCATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGt                                                                                                                                                  <  1:412036/149‑1 (MQ=255)
     cTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGc                                                                                                                                             >  3:520265/1‑149 (MQ=255)
     cTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgt                                                                                                                                             <  1:224360/149‑2 (MQ=255)
         atggcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgtttat                                                                                                                                         <  3:625848/149‑6 (MQ=255)
            gcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtt                                                                                                                                      >  3:686439/1‑149 (MQ=255)
            gcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgt                                                                                                                                             <  3:513151/142‑2 (MQ=255)
            gcatggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGCCGTCGAgt                                                                                                                                             >  4:513151/1‑141 (MQ=255)
             catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAttt                                                                                                                                     >  1:523287/1‑149 (MQ=255)
             catggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGGCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCAtct                                                                                                                                     >  4:592869/1‑147 (MQ=255)
                ggcCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTtgtg                                                                                                                                  >  4:102165/1‑149 (MQ=255)
                     aTGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGtggcgat                                                                                                                             <  1:455691/149‑1 (MQ=255)
                                    gCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGaa                                                                                                              <  2:230216/149‑1 (MQ=255)
                                     ccAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAAt                                                                                                             >  2:567240/1‑149 (MQ=255)
                                     ccAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAAt                                                                                                             <  2:523287/149‑1 (MQ=255)
                                       aaaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCCGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaa                                                                                                           >  4:618285/1‑149 (MQ=255)
                                        aaCCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATaaa                                                                                                          >  4:158947/1‑149 (MQ=255)
                                            tACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCg                                                                                                      >  2:15374/1‑149 (MQ=255)
                                                gAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCt                                                                                                  >  4:208212/1‑149 (MQ=255)
                                                  gggAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTga                                                                                                >  3:384603/1‑149 (MQ=255)
                                                       ttATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAt                                                                                           <  2:341585/149‑1 (MQ=255)
                                                          tGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                        <  4:686439/149‑1 (MQ=255)
                                                          tGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                        <  4:393480/149‑1 (MQ=255)
                                                          tGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCCGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTa                                                                                        <  3:618285/149‑1 (MQ=255)
                                                            tCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAgg                                                                                      >  4:539274/1‑149 (MQ=255)
                                                               gTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                         <  2:490877/127‑1 (MQ=255)
                                                               gTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                         >  1:490877/1‑127 (MQ=255)
                                                                     aTAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGCCTGCTGAc                                                                             >  3:244596/1‑149 (MQ=255)
                                                                         tCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCg                                                                         <  4:520265/149‑1 (MQ=255)
                                                                           aTTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGa                                                                       >  2:46720/1‑149 (MQ=255)
                                                                            ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAg                                                                      >  2:175124/1‑149 (MQ=255)
                                                                            ttAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAg                                                                      >  1:292327/1‑149 (MQ=255)
                                                                               aaCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTg                                                                   >  2:405754/1‑149 (MQ=255)
                                                                               aaCCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTg                                                                   >  1:383540/1‑149 (MQ=255)
                                                                                 ccGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGa                                                                 <  3:102165/149‑1 (MQ=255)
                                                                                         cTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGa                                                         >  1:338840/1‑149 (MQ=255)
                                                                                           aTTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGaaa                                                       >  1:519172/1‑149 (MQ=255)
                                                                                             tttGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGaaaaa                                                     <  1:95792/149‑1 (MQ=255)
                                                                                                 aTTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGc                                                 <  2:519172/149‑1 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                         >  2:546060/1‑92 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAg                                                                                                         <  1:546060/92‑1 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAgg                                                            >  2:78762/1‑137 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAgg                                                            <  1:78762/137‑1 (MQ=255)
                                                                                                  ttCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCg                                                >  1:75923/1‑149 (MQ=255)
                                                                                                    cTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCTGGAGCTGGAGCAGGCGAAAAACGGCGAc                                              <  3:592869/149‑1 (MQ=255)
                                                                                                       tctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtg                                           >  2:432180/1‑149 (MQ=255)
                                                                                                        ctcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtga                                          <  2:75923/149‑1 (MQ=255)
                                                                                                         tcGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtgat                                         <  4:244596/149‑1 (MQ=255)
                                                                                                           gCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                                         <  3:295074/99‑1 (MQ=255)
                                                                                                           gCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAAttt                                                                                         >  4:295074/1‑99 (MQ=255)
                                                                                                             tGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGActgct                                                                                <  2:313661/106‑1 (MQ=255)
                                                                                                             tGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGActgct                                                                                >  1:313661/1‑106 (MQ=255)
                                                                                                             tGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGAtt                                     <  2:630702/149‑1 (MQ=255)
                                                                                                              gATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTgg                                                                  >  2:234434/1‑119 (MQ=255)
                                                                                                              gATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTgg                                                                  <  1:234434/119‑1 (MQ=255)
                                                                                                                   cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAAc                                                    >  3:242317/1‑128 (MQ=255)
                                                                                                                   cTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAAc                                                    <  4:242317/128‑1 (MQ=255)
                                                                                                                       ctcCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATg                           <  2:383540/149‑1 (MQ=255)
                                                                                                                          cACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTa                        >  3:221358/1‑149 (MQ=255)
                                                                                                                           aCGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTaa                       >  4:216883/1‑149 (MQ=255)
                                                                                                                                tGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCgg                  <  1:46720/149‑1 (MQ=255)
                                                                                                                                 gCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGActg                                                                                  >  3:160655/1‑84 (MQ=255)
                                                                                                                                 gCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGActg                                                                                  <  4:160655/84‑1 (MQ=255)
                                                                                                                                   aaaTAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTg               >  2:568928/1‑149 (MQ=255)
                                                                                                                                       ctgcTCGACGCCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGa           >  2:266412/4‑149 (MQ=255)
                                                                                                                                           tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtgatg                                        >  2:644213/1‑116 (MQ=255)
                                                                                                                                           tCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCtgatg                                        <  1:644213/116‑1 (MQ=255)
                                                                                                                                             gACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAgg                                                            >  3:86006/1‑94 (MQ=255)
                                                                                                                                             gACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAgg                                                            <  4:86006/94‑1 (MQ=255)
                                                                                                                                              aCGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAg    >  2:187643/1‑149 (MQ=255)
                                                                                                                                              aCGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAg    <  1:584830/149‑1 (MQ=255)
                                                                                                                                              aCGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAg    <  1:568928/149‑1 (MQ=255)
                                                                                                                                                ggCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCa  <  3:216883/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ACATGCTGTATGGCATGGCCGATGGTTGTGCGCTGAGCCAAACCTACGGAGGGAATTATGTCAGTTAAAATAGTCATTAAACCGAATACCTATTTTGATTCTGTCTCGCTGATGTCTATCTCCACGCGTGCAAATAAACTCGACGGCGTCGAGCAGGCATTTGTGGCGATGGCGACCGAAATGAATAAAGGCGTGCTGAAGAATTTAGGACTGCTGACGCCGGAGCTGGAGCAGGCGAAAAACGGCGACCTGATGATTGTCATCAATGGTAAATCGGGTGCGGACAACGAGCA  >  NC_000913/336721‑337013

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