Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 111,433 G→A *130* (TAG→TAA mutT → dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,4330GA98.7% 237.4 / ‑4.4 76*130* (TAG→TAAmutTdGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/37);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

GCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAC  >  NC_000913/111298‑111567
                                                                                                                                       |                                                                                                                                      
gCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTc                                                                                                                                     >  1:733224/1‑139 (MQ=255)
       gAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAg                                                                                                                              <  1:902732/139‑1 (MQ=255)
               ggAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTc                                                                                                                      >  1:1650885/1‑139 (MQ=255)
               ggAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTAAAACGTCTGTAGGTCAGATAAGGCGTTTTc                                                                                                                      >  2:391841/1‑139 (MQ=255)
                 aaGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGAt                                                                                                                                 >  1:1275842/1‑126 (MQ=255)
                 aaGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGAt                                                                                                                                 <  2:1275842/126‑1 (MQ=255)
                       aGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               >  2:453961/1‑122 (MQ=255)
                       aGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               <  1:453961/122‑1 (MQ=255)
                                        gggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGAt                                                                                             <  1:260361/139‑1 (MQ=255)
                                         ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGt                                                                                                                                      <  2:1586372/97‑1 (MQ=255)
                                         ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGt                                                                                                                                      >  1:1586372/1‑97 (MQ=255)
                                         ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    <  1:668708/115‑1 (MQ=255)
                                         ggCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    >  2:668708/1‑115 (MQ=255)
                                           cAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                >  2:1604007/1‑133 (MQ=255)
                                           cAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                <  1:1604007/133‑1 (MQ=255)
                                               ccGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATcc                                                                                                             >  2:1503241/1‑116 (MQ=255)
                                               ccGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATcc                                                                                                             <  1:1503241/116‑1 (MQ=255)
                                                cGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATg                                                                                     <  2:733224/139‑1 (MQ=255)
                                                  gTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCg                                                                                   >  1:682163/1‑139 (MQ=255)
                                                     aGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACg                                                                                >  1:51097/1‑139 (MQ=255)
                                                        ggATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTg                                                                             >  1:295641/1‑139 (MQ=255)
                                                         gATGTCGCTGGTCGGTCTTAATGCAGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGATTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  1:159835/139‑1 (MQ=255)
                                                          aTGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGgc                                                                           >  2:505372/1‑139 (MQ=255)
                                                             tCGCTGGTCGGTCTTAATGCCGATGAGTTTCCGCGAGCCAATGAGCCGGTAATTGCCAAGCTTAAACGTCATTAAGTCAGTTATTCAGTTTTCGCCGCATCCGCCATTCGCACACGATACCTTGTGGGAGGCTGGcgcg                                                                        >  2:159835/1‑139 (MQ=255)
                                                              cGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                                       <  2:723757/139‑1 (MQ=255)
                                                              cGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGt                                                                       <  2:682163/139‑1 (MQ=255)
                                                               gCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    >  2:1098599/1‑93 (MQ=255)
                                                               gCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    <  1:1098599/93‑1 (MQ=255)
                                                               gCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTc                                                                      >  2:1116893/1‑139 (MQ=255)
                                                                 tGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCtt                                                                    >  2:67842/1‑139 (MQ=255)
                                                                   gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  <  2:1650885/139‑1 (MQ=255)
                                                                   gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  >  1:1222675/1‑139 (MQ=255)
                                                                   gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  <  1:1116893/139‑1 (MQ=255)
                                                                   gtcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  <  1:67842/139‑1 (MQ=255)
                                                                    tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  2:1167875/128‑1 (MQ=255)
                                                                    tcggtcTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            >  1:1167875/1‑128 (MQ=255)
                                                                          ttAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt                                                           <  1:505372/139‑1 (MQ=255)
                                                                             aTGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               >  2:1525165/1‑68 (MQ=255)
                                                                             aTGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATaa                                                                                                                               <  1:1525165/68‑1 (MQ=255)
                                                                             aTGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCACATCCGACATTCGcaca                                                                                                 <  1:580971/98‑1 (MQ=255)
                                                                             aTGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCACATCCGACATTCGcaca                                                                                                 >  2:580971/1‑98 (MQ=255)
                                                                                ccGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAggg                                                     <  2:295641/139‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccg                                                                                                                  <  1:617242/76‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgccg                                                                                                                  >  2:617242/1‑76 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    >  2:1265620/1‑74 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTcgc                                                                                                                    <  1:1265620/74‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGca                                                                                                   >  2:354118/1‑91 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGca                                                                                                   <  1:354118/91‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                >  1:713878/1‑94 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCacac                                                                                                <  2:713878/94‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGTCTGATGCGACGCTgg                                                                            >  2:1066590/1‑114 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  1:1066590/114‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            >  2:837382/1‑114 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTgg                                                                            <  1:837382/114‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAgg                                                      <  1:1276251/136‑1 (MQ=255)
                                                                                  gatgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAgg                                                      >  2:1276251/1‑136 (MQ=255)
                                                                                   atgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCg                                                                                                     >  2:1010826/1‑88 (MQ=255)
                                                                                   atgatTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCg                                                                                                     <  1:1010826/88‑1 (MQ=255)
                                                                                             gccagccaATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCAt                                        >  1:437562/1‑139 (MQ=255)
                                                                                                            ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCg                                                                                                     >  1:79071/1‑63 (MQ=255)
                                                                                                            ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCg                                                                                                     <  2:79071/63‑1 (MQ=255)
                                                                                                            ggTAATTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGtttt                         >  1:1631117/1‑139 (MQ=255)
                                                                                                               aaTTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGt                      <  2:437562/139‑1 (MQ=255)
                                                                                                                aTTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                >  2:414437/1‑96 (MQ=255)
                                                                                                                aTTGCGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCa                                                                <  1:414437/96‑1 (MQ=255)
                                                                                                                    cGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGt                  <  1:468282/138‑1 (MQ=255)
                                                                                                                    cGAAGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGt                  >  2:468282/1‑138 (MQ=255)
                                                                                                                      aaGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCgg               >  1:285051/1‑139 (MQ=255)
                                                                                                                       aGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCt                                               <  2:506497/106‑1 (MQ=255)
                                                                                                                       aGCTTAAACGTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCt                                               >  1:506497/1‑106 (MQ=255)
                                                                                                                                gTCTGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGtt     <  2:891882/139‑1 (MQ=255)
                                                                                                                                   tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGAc                                                                                                          >  2:654052/1‑35 (MQ=14)
                                                                                                                                   tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGAc                                                                                                          <  1:654052/35‑1 (MQ=14)
                                                                                                                                   tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  >  1:318980/1‑75 (MQ=255)
                                                                                                                                   tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTAt                                                                  <  2:318980/75‑1 (MQ=255)
                                                                                                                                   tGTAAGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAc  >  1:1347838/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCAC  >  NC_000913/111298‑111567

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