Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 127,791 A→C 6.2% intergenic (+204/‑121) 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,7910AC 6.2% 345.7 / 13.8 96intergenic (+204/‑121)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 A (44/46);  new base C (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 = 7.01e-01

GGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCAAAACTC  >  NC_000913/127644‑127933
                                                                                                                                                     |                                                                                                                                              
ggTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAg                                                                                                                                               >  4:935087/1‑149 (MQ=255)
  ttattaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTTTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCG                                                                                                                                             >  2:1595614‑M1/1‑149 (MQ=255)
     ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          >  4:719677/1‑149 (MQ=255)
     ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          >  4:1320871/1‑149 (MQ=255)
     ttaGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCa                                                                                                                                          <  3:935087/149‑1 (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                                                                                 <  4:185313‑M1/146‑59 (MQ=255)
                                                                   aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 >  3:185313‑M1/1‑88 (MQ=255)
                                                                   aCCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCGGCGTTAGCTTCACGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgcttgctaccttgc                                                                                 >  3:1340577‑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)
                                                                               aaCGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           <  4:234875/90‑1 (MQ=255)
                                                                               aaCGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           >  3:234875/1‑90 (MQ=255)
                                                                                   aCCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          <  3:1523214/119‑1 (MQ=255)
                                                                                   aCCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          >  4:1523214/1‑119 (MQ=255)
                                                                                                    gTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    >  4:219612/1‑76 (MQ=255)
                                                                                                    gTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACAtaa                                                                                                                    <  3:219612/76‑1 (MQ=255)
                                                                                                     tATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTgg                                                                                                        >  4:1781741/1‑87 (MQ=255)
                                                                                                     tATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTgg                                                                                                        <  3:1781741/87‑1 (MQ=255)
                                                                                                      aTTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGt                                                         <  4:382986/133‑1 (MQ=255)
                                                                                                      aTTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGt                                                         >  3:382986/1‑133 (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                                                                          >  3:575289/1‑111 (MQ=255)
                                                                                                           ttCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCgg                                                                          <  4:575289/111‑1 (MQ=255)
                                                                                                              acaTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           <  4:1060276/59‑1 (MQ=255)
                                                                                                              acaTACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCt                                                                                                                           >  3:1060276/1‑59 (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                                                                                           >  1:1283177‑M1/1‑35 (MQ=255)
                                                                                                                      tATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCctcaccccagtcacttacttatgtaagctcctggggattctggcgctt                                                                                           <  2:1283177‑M1/83‑49 (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                                                                                          <  3:107707/79‑1 (MQ=255)
                                                                                                                           aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGaa                                                                                          >  4:107707/1‑79 (MQ=255)
                                                                                                                           aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCAtgatga                    <  3:984226/149‑1 (MQ=255)
                                                                                                                           aTTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCAtgatga                    >  3:885625/1‑149 (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          >  3:185313‑M2/1‑58 (MQ=255)
                                                                                                                                        gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          <  4:185313‑M2/146‑89 (MQ=255)
                                                                                                                                        gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          >  3:1340577‑M2/1‑58 (MQ=255)
                                                                                                                                        gcaAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgtgaagctaacgccgctgcacccagaaatccatagggtatgtgaagtaataccctaaccaccaccgggtcgttctatccgtcggt          <  4:1340577‑M2/146‑89 (MQ=255)
                                                                                                                                           aaGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCaaaa     <  1:1620956/149‑1 (MQ=255)
                                                                                                                                           aaGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCaaaa     <  2:675159/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)
                                                                                                                                                    aaGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgcttgctaccttgctgcacccagaaatccata                                                               <  2:1283177‑M2/83‑36 (MQ=255)
                                                                                                                                                    aaGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGggcgcttgctaccttgctgcacccagaaatccata                                                               >  1:1283177‑M2/1‑48 (MQ=255)
                                                                                                                                                     aGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGcc                                                                        >  3:1164935/1‑73 (MQ=255)
                                                                                                                                                     aGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGCGTTAGGGCGTGAAGCTAACGCCGCTGCGGcc                                                                        <  4:1164935/73‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACTTCACATACCCTATGGATTTCTGGGTGCAGCAAGGTAGCAAGCGCCAGAATCCCCAGGAGCTTACATAAGTAAGTGACTGGGGTGAGGGCGTGAAGCTAACGCCGCTGCGGCCTGAAAGACGACGGGTATGACCGCCGGAGATAAATATATAGAGGTCATGATGAGTACTGAAATCAAAACTC  >  NC_000913/127644‑127933

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