Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 984,677 T→A 6.3% N345Y (AAC→TAC)  aspC ← aspartate aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913984,6770TA 6.3% 192.9 / 3.7 64N345Y (AAC→TAC) aspCaspartate aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

AATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAGGATGGT  >  NC_000913/984534‑984824
                                                                                                                                               |                                                                                                                                                   
aaTCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGt                                                                                                                                                 <  6:154480/148‑1 (MQ=255)
 aTCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCCTCCCGGCCACAGTTACGCGACAAGAAGCACCCGAATATACGCCAAACTCTTCGCTCAGACGCAGCACTTGTTCTTTTGTCAGGCAAGTGAAGGAGACCATTCCGTTCTGtt                                                                                                                                                >  2:248247/1‑148 (MQ=255)
 aTCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTg                                                                                                                                                  >  5:42727/1‑146 (MQ=255)
 aTCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTg                                                                                                                                                  <  6:42727/146‑1 (MQ=255)
 aTCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGtt                                                                                                                                                >  7:213172/1‑148 (MQ=255)
    gCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTtg                                                                                                                                              <  3:266041/147‑1 (MQ=255)
     cTTCGCACAGCGGAGCCATGTTATCTGCAGTAATCCCGGCCACATTTACGCGAACAGAAGCAACCGCATAAACGCGAAACTCTTCGCGCAGACGCAGCACTTGTAAATTTCACAGGCCACTGAAGGAGAACATGCCGTTCTCTTtg                                                                                                                                              >  5:189206/1‑146 (MQ=255)
                ggAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCTCATATAGGCCAACCTCTTCGCGCAGACGCATCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGa                                                                                                                                 <  1:47839/148‑1 (MQ=255)
                   gCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACAGGCCGTTCTGTTTGATGATAAAGCTGAAGt                                                                                                                              >  4:144237/1‑148 (MQ=255)
                      aTGTTATCTGTGGGCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTcgc                                                                                                                           <  8:139193/148‑1 (MQ=255)
                        gTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCgg                                                                                                                         >  7:18568/1‑148 (MQ=255)
                                            ccACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAg                                                                                                     <  3:286562/148‑1 (MQ=255)
                                                aTTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCGTGCAGCGTa                                                                                                 >  6:61518/1‑148 (MQ=255)
                                                aTTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCC‑TTTTCCTGCAGCGTAt                                                                                                >  3:147413/1‑148 (MQ=255)
                                                    aCGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTg                                                                                              <  7:254879/147‑1 (MQ=255)
                                                     cgcgACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGTAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGa                                                                                             <  7:114730/147‑1 (MQ=255)
                                                     cgcgACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGa                                                                                             <  8:286788/147‑1 (MQ=255)
                                                         aCCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGaa                                                                                         >  2:173650/1‑147 (MQ=255)
                                                         aCCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGaa                                                                                         >  5:313034/1‑147 (MQ=255)
                                                         aCCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTACTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGaa                                                                                         >  6:248186/1‑147 (MQ=255)
                                                         aCCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTACTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGaa                                                                                         >  6:242718/1‑147 (MQ=255)
                                                               aGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGa                                                                                  >  6:158689/1‑148 (MQ=255)
                                                               aGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGa                                                                                  <  1:322240/148‑1 (MQ=255)
                                                                 cAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACg                                                                                >  2:170558/1‑148 (MQ=255)
                                                                 cAACCGCATATACGCCAAACTCTTCGCGCAGACGAAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACg                                                                                <  7:383628/148‑1 (MQ=255)
                                                                  aaCCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGc                                                                               <  8:213172/148‑1 (MQ=255)
                                                                   aCCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTATTTCCTGCAGCGTATTGACGAACAACTGACGCa                                                                              >  5:157760/1‑148 (MQ=255)
                                                                     cGCATATACGCCAAACTCTTCGCGCAGCCGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATa                                                                            >  5:4287/1‑148 (MQ=255)
                                                                          ataCGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCt                                                                        >  4:177808/1‑147 (MQ=255)
                                                                           taCGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTg                                                                       <  8:18568/147‑1 (MQ=255)
                                                                           taCGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGa                                                                      <  4:103922/148‑1 (MQ=255)
                                                                                cAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACg                                                                 <  5:158689/148‑1 (MQ=255)
                                                                                      ctTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATATAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTg                                                                       <  5:53124/136‑1 (MQ=255)
                                                                                      ctTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATATAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTg                                                                       >  6:53124/1‑136 (MQ=255)
                                                                                            gcAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGATAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtc                                                     <  6:303393/148‑1 (MQ=255)
                                                                                            gcAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtc                                                     <  3:340696/148‑1 (MQ=255)
                                                                                                aCGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtcagt                                                  >  5:286807/1‑147 (MQ=255)
                                                                                                 cGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTACTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtcagtca                                                <  5:242718/148‑1 (MQ=255)
                                                                                                 cGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTACTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtcagtca                                                <  5:248186/148‑1 (MQ=255)
                                                                                                  gcagcaCTTGTTCTTTTGTCAGGCCACTGATGGAGAACATGCCGTTGTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGATCAACTGACGCATACGCTGAATTCGCTGGCGCATAtcagtcag                                               <  6:189206/148‑1 (MQ=255)
                                                                                                  gcagcaCTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATAtcagtcag                                               <  8:284027/148‑1 (MQ=255)
                                                                                                   cagcaCTTGTTCTTTTGTCAGGCCACTTAATGATAACATGCCTTTCTGTTTTATGATAAATATTAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACTCTGGCGCATATCAGTCAGc                                              <  4:126965/148‑1 (MQ=255)
                                                                                                   cagcaCTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGTTTAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGc                                              <  4:121714/148‑1 (MQ=255)
                                                                                                   cagcaCTTCTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATAAGCTGAAAACGCTGGCGCATATCAGTCAGc                                              >  2:49891/1‑148 (MQ=255)
                                                                                                      caCTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGGTGGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGctct                                           <  2:20032/148‑1 (MQ=255)
                                                                                                        cTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTAt                                                                                                <  3:166256/93‑1 (MQ=255)
                                                                                                        cTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTAt                                                                                                >  4:166256/1‑93 (MQ=255)
                                                                                                        cTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTg                                         >  4:136882/1‑148 (MQ=255)
                                                                                                              cTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTccc                                    >  5:150774/1‑147 (MQ=255)
                                                                                                                tttGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCaaa                                 >  6:321249/1‑148 (MQ=255)
                                                                                                                tttGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCCGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCaaa                                 >  5:63590/1‑148 (MQ=255)
                                                                                                                    tCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCg                              >  4:174559/1‑147 (MQ=255)
                                                                                                                    tCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGc                             <  3:174559/148‑1 (MQ=255)
                                                                                                                       ggCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACg                          >  1:109632/1‑148 (MQ=255)
                                                                                                                           aCTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAAc                      >  3:184853/1‑148 (MQ=255)
                                                                                                                                     aaCATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTcc                                     <  2:294070/123‑1 (MQ=255)
                                                                                                                                     aaCATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTcc                                     >  1:294070/1‑123 (MQ=255)
                                                                                                                                       cATGCCGTTCTGTTTGATGATAAAGCTGAAGGCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTc          <  5:41989/148‑1 (MQ=255)
                                                                                                                                         tGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAg        <  6:150774/148‑1 (MQ=255)
                                                                                                                                          gCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTACAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAg        >  5:127552/1‑147 (MQ=255)
                                                                                                                                           ccGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAGGa      >  2:2405/1‑148 (MQ=255)
                                                                                                                                             gTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTc          >  3:175625/1‑142 (MQ=255)
                                                                                                                                             gTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTc          <  4:175625/142‑1 (MQ=255)
                                                                                                                                             gTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAGGAtg    >  6:373647/1‑148 (MQ=255)
                                                                                                                                               tCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAGGAtggt  >  1:271969/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AATCGCTTCGCACAGCGGAGCCATGTTATCTGGTGTCATCCCGGCCACATTTACGCGACCAGAAGCAACCGCATATACGCCAAACTCTTCGCGCAGACGCAGCACTTGTTCTTTTGTCAGGCCACTGAAGGAGAACATGCCGTTCTGTTTGATGATAAAGCTGAAGTCGCGGTTTGCGCCTTTTTCCTGCAGCGTATTGACGAACAACTGACGCATACGCTGAATACGCTGGCGCATATCAGTCAGCTCTTGTTCCCAAATCGCACGTAACGCATCGTTGCTCAGGATGGT  >  NC_000913/984534‑984824

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