Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 651,030 C→T 5.5% R276H (CGC→CAC)  citC ← citrate lyase ligase; [citrate [pro‑3S]‑lyase] ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913651,0300CT 5.5% 228.8 / 2.9 72R276H (CGC→CAC) citCcitrate lyase ligase; [citrate [pro‑3S]‑lyase] ligase
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base T (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

CAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTT  >  NC_000913/650885‑651178
                                                                                                                                                 |                                                                                                                                                    
cAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGTGCa                                                                                                                                                   >  1:363069/1‑149 (MQ=255)
cAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCa                                                                                                                                                   >  1:208842/1‑149 (MQ=255)
cAGCGGCGCGATAGCCGTGAGAGCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATATCCTGGGAACGCAGCCGCTCAATTTCACCCAGTTCGAGGGGCGGGGCGGAGATAGTCGGCGTTTCCAGCCAGTAGTGCa                                                                                                                                                   >  1:357219/1‑149 (MQ=255)
        cgATAGCCGTGAGATCGTTTTTCGGCAGCAGTTGACGTACCCGGGAAGCGGTTATCTGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTgg                                                                                                                                           <  5:125656/149‑1 (MQ=255)
                  gagaTCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCCGCCGCTCAATTTCAACCAGGTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGGATATCCTGGTTGTACTggg                                                                                                                                 >  7:391442/1‑149 (MQ=255)
                     aTCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCgg                                                                                                                              >  4:338869/1‑149 (MQ=255)
                     aTCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCgg                                                                                                                              >  3:261182/1‑149 (MQ=255)
                      tCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACACAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCgg                                                                                                                              >  2:262387/1‑148 (MQ=255)
                           tttCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAAcg                                                                                                                         <  4:385713/148‑1 (MQ=255)
                            ttCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAAcgc                                                                                                                        <  1:14650/148‑1 (MQ=255)
                                    cagTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAgg                                                                                                               <  4:97477/149‑1 (MQ=255)
                                    cagTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTAAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAgg                                                                                                               <  2:365109/149‑1 (MQ=255)
                                      gTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGt                                                                                                             <  3:338869/149‑1 (MQ=255)
                                          aCGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCa                                                                                                          <  7:409598/148‑1 (MQ=255)
                                          aCGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAg                                                                                                         >  5:403540/1‑149 (MQ=255)
                                              aCCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTAcc                                                                                                     >  6:53429/1‑149 (MQ=255)
                                               cccGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTAcc                                                                                                     >  8:379500/1‑148 (MQ=255)
                                               cccGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTAcc                                                                                                     >  4:16194/1‑148 (MQ=255)
                                               cccGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTAcc                                                                                                     >  8:373628/1‑148 (MQ=255)
                                                    gAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACaaa                                                                                               <  8:108098/149‑1 (MQ=255)
                                                       gCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCg                                                                                            >  6:402296/1‑149 (MQ=255)
                                                        cGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCgg                                                                                           >  3:390795/1‑149 (MQ=255)
                                                         ggATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGt                                                                                          <  7:373628/149‑1 (MQ=255)
                                                         ggATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGt                                                                                          <  7:379500/149‑1 (MQ=255)
                                                         ggATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGt                                                                                          >  1:217749/1‑149 (MQ=255)
                                                          gATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGa                                                                                                                      >  3:147475/1‑120 (MQ=255)
                                                          gATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGa                                                                                                                      <  4:147475/120‑1 (MQ=255)
                                                           atatCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGt                                                                                          >  4:25789/1‑147 (MQ=255)
                                                             atCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGAc                                                                                                                     <  4:343582/118‑1 (MQ=255)
                                                             atCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGAc                                                                                                                     >  3:343582/1‑118 (MQ=255)
                                                              tCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCACCCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGtt                                                                                     <  3:299412/149‑1 (MQ=255)
                                                              tCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGt                                                                                      <  3:16194/148‑1 (MQ=255)
                                                                   aTCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACAcc                                                                                <  8:401229/149‑1 (MQ=255)
                                                                      tcCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGAc                                                                                                                     <  8:386424/109‑1 (MQ=255)
                                                                      tcCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGAc                                                                                                                     >  7:386424/1‑109 (MQ=255)
                                                                      tcCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTagctcgg                                                                              >  3:52983/1‑142 (MQ=255)
                                                                      tcCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTGCAGTACCGACAAAGCGGTGAGt                                                                                      <  4:52983/140‑1 (MQ=255)
                                                                         tGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAgcgc                                                                          >  7:220022/1‑149 (MQ=255)
                                                                                gCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATACTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGc                                                                    >  3:119698/1‑148 (MQ=255)
                                                                                  atgcGCTCAATTTCAACCAGTTCGTTGGGCGGTGCGGAGATTGTCGGCTTTTCGAGGAAGTTGCGCATATGCTGGTTGTACTGGGCGGTAACGCGACAATAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAg                                                                 <  4:261182/146‑1 (MQ=255)
                                                                                  aGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAg                                                                 <  5:402296/149‑1 (MQ=255)
                                                                                    ccGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAgg                                                                >  4:205591/1‑148 (MQ=255)
                                                                                     cGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTa                                                              >  4:388133/1‑149 (MQ=255)
                                                                                        tCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTg                                                           >  8:278987/1‑149 (MQ=255)
                                                                                          aaTTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGAc                                                         >  5:366457/1‑149 (MQ=255)
                                                                                            tttCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACg                                                        >  1:393950/1‑148 (MQ=255)
                                                                                                     gTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTc                                              <  2:217749/149‑1 (MQ=255)
                                                                                                      ttCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCa                                             <  4:390795/149‑1 (MQ=255)
                                                                                                      ttCGATGGGCGGTGAGGAGATAGTCGGCGTTTACAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCa                                             <  8:414553/149‑1 (MQ=255)
                                                                                                         gATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGa                                           <  3:332506/148‑1 (MQ=255)
                                                                                                                 gTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCg                                  <  8:403348/149‑1 (MQ=255)
                                                                                                                        gaTAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTaa                            >  8:161661/1‑148 (MQ=255)
                                                                                                                         aTAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACa                          <  8:220022/149‑1 (MQ=255)
                                                                                                                            gTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAAt                        <  6:366457/148‑1 (MQ=255)
                                                                                                                            gTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATg                       >  2:188145/1‑149 (MQ=255)
                                                                                                                             tCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATgg                      <  6:403540/149‑1 (MQ=255)
                                                                                                                               ggCGTTTCAAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCTACAAAGCTGTGAGTTACACCCAGCCCGGGAGAGAGGTACTGACGGAAAATCTTCAGATCAATATAGGTGTAACAATgct                     >  3:239524/1‑146 (MQ=255)
                                                                                                                                gCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGtt                    >  8:110125/1‑148 (MQ=255)
                                                                                                                                 cGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTaa                  <  2:393950/149‑1 (MQ=255)
                                                                                                                                  gTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGCGCTAGGTACTGACGGAAAATCTTCAGCTCAATTTCGGTGTAACACTGGTTAAt                 >  3:105335/1‑149 (MQ=255)
                                                                                                                                   tttCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCCGTGTAACAATGGTTAAt                 >  5:338012/1‑148 (MQ=255)
                                                                                                                                       cagccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGc            >  8:55954/1‑149 (MQ=255)
                                                                                                                                       cagccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGc            >  7:78608/1‑149 (MQ=255)
                                                                                                                                       cagccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTATGACGct           >  5:176929/1‑149 (MQ=255)
                                                                                                                                         gccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTc                                              <  4:350351/113‑1 (MQ=255)
                                                                                                                                         gccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTc                                              >  3:350351/1‑113 (MQ=255)
                                                                                                                                         gccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGctc          >  5:298422/1‑149 (MQ=255)
                                                                                                                                          ccagTCGTGCATATCCTGTGTGTACTGGGCGGTAACGCGACAAAAGGGGTCAGTACCGAAAAAGCGGTGATTTCCACCCAGCGCGGGAGCGAGGTACTGAAGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGctct         <  2:357219/149‑1 (MQ=255)
                                                                                                                                          ccagTAGTGCATATCCTTGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCATCGCGTGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGctct         <  2:363069/149‑1 (MQ=255)
                                                                                                                                          ccagTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGctct         <  3:25789/149‑1 (MQ=255)
                                                                                                                                              tAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCtt                                               <  6:5784/107‑1 (MQ=255)
                                                                                                                                              tAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCtt                                               >  5:5784/1‑107 (MQ=255)
                                                                                                                                                gcgcATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCtt   <  4:119698/149‑1 (MQ=255)
                                                                                                                                                 cgcATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCttt  <  3:205591/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAGCGGCGCGATAGCCGTGAGATCGTTTTTCGCCAGCAGTTGACGTACCCGGGAAGCGGATATCGGCATCTCCTGGTAACGCAGCCGCTCAATTTCAACCAGTTCGATGGGCGGTGCGGAGATAGTCGGCGTTTCCAGCCAGTAGCGCATATCCTGGTTGTACTGGGCGGTAACGCGACAAAAGGGTTCAGTACCGACAAAGCGGTGAGTTACACCCAGCGCGGGAGCGAGGTACTGACGGAAAATCTTCAGATCAATTTCGGTGTAACAATGGTTAATGACGCTCTGTTCTTT  >  NC_000913/650885‑651178

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