Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 127,603 C→G 7.6% intergenic (+16/‑309) aceF → / → lpd pyruvate dehydrogenase, E2 subunit/lipoamide dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913127,6030CG 7.6% 348.5 / 15.7 105intergenic (+16/‑309)aceF/lpdpyruvate dehydrogenase, E2 subunit/lipoamide dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (44/53);  new base G (5/3);  total (49/56)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.67e-01

Only 100 of 107 total aligned reads displayed.

AAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACT  >  NC_000913/127454‑127751
                                                                                                                                                     |                                                                                                                                                     
aaGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGccc                                                                                                                                                       >  1:607454/1‑150 (MQ=255)
      tCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggc                                                                                                                                                 <  1:589204/150‑1 (MQ=255)
         tGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggccgg                                                                                                                                              <  1:184769/150‑1 (MQ=255)
          gCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggccggc                                                                                                                                             <  2:757451/150‑1 (MQ=255)
          gCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggccggc                                                                                                                                             >  2:726353/1‑150 (MQ=255)
               gTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCACCGTCTGGTGATGTAAGTAAAAAAGCCGGCCGTTGGGCCGGCtcttt                                                                                                                                       <  2:913226/150‑1 (MQ=255)
                  tGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggcc                                                                                                                                                >  1:180347/1‑139 (MQ=255)
                  tGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAAcggcc                                                                                                                                                <  2:180347/139‑1 (MQ=255)
                   taTGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTg                                                                                                                                   <  2:469719/149‑1 (MQ=255)
                   gATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTg                                                                                                                                   >  1:1645772/1‑150 (MQ=255)
                   gATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTg                                                                                                                                   <  1:1266345/150‑1 (MQ=255)
                   gATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGCAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTg                                                                                                                                   >  1:900909/1‑150 (MQ=255)
                       ctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCGTTGGGCCGGCtcttttactt                                                                                                                                  >  2:1569085/1‑143 (MQ=255)
                       ctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCGTTGGGCCGGCtcttttactt                                                                                                                                  <  1:1569085/147‑5 (MQ=255)
                       ctgcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTaa                                                                                                                               <  1:1094900/150‑1 (MQ=255)
                         gcCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAAtc                                                                                                                             >  2:84513/1‑150 (MQ=255)
                            gATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCa                                                                                                                          <  1:1363646/150‑1 (MQ=255)
                            gATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCa                                                                                                                          <  2:1520486/150‑1 (MQ=255)
                              tttCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTTTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATg                                                                                                                        <  2:1447554/150‑1 (MQ=255)
                              tttCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATg                                                                                                                        <  2:743014/150‑1 (MQ=255)
                                 ctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAAt                                                                                                                     <  1:1506165/150‑1 (MQ=255)
                                 ctctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAAt                                                                                                                     >  1:897778/1‑150 (MQ=255)
                                   ctctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGt                                                                                                                   >  1:615356/1‑150 (MQ=255)
                                     ctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTAt                                                                                                                 >  2:811374/1‑150 (MQ=255)
                                     ctcCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTAt                                                                                                                 >  2:636575/1‑150 (MQ=255)
                                       ccTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTg                                                                                                               >  2:1381786/1‑150 (MQ=255)
                                        cTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGa                                                                                                              <  1:721352/150‑1 (MQ=255)
                                        cTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGa                                                                                                              >  1:1018652/1‑150 (MQ=255)
                                          tCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAgg                                                                                                            <  1:1643040/150‑1 (MQ=255)
                                            gACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGtt                                                                                                          >  2:1531423/1‑150 (MQ=255)
                                              ccagcGCGTGATCGACTGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGttat                                                                                                        <  1:1129493/150‑1 (MQ=255)
                                                 ccGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAg                                                                                                     >  2:1073605/1‑150 (MQ=255)
                                                   gcgTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCg                                                                                                   >  2:335379/1‑150 (MQ=255)
                                                    cgTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGa                                                                                                  >  1:569621/1‑150 (MQ=255)
                                                      tGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAAt                                                                                                <  1:509589/150‑1 (MQ=255)
                                                       ggtCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAAt                                                                                                <  2:128148/147‑1 (MQ=255)
                                                       ggtCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAAt                                                                                                >  1:128148/3‑149 (MQ=255)
                                                       gATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                               >  2:76005/1‑150 (MQ=255)
                                                       gATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                               <  1:811374/150‑1 (MQ=255)
                                                       gATCGACGGTGCTGATGGTGCCCGATTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATa                                                                                               <  2:1547341/150‑1 (MQ=255)
                                                                 gCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCgg                                                                                     >  1:1122119/1‑150 (MQ=255)
                                                                   tGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGtt                                                                                   >  2:1609598/1‑150 (MQ=255)
                                                                   tGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGtt                                                                                   >  2:108706/1‑150 (MQ=255)
                                                                     aTGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGc                                                                                 <  1:1312226/150‑1 (MQ=255)
                                                                       ggTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCg                                                                               >  1:1580424/1‑150 (MQ=255)
                                                                       ggTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCg                                                                               <  1:726353/150‑1 (MQ=255)
                                                                       ggTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTCCCg                                                                               <  2:1635573/150‑1 (MQ=255)
                                                                           cccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttg                                                                           <  1:281439/150‑1 (MQ=255)
                                                                            ccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttgt                                                                          >  2:87601/1‑150 (MQ=255)
                                                                            ccGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTttgt                                                                          >  1:232380/1‑150 (MQ=255)
                                                                                         atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                             <  1:84513/150‑1 (MQ=255)
                                                                                         atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                             >  1:1096682/1‑150 (MQ=255)
                                                                                         atcatTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTg                                                             >  1:263042/1‑150 (MQ=255)
                                                                                           catTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGtt                                                           <  2:268739/150‑1 (MQ=255)
                                                                                              tAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAAc                                                        <  2:1647767/150‑1 (MQ=255)
                                                                                               aacaacACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATc                                                                                       <  2:1078390/118‑1 (MQ=255)
                                                                                               aacaacACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATc                                                                                       >  1:1078390/1‑118 (MQ=255)
                                                                                               aacaacACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACa                                                       >  2:1191431/1‑150 (MQ=255)
                                                                                                    cacGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAAtttt                                                  <  2:607454/150‑1 (MQ=255)
                                                                                                     acGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTg                                                 <  2:569621/150‑1 (MQ=255)
                                                                                                      cGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGt                                                <  1:1277194/150‑1 (MQ=255)
                                                                                                      cGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGt                                                <  1:1073605/150‑1 (MQ=255)
                                                                                                       gCTGTCTGACATTCGCCGTCTGGTGATGCAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTa                                               <  2:900909/150‑1 (MQ=255)
                                                                                                          gTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTaaaa                                            <  1:1451998/150‑1 (MQ=255)
                                                                                                                cATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGa                                      >  2:769109/1‑150 (MQ=255)
                                                                                                                cATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGa                                      <  1:76005/150‑1 (MQ=255)
                                                                                                                    cGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGa                                  <  1:1531423/150‑1 (MQ=255)
                                                                                                                     gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGAt                                 >  2:1226194/1‑150 (MQ=255)
                                                                                                                     gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGAt                                 <  1:1171952/150‑1 (MQ=255)
                                                                                                                     gCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGAt                                 <  2:233514/150‑1 (MQ=255)
                                                                                                                      ccGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                                <  2:1337163/150‑1 (MQ=255)
                                                                                                                      ccGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATa                                <  2:1096682/150‑1 (MQ=255)
                                                                                                                       cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAg                               <  2:897778/150‑1 (MQ=255)
                                                                                                                       cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAg                               <  1:1609598/150‑1 (MQ=255)
                                                                                                                       cGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAg                               <  2:74201/150‑1 (MQ=255)
                                                                                                                          cTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAAc                            <  1:1181455/150‑1 (MQ=255)
                                                                                                                           tGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACg                           >  2:1589726/1‑150 (MQ=255)
                                                                                                                            ggTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGa                          >  1:979150/1‑150 (MQ=255)
                                                                                                                            ggTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGa                          >  1:1486758/1‑150 (MQ=255)
                                                                                                                            ggTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGa                          >  1:232779/1‑150 (MQ=255)
                                                                                                                              tGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGAcc                        >  1:432530/1‑150 (MQ=255)
                                                                                                                                 tGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                     >  1:386110/1‑150 (MQ=255)
                                                                                                                                 tGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCgg                     >  1:485137/1‑150 (MQ=255)
                                                                                                                                   accagaaaaaaGCCGGCCGTTGGGCCGGCTCTTTTACTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggt                    >  1:1005565/6‑150 (MQ=38)
                                                                                                                                    ccagaaaaaaGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                  >  2:563861/5‑150 (MQ=255)
                                                                                                                                    ccagaaaaaaGCCGGCCGTTGGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                  >  1:691937/5‑150 (MQ=255)
                                                                                                                                    ccagaaaaaaGCCGGCCGTTGGGCCGGCTCTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                  >  2:1312044/5‑150 (MQ=255)
                                                                                                                                    aagtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                  <  1:769109/150‑1 (MQ=255)
                                                                                                                                    aactaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtgg                  <  2:420372/147‑1 (MQ=255)
                                                                                                                                     agtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggt                 >  2:348739/1‑150 (MQ=255)
                                                                                                                                      gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACACggtggtg                >  2:1078865/1‑150 (MQ=255)
                                                                                                                                      gtaaAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATCTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCggtggtg                >  2:1084532/1‑150 (MQ=255)
                                                                                                                                         aaaGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGtt             >  1:1035446/1‑150 (MQ=255)
                                                                                                                                            gagCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAgg          <  2:1265232/150‑1 (MQ=255)
                                                                                                                                             agCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAggg         <  1:87601/150‑1 (MQ=255)
                                                                                                                                             agCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAggg         <  2:432530/150‑1 (MQ=255)
                                                                                                                                              gCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGt        <  2:1035446/150‑1 (MQ=255)
                                                                                                                                                 ggCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACg                                    >  1:638517/1‑119 (MQ=255)
                                                                                                                                                 ggCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACg                                    <  2:638517/119‑1 (MQ=255)
                                                                                                                                                    ccAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATATAAAGTACGGGTGGTGGTTATTGTAtttat  <  2:1580424/150‑4 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGAGTTCGTGCCGCGTCTGATGCTGCCGATTTCTCTCTCCTTCGACCACCGCGTGATCGACGGTGCTGATGGTGCCCGTTTCATTACCATCATTAACAACACGCTGTCTGACATTCGCCGTCTGGTGATGTAAGTAAAAGAGCCGGCCCAACGGCCGGCTTTTTTCTGGTAATCTCATGAATGTATTGAGGTTATTAGCGAATAGACAAATCGGTTGCCGTTTGTTGTTTAAAAATTGTTAACAATTTTGTAAAATACCGACGGATAGAACGACCCGGTGGTGGTTAGGGTATTACT  >  NC_000913/127454‑127751

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