Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,791,526 G→A A342A (GCC→GCT kbl ← 2‑amino‑3‑ketobutyrate CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,791,5260GA98.3% 186.7 / NA 59A342A (GCC→GCTkbl2‑amino‑3‑ketobutyrate CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (34/24);  minor base T (0/1);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

CGGTACAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTAAATGCTTCTACTGCACGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTAACGTAAATGCCCTCTTTTTGCAGCTCACGGGCAAATTTCTGCGCCACTACCGCATCACCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCA  >  NC_000913/3791283‑3791769
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
cggctcAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGGCATCTGGGTACGCATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGAAGTAAGTGCCCTCTTTTTGCAGCTCCCGAGCAAAtt                                                                                                                                                                                                                                              >  2:425578/6‑251 (MQ=255)
     cAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTgc                                                                                                                                                                                                                                         >  1:204681/1‑251 (MQ=255)
              cGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTAcc                                                                                                                                                                                                                                >  2:387391/1‑251 (MQ=255)
              cGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTAcc                                                                                                                                                                                                                                >  2:137151/1‑251 (MQ=255)
              cGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTAc                                                                                                                                                                                                                                 <  2:21008/250‑1 (MQ=255)
               ggTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCg                                                                                                                                                                                                                               >  1:305744/1‑251 (MQ=255)
                gTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGc                                                                                                                                                                                                                              >  2:121210/1‑251 (MQ=255)
                   aTCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCAGTTCACATCCTCAGGCGATAACTCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATc                                                                                                                                                                                                                           <  1:391556/251‑1 (MQ=255)
                   aTCCAGATGCCCTCTGCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATc                                                                                                                                                                                                                           <  1:217837/251‑1 (MQ=255)
                        gATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGCCATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGCAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCa                                                                                                                                                                                                                       >  2:111664/1‑250 (MQ=255)
                         aTGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCa                                                                                                                                                                                                                       >  1:351699/1‑249 (MQ=255)
                           gCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCa                                                                                                                                                                                                                   <  1:105841/251‑1 (MQ=255)
                            cccTCTTCCGCTGTCAGTATGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCTCAGTTGTTTACCGATACGCGGGAATGCTTCTACTGCGCGCGTAATGTGCTCACGTGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACCTAAATGCCCTCGTTTTGCAGCGCACGTGCAAATTTCTGCGCCCCATCCGCATCGCTAAGCAt                                                                                                                                                                                                                  <  2:190518/251‑1 (MQ=255)
                             ccTCTTCCGCTTTCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATg                                                                                                                                                                                                                 >  1:355596/1‑251 (MQ=255)
                                ctTCCGCTTTCAGTTTGGATAACGCTTTCAGTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTTCTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGTCCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGa                                                                                                                                                                                                                <  1:295608/249‑1 (MQ=255)
                                        ttCAGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGaata                                                                                                                                                                                                        <  1:71875/249‑1 (MQ=255)
                                          cAGTTTGGATAACGCGTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGTCATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATg                                                                                                                                                                                                     <  2:342062/250‑1 (MQ=255)
                                           aGTTTGGATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACGGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCGGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCGTTTTGCAGCTCGCGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCAGGACCGGAATCATGGc                                                                                                                                                                                                   <  2:437765/251‑1 (MQ=255)
                                           aGTTTGGATAACGCTTTCATTTCACATCATCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGc                                                                                                                                                                                                   >  2:29643/1‑251 (MQ=255)
                                            gTTTGGATAACGCTTTCATTTCACATCCTCATGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCGTTTTGCAGCTCACGAGCAAATTTCTGCGCCACGACCGCATCGCCAAGCATGACCGGAATAATGGCg                                                                                                                                                                                                  <  2:397744/251‑1 (MQ=255)
                                                ggATAACGCTTTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATACTGGCGTGAt                                                                                                                                                                                              >  1:346802/1‑251 (MQ=255)
                                                 gATAACGCTCTCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATc                                                                                                                                                                                             <  2:188864/251‑1 (MQ=255)
                                                          ttCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATCCGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTACATGCCCTCTTTTTGCAGATCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGc                                                                                                                                                                                    >  1:25521/1‑251 (MQ=255)
                                                          ttCATTTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGc                                                                                                                                                                                    >  2:354155/1‑251 (MQ=255)
                                                                  acaTCCTCAGGCGATAACGCCCCGTTGTTTACCAATACGCGGGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTaa                                                                                                                                                                            <  1:435462/251‑1 (MQ=255)
                                                                  acaTCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTa                                                                                                                                                                             >  1:94424/1‑250 (MQ=255)
                                                                              gATAACGCCCAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCGAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGTGTACAGCccgccgccg                                                                                                                                                                >  2:559/1‑251 (MQ=255)
                                                                                    acccAGTTGTTTGCCACTACGCGAGAATGCCTTTACCGCGCGCGTAAATTGCTCAGGGGTATGAGCCGCAGCCAGATGGGTACGAATACGCGGCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCACTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCAAAGCATGACCGGAATAATGGCGAGAACGGCTACAGCCAGGGGCAAGACCGCCGCCGAAAttt                                                                                                                                                          >  2:51157/2‑251 (MQ=255)
                                                                                     cccAGTTGTTTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTg                                                                                                                                                         <  2:351699/251‑1 (MQ=255)
                                                                                            gttTACCAATACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCGTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCTCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACgg                                                                                                                                                  <  1:111664/251‑1 (MQ=255)
                                                                                               ttccAATCCGCGTGAATGCTTCTACTGCGCGCTTAAGTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAACCCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCTTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGaa                                                                                                                                                <  2:351726/248‑1 (MQ=255)
                                                                                                    aTACGCGTGAATGCTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGAcg                                                                                                                                          >  1:424432/1‑251 (MQ=255)
                                                                                                                 cTTCTACTGCGCGCGTAATTTGCTCAGGGGTATGAGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCcac                                                                                                                             >  1:126536/1‑251 (MQ=255)
                                                                                                                     tACTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGAc                                                                                                                         >  1:371427/1‑251 (MQ=255)
                                                                                                                      aCTGCGCGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAACGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGGTCGCCCACCGACg                                                                                                                        >  2:358402/1‑251 (MQ=255)
                                                                                                                                 aaTTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTAACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAgct                                                                                                             >  1:310368/1‑249 (MQ=255)
                                                                                                                                  aTTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCAGCGCCAAGCATGACCGGAATAATGGCGGGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAgctc                                                                                                            <  1:78469/251‑4 (MQ=255)
                                                                                                                                         ccgggtttttcatctcAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCACCACTACCGCATCGCCAAGCATGACCTGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCCCGCAGCTCGCTGccc                                                                                                     >  1:370345/16‑251 (MQ=255)
                                                                                                                                                 tGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGATGCCCGCTTcaac                                                                                             >  1:295961/1‑248 (MQ=255)
                                                                                                                                                       gcAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCCCGCAGCTCGCTGCCCGCTTCCACCAtctc                                                                                       >  1:297560/1‑251 (MQ=255)
                                                                                                                                                                tGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACttt                                                                              >  1:232240/1‑251 (MQ=255)
                                                                                                                                                                       gAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGa                                                                        >  1:184645/1‑250 (MQ=255)
                                                                                                                                                                       gAATACGCGCCTGACCTCTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACAGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGATGCCCGCTTCAACCCTCTCCAGCACTTTAATGGAc                                                                       >  2:224452/1‑251 (MQ=255)
                                                                                                                                                                         aTACGCGCCTGACCTTTCGGAACGCCCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTTCAGCTCACTATCAAATTTCTGCGCCACTACCGCATCGCCAATCATGACCGGAATAATTGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGaagc                                                                     <  1:280897/251‑4 (MQ=255)
                                                                                                                                                                             gcgcCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGcg                                                                 <  2:355596/251‑2 (MQ=255)
                                                                                                                                                                                     acTTTCGGAACGACCGGATAAAAGAAACCGGTGACGCAAATGCCCTCTTTTTGCAACTCACGAGCAAATTTCTGCGCCACTACCGCCTCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGATCGCTGCCCGCTTCAACCATCTCCAGCCCTTTAATGGAAACGGcaacaaaagc                                                         >  2:166194/2‑247 (MQ=255)
                                                                                                                                                                                      cTTTCGGAACGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGcc                                                        >  1:305663/1‑251 (MQ=255)
                                                                                                                                                                                             aaCGACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGcc                                                  >  1:160718/1‑250 (MQ=255)
                                                                                                                                                                                                gACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTg                                                                                                        >  1:395978/1‑193 (MQ=255)
                                                                                                                                                                                                gACCGGATAGAAGAAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTGCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAgcgcgctg                                                                                                        <  2:395978/193‑8 (MQ=255)
                                                                                                                                                                                                        agaagaAACCGGTGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCCCCCCTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGACCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGga                                       <  2:424432/250‑1 (MQ=255)
                                                                                                                                                                                                        agaagaAACCGGTGACGGAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGgag                                      <  1:358402/251‑1 (MQ=255)
                                                                                                                                                                                                         gaagaaACCGTGGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCCCTACCGCACCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGgaga                                     <  2:305744/251‑1 (MQ=255)
                                                                                                                                                                                                             aaaCCGGCGACGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCCCGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGgagcaccg                                 >  2:186026/1‑246 (MQ=255)
                                                                                                                                                                                                                       cGTAAATGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGGAGAACCGGTACgggcg                        >  2:36878/1‑247 (MQ=255)
                                                                                                                                                                                                                          aaaTGCCCTCTTTTTGCAGCTCACGAGCAAATTTCTGCGCCACTAACGCATCGCCAAGCATGACCGGAAGAATGGCGTGATCGGCACCCGCCAGGGAAAAGCACGCCGCCGACATTGGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGGGACAAGAGCCGGTGCAAGCGAGTTGGAGAACAGGGACGGACgaga                     >  1:143653/1‑250 (MQ=255)
                                                                                                                                                                                                                                cgtccttttGCAGCTACCGAGCAATTTTCTGCGCCACCACCGCACCTCCAAGCATGACCGGAATAATGGCGCGATCGGCTTCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTACGCCCACAGACGGTCACGCAGCGCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGGGCCAGCTAGTTGGAGAACAGGTACGGACGAGAACGCTgg              <  1:14101/246‑1 (MQ=255)
                                                                                                                                                                                                                                       ttGCAGCTCACGAGCAAATTTCTGCGCCACTCCCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGc        <  1:317588/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gCTCACGAGCAAATTTCTGCGCCACTACCGCATCGCCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGCTCGCTGCCCGCTTCAACCATCTCCAGCACTTTAATGGAAGCGGCGACAATAGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCa  <  1:121210/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CGGTACAGGAACGTCGGTCATCCAGATGCCCTCTTCCGCTTTCAGTTTGGATAACGCTTTCATCTCACATCCTCAGGCGATAACGCCCAGTTGTTTACCAATACGCGTAAATGCTTCTACTGCACGCGTAATTTGCTCAGGGGTATGCGCCGCAGACATCTGGGTACGAATACGCGCCTGACCTTTCGGAACGACCGGATAGAAGAAACCGGTAACGTAAATGCCCTCTTTTTGCAGCTCACGGGCAAATTTCTGCGCCACTACCGCATCACCAAGCATGACCGGAATAATGGCGTGATCGGCTCCCGCCAGGGTAAAGCCCGCCGCCGACATTTGCTCACGGAACTGACGCGCGTTCGCCCACAGACGGTCACGCAGTTCGCTGCCCGCTTCGACCATCTCCAGTACTTTGATGGACGCGGCAACAATGGCCGGTGCCAGCGAGTTGGAGAACAGGTACGGACGAGAACGCTGGCGCAGCCACTCA  >  NC_000913/3791283‑3791769

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