Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 127,789:1 +T 6.2% intergenic (+202/‑123) aceF → / → lpd pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2/dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913127,7891.T 6.2% 411.3 / 22.3 96intergenic (+202/‑123)aceF/lpdpyruvate dehydrogenase, dihydrolipoyltransacetylase component E2/dihydrolipoyl dehydrogenase; E3 component of pyruvate and 2‑oxoglutarate dehydrogenases complexes; glycine cleavage system L protein; dihydrolipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (44/46);  new base T (3/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-02

GAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCAAAACTC  >  NC_000913/127642‑127933
                                                                                                                                                     |                                                                                                                                                
gAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGc                                                                                                                                                   <  2:505240/148‑1 (MQ=255)
gAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGc                                                                                                                                                   >  1:505240/1‑148 (MQ=255)
  ggTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAg                                                                                                                                               >  4:935087/1‑149 (MQ=255)
    ttattaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTTTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCG                                                                                                                                             >  2:1595614‑M1/1‑149 (MQ=255)
       ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          <  3:935087/149‑1 (MQ=255)
       ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          >  4:1320871/1‑149 (MQ=255)
       ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          >  4:719677/1‑149 (MQ=255)
           cGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGTGTGCAGCAAGGTAGCAAGCGCCAGAAt                                                                                                                                      <  1:781602/149‑1 (MQ=255)
             aaTAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATcc                                                                                                                                    <  2:1320425/149‑1 (MQ=255)
             aaTAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCCCCAGAATcc                                                                                                                                    <  3:1320871/149‑1 (MQ=255)
                  aCAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCTTCACGCCctcaccccagt                                                                                                                                >  2:1078765‑M1/1‑138 (MQ=255)
                  aCAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAgg                                                                                                                               >  1:1283575/1‑149 (MQ=255)
                    aaaTCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAg                                                                                                                             >  4:1368188/1‑149 (MQ=255)
                     aaTCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGc                                                                                                                            <  1:1476247/149‑1 (MQ=255)
                      aTCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           <  1:1233387/149‑1 (MQ=255)
                             gCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    <  2:1283575/149‑1 (MQ=255)
                                   tgttgtTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATaagtaagt                                                                                                              <  1:1595614/149‑1 (MQ=255)
                                   tgttgtTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATaagtaagt                                                                                                              <  3:1368188/149‑1 (MQ=255)
                                    gttgttTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATcc                                                                                                                                    >  4:320374/1‑126 (MQ=255)
                                    gttgttTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATcc                                                                                                                                    <  3:320374/126‑1 (MQ=255)
                                     ttgttTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCTTCACGCCctcaccccagtcacttacttatgtaagctc                                                                                                             <  1:1078765‑M1/149‑31 (MQ=255)
                                        tttAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCctcaccccagtcacttacttatgtaagctcctgg                                                                                                         >  1:709763‑M1/1‑115 (MQ=255)
                                        tttAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTg                                                                                                         >  2:1257264/1‑149 (MQ=255)
                                              aaTTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCctcaccccagtcacttacttatgtaagctcctggggattc                                                                                                   <  2:709763‑M1/149‑41 (MQ=255)
                                                   ttAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACa                                                                                                                       <  3:51494/124‑1 (MQ=255)
                                                   ttAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACa                                                                                                                       >  4:51494/1‑124 (MQ=255)
                                                   ttAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAg                                                                                                  >  2:1284240/1‑145 (MQ=255)
                                                   ttAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAg                                                                                                  <  1:1284240/145‑1 (MQ=255)
                                                       cAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaagt                                                                                                                  >  4:1103826/1‑125 (MQ=255)
                                                       cAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaagt                                                                                                                  <  3:1103826/125‑1 (MQ=255)
                                                        aaTTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATcccc                                                                                                                                  <  1:713219/108‑1 (MQ=255)
                                                        aaTTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATcccc                                                                                                                                  >  2:713219/1‑108 (MQ=255)
                                                          ttttGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCt                                                                                       <  3:719677/149‑1 (MQ=255)
                                                             tGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCa                                                                                                                                 <  3:565905/104‑1 (MQ=255)
                                                             tGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCa                                                                                                                                 >  4:565905/1‑104 (MQ=255)
                                                               tAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGa                                                                                                            >  4:198709/1‑123 (MQ=255)
                                                               tAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGa                                                                                                            <  3:198709/123‑1 (MQ=255)
                                                                  aaTACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTg                                                                                            <  3:353194/136‑1 (MQ=255)
                                                                  aaTACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTg                                                                                            >  4:353194/1‑136 (MQ=255)
                                                                     aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 <  4:1340577‑M1/146‑59 (MQ=255)
                                                                     aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 >  3:1340577‑M1/1‑88 (MQ=255)
                                                                     aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 <  4:185313‑M1/146‑59 (MQ=255)
                                                                     aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 >  3:185313‑M1/1‑88 (MQ=255)
                                                                        gACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGc                                                                         <  1:1257264/149‑1 (MQ=255)
                                                                         aCGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATccc                                                                                                                                   >  1:765322/1‑90 (MQ=255)
                                                                         aCGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATccc                                                                                                                                   <  2:765322/90‑1 (MQ=255)
                                                                             aTAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGa                                                                                                                              >  3:143562/1‑91 (MQ=255)
                                                                             aTAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAAGa                                                                                                                              <  4:143562/91‑1 (MQ=255)
                                                                               aGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGc                                                                                                                                                   <  4:1017779/69‑1 (MQ=255)
                                                                               aGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGc                                                                                                                                                   >  3:1017779/1‑69 (MQ=255)
                                                                                 aaCGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           >  3:234875/1‑90 (MQ=255)
                                                                                 aaCGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           <  4:234875/90‑1 (MQ=255)
                                                                                     aCCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          <  3:1523214/119‑1 (MQ=255)
                                                                                     aCCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          >  4:1523214/1‑119 (MQ=255)
                                                                                                      gTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    <  3:219612/76‑1 (MQ=255)
                                                                                                      gTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    >  4:219612/1‑76 (MQ=255)
                                                                                                       tATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTgg                                                                                                        <  3:1781741/87‑1 (MQ=255)
                                                                                                       tATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTgg                                                                                                        >  4:1781741/1‑87 (MQ=255)
                                                                                                        aTTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGt                                                         >  3:382986/1‑133 (MQ=255)
                                                                                                        aTTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGt                                                         <  4:382986/133‑1 (MQ=255)
                                                                                                          tACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATaaa                                       >  2:1274204/1‑149 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    <  3:363322/69‑1 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    >  4:363322/1‑69 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTg                                                                                                    <  2:1401374/85‑1 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTg                                                                                                    >  1:1401374/1‑85 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCgg                                                                          <  4:575289/111‑1 (MQ=255)
                                                                                                             ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCgg                                                                          >  3:575289/1‑111 (MQ=255)
                                                                                                                acaTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           >  3:1060276/1‑59 (MQ=255)
                                                                                                                acaTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           <  4:1060276/59‑1 (MQ=255)
                                                                                                                  aTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATaga                               >  1:675159/1‑149 (MQ=255)
                                                                                                                  aTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATaga                               >  2:1620956/1‑149 (MQ=255)
                                                                                                                     cccTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGt                            >  4:984226/1‑149 (MQ=255)
                                                                                                                      ccTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCtt                                                                                                                          <  4:664069/54‑1 (MQ=255)
                                                                                                                      ccTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCtt                                                                                                                          >  3:664069/1‑54 (MQ=255)
                                                                                                                       cTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAg                                                                                                                                <  4:916563/47‑1 (MQ=255)
                                                                                                                       cTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAg                                                                                                                                >  3:916563/1‑47 (MQ=255)
                                                                                                                        tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgctt                                                                                           <  2:1283177‑M1/83‑49 (MQ=255)
                                                                                                                        tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgctt                                                                                           >  1:1283177‑M1/1‑35 (MQ=255)
                                                                                                                        tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTc                           <  1:128318/147‑1 (MQ=255)
                                                                                                                        tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTc                           >  2:128318/1‑147 (MQ=255)
                                                                                                                        tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCat                         >  4:1206666/1‑149 (MQ=255)
                                                                                                                          tGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCatga                       >  4:345781/1‑149 (MQ=255)
                                                                                                                             aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          >  4:107707/1‑79 (MQ=255)
                                                                                                                             aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          <  3:107707/79‑1 (MQ=255)
                                                                                                                             aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCAtgatga                    >  3:885625/1‑149 (MQ=255)
                                                                                                                             aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCAtgatga                    <  3:984226/149‑1 (MQ=255)
                                                                                                                                tCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTa                 <  4:885625/149‑1 (MQ=255)
                                                                                                                                         agcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAata                                     <  1:1003126/120‑1 (MQ=255)
                                                                                                                                         agcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAata                                     >  2:1003126/1‑120 (MQ=255)
                                                                                                                                         agcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCatga                       >  3:802023/1‑134 (MQ=255)
                                                                                                                                         agcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCatga                       <  4:802023/134‑1 (MQ=255)
                                                                                                                                          gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          <  4:185313‑M2/146‑89 (MQ=255)
                                                                                                                                          gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          >  3:185313‑M2/1‑58 (MQ=255)
                                                                                                                                          gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          >  3:1340577‑M2/1‑58 (MQ=255)
                                                                                                                                          gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          <  4:1340577‑M2/146‑89 (MQ=255)
                                                                                                                                             aaGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCaaaa     <  2:675159/149‑1 (MQ=255)
                                                                                                                                             aaGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCaaaa     <  1:1620956/149‑1 (MQ=255)
                                                                                                                                               ggTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGa                                                                     >  3:1419205/1‑83 (MQ=255)
                                                                                                                                               ggTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGa                                                                     <  4:1419205/83‑1 (MQ=255)
                                                                                                                                                gTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCAAAACTc  <  3:1206666/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCAAAACTC  >  NC_000913/127642‑127933

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