Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 127,603 C→G 7.5% intergenic (+16/‑309) 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,6030CG 7.5% 265.3 / 12.3 80intergenic (+16/‑309)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 C (37/37);  new base G (2/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

TTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTAC  >  NC_000913/127459‑127750
                                                                                                                                                |                                                                                                                                                    
ttCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGAGCCCGTTTCATTACCATCATTAACAACACGCGGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcg                                                                                                                                                  <  4:83707/149‑1 (MQ=255)
    tGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggccg                                                                                                                                              >  4:196517/1‑149 (MQ=255)
            cTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTCACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtttttt                                                                                                                                      <  1:55007/149‑1 (MQ=255)
               aTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCtttttt                                                                                                                                      <  4:272200/146‑1 (MQ=255)
               aTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTg                                                                                                                                  <  5:31287/149‑1 (MQ=255)
                 gctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTgg                                                                                                                                 >  4:11457/1‑148 (MQ=255)
                       gATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAAtct                                                                                                                           >  1:25613/1‑148 (MQ=255)
                           tctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATg                                                                                                                       <  7:244926/148‑1 (MQ=255)
                             tctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAAt                                                                                                                    >  8:292265/1‑149 (MQ=255)
                                 tcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCCGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGt                                                                                                                  <  4:202658/147‑1 (MQ=255)
                                   cTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTAtt                                                                                                               >  8:246348/1‑148 (MQ=255)
                                       gACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAgg                                                                                                           <  3:258967/148‑1 (MQ=255)
                                       gACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGt                                                                                                          >  6:328726/1‑149 (MQ=255)
                                             cgcgTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACACCACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAg                                                                                                    <  4:158319/149‑1 (MQ=255)
                                                  gATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGaa                                                                                                >  7:122394/1‑148 (MQ=255)
                                                   aTCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTACAAGAGCCGGCCCAACGGCCGGCGTTTTTCTGTTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                              <  6:297129/149‑1 (MQ=255)
                                                   aTCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                              <  7:50999/149‑1 (MQ=255)
                                                     cGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGa                                                                                            >  8:298487/1‑149 (MQ=255)
                                                       aCGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGAc                                                                                           >  7:158732/1‑148 (MQ=255)
                                                           aatttttaaacaacaaacggcaaccgatttgtctattcgctaattacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGttatta                                                                                                     <  7:162021‑M1/57‑1 (MQ=255)
                                                           aatttttaaacaacaaacggcaaccgatttgtctattcgctaattacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGttatta                                                                                                     >  8:162021‑M1/78‑134 (MQ=255)
                                                                aTGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGtt                                                                                  <  7:246348/148‑1 (MQ=255)
                                                                aTGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTg                                                                                 >  4:278348/1‑149 (MQ=255)
                                                                    tGCCCGTTTCATTAACATCATTAACAACACGCTGTCTGAAATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCAGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGACGt                                                                             >  7:95958/1‑149 (MQ=255)
                                                                          tttCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttgttg                                                                       >  2:89983/1‑149 (MQ=255)
                                                                              aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCgg                                                                                    <  7:96376/132‑1 (MQ=255)
                                                                              aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCgg                                                                                    >  8:96376/1‑132 (MQ=255)
                                                                              aTTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGttt                                                                    <  3:196517/148‑1 (MQ=255)
                                                                               ttACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGtt                                                                                                         <  6:7697/110‑1 (MQ=255)
                                                                                 aCCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTaaa                                                                 <  3:212136/148‑1 (MQ=255)
                                                                                   catcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAt                                                              >  3:296931/1‑149 (MQ=255)
                                                                                    atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAtt                                                             >  6:330099/1‑149 (MQ=255)
                                                                                    atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAAtt                                                             <  5:328726/149‑1 (MQ=255)
                                                                                       atTAACAACACGCTGTCTGACAATCGCCGTCTGGTCATGTAAGTAAAAGAGCCGGCCCCACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGtt                                                          <  4:158725/149‑1 (MQ=255)
                                                                                        ttAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGtt                                                          >  1:11233/1‑148 (MQ=255)
                                                                                            caacaCGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACaa                                                     >  6:304094/1‑149 (MQ=255)
                                                                                             attcgctaataacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTCTTTTACTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACaa                                                     <  6:101030‑M1/106‑1 (MQ=255)
                                                                                             aacaCGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGGCGTTTGTTGTTTAAAAATTGTTTACAAt                                                    >  3:181552/1‑149 (MQ=255)
                                                                                                acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAtttt                                                 <  6:89045/149‑1 (MQ=255)
                                                                                                acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAtttt                                                 <  2:290188/149‑1 (MQ=255)
                                                                                                  gCTGTCTGACATTCGTCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                >  7:15608/1‑148 (MQ=255)
                                                                                                  gCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGt                                               <  3:278348/149‑1 (MQ=255)
                                                                                                  ctaataacctcaatacattcatgagattaccagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                <  8:2984‑M1/110‑1 (MQ=255)
                                                                                                        tGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTTACAATTTTGTAAAATa                                         >  5:222/1‑149 (MQ=255)
                                                                                                        tGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTCGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATa                                         >  1:51338/1‑149 (MQ=255)
                                                                                                         gACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAc                                        >  1:270716/1‑149 (MQ=255)
                                                                                                          aCATTCGCCGTCTGGTTATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATTTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAc                                        <  3:107089/148‑1 (MQ=255)
                                                                                                            aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTa                                              >  4:79523/1‑140 (MQ=255)
                                                                                                            aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTa                                              <  3:79523/140‑1 (MQ=255)
                                                                                                            aTTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGa                                     >  7:152620/1‑149 (MQ=255)
                                                                                                             ttCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                    <  7:36978/149‑1 (MQ=255)
                                                                                                             ttCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                    >  3:243031/1‑149 (MQ=255)
                                                                                                               cGCCGTCTGGCGATGTAAGTAAACGAGCCGGCCCACCGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGCTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGAc                                    <  4:181552/147‑1 (MQ=255)
                                                                                                                gCCGTCTGTTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                 <  7:292265/149‑1 (MQ=255)
                                                                                                                gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATAACGACGGa                                 >  4:134158/1‑149 (MQ=255)
                                                                                                                gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGACGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACgg                                  <  8:95958/148‑1 (MQ=255)
                                                                                                                  cGTCTTGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                               <  8:229412/149‑1 (MQ=255)
                                                                                                                  cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                               <  5:297543/149‑1 (MQ=255)
                                                                                                                    tCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGa                             <  4:243031/149‑1 (MQ=255)
                                                                                                                        gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                            <  8:158732/146‑1 (MQ=255)
                                                                                                                        gTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                            <  1:332709/146‑1 (MQ=255)
                                                                                                                             gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                 >  1:16583/1‑136 (MQ=255)
                                                                                                                             gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                 <  2:16583/136‑1 (MQ=255)
                                                                                                                             gtaagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                    >  3:109517/1‑149 (MQ=255)
                                                                                                                                cagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                            >  6:65506‑M1/9‑138 (MQ=255)
                                                                                                                                cagaaaaaAGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGaa                            <  5:65506‑M1/130‑1 (MQ=255)
                                                                                                                                agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTGTGTTGTTTAAAAATTGTTATCAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                 >  5:16977/1‑149 (MQ=255)
                                                                                                                                 gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                >  1:270779/1‑149 (MQ=255)
                                                                                                                                 gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                >  3:159138/1‑149 (MQ=255)
                                                                                                                                  taaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg               >  7:59028/1‑149 (MQ=255)
                                                                                                                                  taaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg               <  8:152620/149‑1 (MQ=255)
                                                                                                                                  taaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg               >  7:102620/1‑149 (MQ=255)
                                                                                                                                   aaaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg               <  2:270716/148‑1 (MQ=255)
                                                                                                                                    aaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                            >  2:134040/1‑102 (MQ=255)
                                                                                                                                    aaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                            <  1:134040/102‑1 (MQ=255)
                                                                                                                                            GGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTAt     >  2:25613‑M1/1‑149 (MQ=255)
                                                                                                                                             GCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTAt     <  7:14248‑M1/148‑1 (MQ=255)
                                                                                                                                              CCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTa   >  1:158180‑M1/1‑149 (MQ=255)
                                                                                                                                              CCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTa   >  4:296931‑M1/1‑149 (MQ=255)
                                                                                                                                               CCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTAc  <  3:32282‑M1/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTAC  >  NC_000913/127459‑127750

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