Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,460,345 C→T A431V (GCC→GTC)  paaH → 3‑hydroxyadipyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,460,3450CT95.9% 147.0 / ‑1.3 52A431V (GCC→GTC) paaH3‑hydroxyadipyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (26/21);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

GTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAC  >  NC_000913/1460104‑1460586
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
gtggtgACCCTTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCACAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGATTTATCCATATGGCCCACTTGTCTGGGGAgc                                                                                                                                                                                                                                          >  2:229798/1‑251 (MQ=255)
  ggtgACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAgcg                                                                                                                                                                                                                                         <  1:170569/250‑1 (MQ=255)
  ggtgACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGACTGTGGAgcgc                                                                                                                                                                                                                                        >  1:453782/1‑251 (MQ=255)
       ccATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACtt                                                                                                                                                                                                                                   <  1:192119/251‑1 (MQ=255)
         aTTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGGGAATTATCCATATGGCCCACTTGTCTGGGGGGCGCAACTtgg                                                                                                                                                                                                                                 >  1:102848/1‑251 (MQ=255)
                    aGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGaa                                                                                                                                                                                                                      <  1:8701/251‑1 (MQ=255)
                              ccGGACTAAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCGGGCAGCGAATATGAAGGCt                                                                                                                                                                                                            >  2:117955/1‑251 (MQ=255)
                                ggACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTcc                                                                                                                                                                                                          <  1:263765/251‑1 (MQ=255)
                                     cAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTg                                                                                                                                                                                                       <  1:344032/249‑1 (MQ=255)
                                        cGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGAAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAAt                                                                                                                                                                                                  <  1:386971/251‑1 (MQ=255)
                                         gACCCGCAAGGCCATGTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGCATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGAAATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCACCTTGTCTGGGGGGCGCAACTTGGCTGGAAGCGAATATTAAGGATGATTGAAAATc                                                                                                                                                                                                 >  2:359347/1‑251 (MQ=255)
                                            ccGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTAc                                                                                                                                                                                              <  1:437048/251‑1 (MQ=255)
                                                    ccATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCAc                                                                                                                                                                                      <  1:312868/251‑1 (MQ=255)
                                                        ttattaCCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATg                                                                                                                                                                                  >  1:433507/1‑251 (MQ=255)
                                                            taCCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCg                                                                                                                                                                               >  1:198124/1‑250 (MQ=255)
                                                                   aaCAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGc                                                                                                                                                                       >  1:314407/1‑251 (MQ=255)
                                                                                    gggctgcAAAGTGCATATTACCCAGGAATGCTGATTGGGAGAACGGTAGCAAGGATCATCAATGATGCCATTGATGCGCTTCAAAAAGCCGTGGCCTCGTAACAGGATTTCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCATCTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACTCTATCGCCCATGTTCAtt                                                                                                                                                      <  2:167263/249‑1 (MQ=255)
                                                                                           aaaTTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTgcgc                                                                                                                                               <  1:103599/251‑1 (MQ=255)
                                                                                             aTTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGcc                                                                                                                                              >  1:314192/1‑250 (MQ=255)
                                                                                                     ttACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTtggct                                                                                                                                                                                                                               <  2:235633/161‑1 (MQ=255)
                                                                                                     ttACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTtggct                                                                                                                                                                                                                               >  1:235633/1‑161 (MQ=255)
                                                                                                     ttACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCt                                                                                                                                                                      >  1:421800/1‑218 (MQ=255)
                                                                                                     ttACCCAGGAATGCTGATTTGGCGAACGGTAGCAACGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCt                                                                                                                                                                      <  2:421800/218‑1 (MQ=255)
                                                                                                        cccAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTc                                                                                                                                  >  1:484326/1‑251 (MQ=255)
                                                                                                           aGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGCTACCGCCATGCGTCTTGGGGTGAATTATACATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAACGCATATTAAGGCTCCTTGATAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCCTTGCTGCGCCAACGGGAGCTTCTgg                                                                                                                               >  2:143420/1‑251 (MQ=255)
                                                                                                                   tGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGACCAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTACGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTgcg                                                                                                                                                >  2:427884/1‑226 (MQ=255)
                                                                                                                   tGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTgcg                                                                                                                                                <  1:427884/226‑1 (MQ=255)
                                                                                                                    gATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGtt                                                                                                                      >  1:353776/1‑251 (MQ=255)
                                                                                                                     aTTTGGCGAAGGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGATGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTAATTACTTCGCCAACGGGCGCTTCTGGAGAGCggat                                                                                                                      >  2:275307/1‑248 (MQ=255)
                                                                                                                           cGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTc                                                                                                               <  1:533265/251‑1 (MQ=255)
                                                                                                                            gAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGAGCCAACGGGCGCTTCTGGAGAGCGGTTGTGAGTc                                                                                                               >  2:488016/1‑250 (MQ=255)
                                                                                                                                      aaTGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTgg                                                                                                    >  1:157829/1‑251 (MQ=255)
                                                                                                                                              ccaATGAATGCCTTGATGCGCGTAAAAAAGGCATGGACTCGTAACATGATAGCGGTGCCGCCATGCGACGTGGGGTGAATTATCCATATGGCCCACTTGGCTGGGGAGCGCAACGTGGCGTGCAGCTAATATGAAGGCTCCTTGAAAATCTTCAGCATCACTTTGGCGAAGAACGCTATCGCCCATGTTCATTGCTGTGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGc                                                                                            <  2:531930/250‑1 (MQ=255)
                                                                                                                                                aaTGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGccc                                                                                          >  2:92379/1‑251 (MQ=255)
                                                                                                                                                  tGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGACCGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGCGTACTAAGGCCTGGCACAATGcccc                                                                                         >  2:23706/1‑249 (MQ=255)
                                                                                                                                                    aggcccttgatagtccgaaaaaaGGCGTGGCGTCAGATCAGGATATATATGCCACACTGCGTCTGGTGGTGAATTAGCCATAGTGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGTAGCGAATTTTAAGGCTCCTTGAAAATCTACAGCATCACTATGGTGAAGAACGCTATCGCGCAGGTTCATTGCTGCGCCTACGTGCGCTTCTGGAGAGCTGTTATGAGTCATAAGGCCTGGCAAAATGCCCAttc                                                                                      <  2:69527/233‑3 (MQ=255)
                                                                                                                                                    aaGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGc                                                                                      >  1:315272/1‑251 (MQ=255)
                                                                                                                                                        ccTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGACATGGGGTGAATTATCCATATGGCCAACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGACAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTACTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAATGCCTGGCAAAAAGACCATGCAATg                                                                                  >  2:40963/1‑251 (MQ=255)
                                                                                                                                                               gcgcTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATgaga                                                                           <  1:457325/251‑1 (MQ=255)
                                                                                                                                                                cgcTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGaa                                                                          <  1:92379/251‑1 (MQ=255)
                                                                                                                                                                 gcTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAAc                                                                         >  1:242026/1‑251 (MQ=255)
                                                                                                                                                                      aaaaaGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGCGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGc                                                                    >  1:289361/1‑251 (MQ=255)
                                                                                                                                                                                ggCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTAACGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGccca                                                           >  2:128307/1‑248 (MQ=255)
                                                                                                                                                                                       gAACAGGATATCGATACCGCCAGGTGTCTTGGGGGGAATTATCCATATGGCCCACTTGTCTGGGTAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAgcgc                                                      <  2:314192/248‑1 (MQ=255)
                                                                                                                                                                                        aaCAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGGCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGATTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGACTGCGCCAAAGCGCTTga                                                  <  2:57005/251‑2 (MQ=255)
                                                                                                                                                                                             gATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCg                                             >  1:230296/1‑251 (MQ=255)
                                                                                                                                                                                                atCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGAc                                           >  1:471892/1‑250 (MQ=255)
                                                                                                                                                                                                                ggcgtGGGGTGAATTATTCAATTGGCCCCCTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACTGGCGCTTCTGGAGAGCGGTTATGAGTCATTAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGACTGCGCCAAAGCGCTTGGCATTGACATTATCTCAATGGATGa                          <  2:102848/247‑1 (MQ=255)
                                                                                                                                                                                                                gTCTTGGGGTGAATTATCCATATGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATAGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATg                           >  1:31802/1‑250 (MQ=255)
                                                                                                                                                                                                                            aTTATCCATATGGCCCTCTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCCCCTTGAAAATCTACAGCATCACTATGTCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCGGGCAAAATGCCCATGCAATGTATGAGATCGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGt              <  2:315272/251‑1 (MQ=255)
                                                                                                                                                                                                                                     atGGCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCag     <  2:198124/251‑2 (MQ=255)
                                                                                                                                                                                                                                        gCCCACTTGTCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGGCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAc  <  1:524700/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GTGGTGACCATTGCTGCTGCAGCGGTGAACCCGGACTCAGCGACCCGCAAGGCCATTTATTACCTGCAACAGCAGGGCAAAACAGTGCTGCAAATTGCAGATTACCCAGGAATGCTGATTTGGCGAACGGTAGCAATGATCATCAATGAAGCCCTTGATGCGCTTCAAAAAGGCGTGGCCTCTGAACAGGATATCGATACCGCCATGCGTCTTGGGGTGAATTATCCATATGGCCCACTTGCCTGGGGAGCGCAACTTGGCTGGCAGCGAATATTAAGGCTCCTTGAAAATCTACAGCATCACTATGGCGAAGAACGCTATCGCCCATGTTCATTGCTGCGCCAACGGGCGCTTCTGGAGAGCGGTTATGAGTCATAAGGCCTGGCAAAATGCCCATGCAATGTATGAGAACGATGCCTGCGCCAAAGCGCTTGGCATCGACATTATCTCAATGGATGAAGGCTTTGCTGTAGTGACCATGAC  >  NC_000913/1460104‑1460586

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