Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 1,544,022 1 bp→AA coding (39/1389 nt) narU ← nitrate/nitrite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA1,544,0211.A82.3% 159.1 / 30.7 62coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (8/3);  new base A (31/20);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA1,544,0220CA96.8% 176.0 / NA 62L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (38/22);  total (39/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGA  >  NC_000913_3_pae_tpiA/1543879‑1544168
                                                                                                                                                  |||                                                                                                                                                  
ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCAAGTCTCGCaaaaa                                                                                                                                               >  1:768205/1‑148 (MQ=255)
   agaTTAACGGTAACTGCGCTAAATAGCATCCAGGCACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaagat                                                                                                                                            >  2:543823/1‑145 (MQ=255)
   agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaagat                                                                                                                                            >  2:564796/1‑145 (MQ=255)
           ggTAACTGCGCTAAATAGACTACAGACACAGAAGGCAACAACTCTGGTGCTGAACCATCTCAACTGAAGTTGTCTAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTATGGTTTCCAGTCGCAAAAAAAGATAACGACTa                                                                                                                                     >  2:220363/1‑149 (MQ=255)
                cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAttttt                                                                                                                                <  2:387942/149‑1 (MQ=255)
                  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAtttt                                                                                                                                 <  2:14959/146‑1 (MQ=255)
                  gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAtttt                                                                                                                                 >  1:14959/1‑146 (MQ=255)
                    gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCATTTATTTTACAAGAAGGCCGGATTTTCTGGTTTCGAGTCACGAAAAAAGATAAAGACTAtttttacta                                                                                                                            >  1:397820/1‑145 (MQ=255)
                          tAGCATCCAGACACAGAAGGCAAGAAGGAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGc                                                                                                                      >  2:552499/1‑149 (MQ=255)
                           aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                     <  1:663775/149‑1 (MQ=255)
                                ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCTCGAAAAAAGATAACGACTATTTTTCTAATTTTGCAGGGcc                                                                                                                >  2:70446/1‑149 (MQ=255)
                                      acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGGGCCATATTg                                                                                                          >  2:829908/1‑149 (MQ=255)
                                        aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                        >  2:738779/1‑149 (MQ=255)
                                        aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTATGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTTATTTTGCAGTTCCATATTGtt                                                                                                        >  2:1011560/1‑149 (MQ=255)
                                         gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                       <  1:543823/149‑1 (MQ=255)
                                         gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAACGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGTCTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                       <  2:990336/149‑1 (MQ=255)
                                         caaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCACTCGAGACAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                       <  2:527416/148‑1 (MQ=255)
                                             gCAAGGAGTAGACAACTGACCTATATCCAGAGACTTCTGCGAGCAATATGTTTTCCTTTATTTGCCCAGAAGGCCGGATTTGCTGTTTAACAGTAGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                   <  2:56790/149‑1 (MQ=255)
                                             gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATCACGACTATTTTTGCCAGCTTGCAGTGCCATATTGTTCCTca                                                                                                   >  1:661146/1‑149 (MQ=255)
                                              cAAGAAGTAGACCACTGACTGATATCCAGAGATTTCTACGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTAATTTTCCAGTGCCATATTGTTCCTcac                                                                                                  >  2:616478/1‑149 (MQ=255)
                                               aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGGTTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTACCATATTGTTCCTaaca                                                                                                 >  2:897871/1‑149 (MQ=255)
                                                 gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCGGGATTTTCTGGTTTCCCGTCGCGAAAAAAGATAACGACTATTTTTATCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                               >  2:143760/1‑149 (MQ=255)
                                                   agaaGACAACTAACTGATATCCAGAGATTCCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                             <  1:865086/146‑1 (MQ=255)
                                                   aataaaCAACTAACTGATATCCAAAAATTTCTTCAAGCAAAAACATACACGTTATTTTCCCAGAAGGCCGGATTTTCTGGTTGCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                             <  1:338250/144‑1 (MQ=255)
                                                     tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                           >  1:222555/1‑149 (MQ=255)
                                                        aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                        >  2:7844/1‑149 (MQ=255)
                                                            ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                    <  1:738779/149‑1 (MQ=255)
                                                            ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                    >  1:280109/1‑149 (MQ=255)
                                                                 tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                               >  1:292459/1‑149 (MQ=255)
                                                                   atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCGAGAGTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTTGTAATGaaa                                                                             >  1:464522/1‑149 (MQ=255)
                                                                   atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTGAGGaaa                                                                             >  2:397639/1‑149 (MQ=255)
                                                                       ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACCCATTGGTAATGaaaaaaa                                                                         >  2:757084/1‑149 (MQ=255)
                                                                       caagagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGTTTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                         <  1:994268/147‑1 (MQ=255)
                                                                        cAGAGATTTCTTCGAGCAATATGTTTTCCTTTCTTTTCACAGCAATCCCGATTTTCTGGTTTCCAGTAGCGAAAAAAGATAACGACTATTTTTCTCCTTTTGCAGTGCCATATTGTTCCTCACATTCACACATTGTTAATGAAAAAAAg                                                                        >  2:174279/1‑149 (MQ=255)
                                                                         agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGTCTATTTTTCTCATTTTTCGGTGCCATATTGTTTCTAACATGCACAAATTGCTATTGAAAATAAGa                                                                       >  1:10544/1‑149 (MQ=255)
                                                                          gagaTTTCTGCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                      >  2:10683/1‑149 (MQ=255)
                                                                                ttttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATGTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAaca                                                                <  1:788984/147‑1 (MQ=255)
                                                                                tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATAGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaaaa                                                                >  1:648104/1‑147 (MQ=255)
                                                                                 cttcGAGCAATTTGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTTGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                               <  1:1004694/149‑1 (MQ=255)
                                                                                  ttcGAGGAATAGGTTTTCCTTTATGTACCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                              <  2:806674/149‑1 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTTGCGAAAAAAGATAACGACTATTTTTTTTTGGTTGCAGTGTCATATTGTTTGTCACATGTACAGTTTGGTAATGAAAATATAACAAAACACg                                                              >  1:502199/1‑149 (MQ=255)
                                                                                  ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGCGGCGAAAAAAGATAACGACTCTTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAAAGACAAAACACg                                                              >  1:379336/1‑149 (MQ=255)
                                                                                   tggAGCAATATTGTTTCCTTTATTGTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAATAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCATCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                             <  2:927300/147‑1 (MQ=255)
                                                                                       gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCAAATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTa                                                         >  1:741260/1‑149 (MQ=255)
                                                                                         aaTATTTTTTCCTTTATTTGCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                       <  1:650959/149‑1 (MQ=255)
                                                                                            atGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                    <  1:10683/149‑1 (MQ=255)
                                                                                               ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                 >  2:62756/1‑149 (MQ=255)
                                                                                               ttttCCTTGATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                 <  2:940655/149‑1 (MQ=255)
                                                                                               tttaCCTTTATTTTCCCAGAATTCCGGATTTTTTGGTTTCCAATCGTTAAAAAAGATAATGACTTTTTTTTTCTTTTTGCAGTGCCACATTGTCCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                 <  1:295882/149‑1 (MQ=255)
                                                                                                  tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGATAATGAACAAAAGACAAAACAAGGGATAAGGCGCAagag                                              >  2:282707/1‑149 (MQ=255)
                                                                                                     tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTACTGAAAAAAAGACAAAACACGAGGCAAGGCGCAATAGCCa                                           >  2:543496/1‑149 (MQ=255)
                                                                                                           ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                     <  1:62756/149‑1 (MQ=255)
                                                                                                              ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                  >  1:627793/1‑149 (MQ=255)
                                                                                                              ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAAAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                  >  1:627852/1‑149 (MQ=255)
                                                                                                                   aGGCCGGATTTTCTGGTACCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                             <  2:685851/149‑1 (MQ=255)
                                                                                                                       cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAAGAACTATTTTTCTCATTTTACAGTGCCATATTGTTCATCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGGTGGAAGGCGCAATAGCAAGTTACTAAAATTAAGGaa                         >  2:687505/1‑148 (MQ=255)
                                                                                                                       cGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGCCAAAACACGAGGTAAGACGCAAAAGCCAGTTATTAGAATTAAGAAt                         >  2:451201/1‑149 (MQ=255)
                                                                                                                        ggATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg                        <  1:564796/149‑1 (MQ=255)
                                                                                                                        ggATTTTCTGGTTTCCAGTAGCTAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAACACACGCGGTAAGGCGCAAAAGCCAGTTATTAGAATTAAGGAtc                        >  2:739334/1‑148 (MQ=255)
                                                                                                                                 ggTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACACGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAGTAGCCAGTCATTAGGATTAAGGATGAATTGGGTg               >  1:944179/1‑149 (MQ=255)
                                                                                                                                    ttCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg            <  1:991975/149‑1 (MQ=255)
                                                                                                                                     tCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGt           >  2:842909/1‑149 (MQ=255)
                                                                                                                                             tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGGGCTGa      >  1:492264/8‑149 (MQ=255)
                                                                                                                                                 gaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa  <  2:292459/146‑1 (MQ=255)
                                                                                                                                                  |||                                                                                                                                                  
TTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGA  >  NC_000913_3_pae_tpiA/1543879‑1544168

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