Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,626,189 Δ1 bp 6.0% coding (1192/2046 nt) dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,626,1890A. 6.0% 205.3 / 6.1 67coding (1192/2046 nt)dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base . (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

ACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAG  >  NC_000913/1626044‑1626330
                                                                                                                                                 |                                                                                                                                             
aCGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCATTTTGAATCACCGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCcaca                                                                                                                                             <  3:377882/148‑1 (MQ=255)
      taAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGaa                                                                                                                                       >  4:37228/1‑148 (MQ=255)
      taAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGa                                                                                                                                        <  7:291726/147‑1 (MQ=255)
                gATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATg                                                                                                                             <  8:317466/148‑1 (MQ=255)
                  tgtgTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCAACGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATgg                                                                                                                            >  1:370114/1‑147 (MQ=255)
                    tgtTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTa                                                                                                                          <  3:37228/147‑1 (MQ=255)
                      ttttATTAAGCGTATATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCATCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGAc                                                                                                                       <  4:216676/148‑1 (MQ=255)
                       tttATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCCCGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGa                                                                                                                        >  1:189387/1‑146 (MQ=255)
                       tttATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCca                                                                                                                                               >  2:427037/1‑123 (MQ=255)
                       tttATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCca                                                                                                                                               >  2:420880/1‑123 (MQ=255)
                       tttATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCca                                                                                                                                               <  1:427037/123‑1 (MQ=255)
                       tttATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCca                                                                                                                                               <  1:420880/123‑1 (MQ=255)
                        ttattaAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACATGAACGTCAGGATGGTAGACa                                                                                                                      <  1:231078/147‑1 (MQ=255)
                              aGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTTAAATACCGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACCCGAACGTCAGGATGGTAGACAGGAatat                                                                                                               <  5:428157/148‑1 (MQ=255)
                                   tGATTGCTCAACAAAATTGCCCATCCATGCACCTACGCTTTTTGAATCACGTGCGAAGAAATCACCGTAAAATAACGCCAGCACCACGCCATTATGATCAAAAATTTACCACACACGAACTTCAGGATGGTAGACAGGAATATCaaaa                                                                                                          >  7:208156/1‑148 (MQ=255)
                                           caacaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTc                                                                                                                                   >  4:111830/1‑115 (MQ=255)
                                           caacaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTc                                                                                                                                   <  3:111830/115‑1 (MQ=255)
                                             acaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCaaaa                                                                                                          <  2:317783/138‑1 (MQ=255)
                                             acaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCaaaa                                                                                                          >  1:317783/1‑138 (MQ=255)
                                             acaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTc                                                                                                     <  7:45335/143‑1 (MQ=255)
                                             acaaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTc                                                                                                     >  8:45335/1‑143 (MQ=255)
                                              caaAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACaaa                                                                                               <  6:411754/148‑1 (MQ=255)
                                               aaaaTTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGt                                                                                                       <  3:228526/139‑1 (MQ=255)
                                               aaaaTTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGt                                                                                                       >  4:228526/1‑139 (MQ=255)
                                                aaaTTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGAAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGCAATATCAAACCGTTCGACAAACt                                                                                             >  7:180158/1‑148 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATTATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATa                                                                                        <  5:165102/142‑1 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGa                                                                                                                                        >  7:206337/1‑94 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGa                                                                                                                                        <  8:206337/94‑1 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGt                                                                                                                           <  2:49823/107‑1 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGt                                                                                                                           >  1:49823/1‑107 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATa                                                                                        >  6:165102/1‑142 (MQ=255)
                                                           ccatGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCg                                                                                     >  7:316045/1‑145 (MQ=255)
                                                             attcACCGCCGCGTATTGAATCACGGGCGAATAAATCACCGTAAAATAACTACATACCCAATACATTATGATCTAAAATTTCCCACACAAGAACGTCAGGATGGTTGACAGGAATATCAAAACGTTCTACAAACTTAATACCGAatag                                                                                <  8:208156/145‑1 (MQ=255)
                                                                    gccgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTc                                                                         >  4:135371/1‑148 (MQ=255)
                                                                    gccgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTc                                                                         >  7:400837/1‑148 (MQ=255)
                                                                     ccgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTcg                                                                        <  3:135371/148‑1 (MQ=255)
                                                                      cgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTa                                                                                           <  2:122862/128‑1 (MQ=255)
                                                                      cgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTa                                                                                           >  1:122862/1‑128 (MQ=255)
                                                                      cgcTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTcgc                                                                       >  1:72874/1‑148 (MQ=255)
                                                                          ttttGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACt                                                                                             >  2:322104/1‑122 (MQ=255)
                                                                          ttttGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACt                                                                                             <  1:322104/122‑1 (MQ=255)
                                                                                 tCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGa                                                                                    >  8:249478/1‑124 (MQ=255)
                                                                                 tCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGa                                                                                    <  7:249478/124‑1 (MQ=255)
                                                                                   aCGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATaca                                                          <  8:400837/148‑1 (MQ=255)
                                                                                   aCGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATac                                                           <  8:316045/147‑1 (MQ=255)
                                                                                            gaaATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCGAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCAtt                                                  >  5:224559/1‑147 (MQ=255)
                                                                                              aaTCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTa                                                <  2:370114/147‑1 (MQ=255)
                                                                                                     gTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCgtg                                        >  4:174712/1‑148 (MQ=255)
                                                                                                     gTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCC‑‑CACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGtt                                      >  2:434447/1‑148 (MQ=255)
                                                                                                     gTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCC‑‑CACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGtt                                      >  2:440670/1‑148 (MQ=255)
                                                                                                        aaaTAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGtt                                      <  3:174712/147‑1 (MQ=255)
                                                                                                         aaTAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTa                                    <  2:189387/148‑1 (MQ=255)
                                                                                                          aTAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTcc                                                                  >  7:365458/1‑117 (MQ=255)
                                                                                                          aTAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTcc                                                                  <  8:365458/117‑1 (MQ=255)
                                                                                                                    ccccACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGa                                                                            >  2:349792/1‑97 (MQ=255)
                                                                                                                    ccccACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGa                                                                            <  1:349792/97‑1 (MQ=255)
                                                                                                                        aCGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCtt                                                      >  6:18176/1‑115 (MQ=255)
                                                                                                                        aCGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCtt                                                      <  5:18176/115‑1 (MQ=255)
                                                                                                                        aCGCCATTATGATCAAAAATTTCC‑‑CACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTg                   <  1:434447/148‑1 (MQ=255)
                                                                                                                        aCGCCATTATGATCAAAAATTTCC‑‑CACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTg                   <  1:440670/148‑1 (MQ=255)
                                                                                                                               tATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGttt                                     <  8:101858/125‑1 (MQ=255)
                                                                                                                               tATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGttt                                     >  7:101858/1‑125 (MQ=255)
                                                                                                                               tATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTg              >  6:120869/1‑148 (MQ=255)
                                                                                                                                       aaaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATaca                                                          >  7:189835/1‑96 (MQ=255)
                                                                                                                                       aaaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATaca                                                          <  8:189835/96‑1 (MQ=255)
                                                                                                                                       aaaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCaaa                             >  4:88013/1‑125 (MQ=255)
                                                                                                                                       aaaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCaaa                             <  3:88013/125‑1 (MQ=255)
                                                                                                                                        aaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTcgc                                                                       <  6:184639/82‑1 (MQ=255)
                                                                                                                                        aaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTcgc                                                                       >  5:184639/1‑82 (MQ=255)
                                                                                                                                        aaaTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCtcat      >  4:411807/1‑147 (MQ=255)
                                                                                                                                          aTTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCtcatca    >  2:432297/1‑147 (MQ=255)
                                                                                                                                           tttCCAACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAg  <  5:120869/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
ACGTTATAAATTACCGGATGTGTTTTATTAAGCGTTGATTGCTCAACAAAATTGCCCATCCATGCACCGCCGCTTTTTGAATCACGGGCGAAGAAATCACCGTAAAATAACGCCAGCCCCACGCCATTATGATCAAAAATTTCCCACACACGAACGTCAGGATGGTAGACAGGAATATCAAAACGTTCGACAAACTTAATACCGAAGAGCTGATTCGCGGTCCAGAATACACCTTCATTTAACACCGTGTTTAATTCAAAATATGGCTTGAGCTGCGCCTCATCAAG  >  NC_000913/1626044‑1626330

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: