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.A87.0% 208.6 / 25.4 69coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (8/1);  new base A (21/39);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 3.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 229.8 / NA 69L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/40);  total (29/40)

AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543882‑1544170
                                                                                                                                              |||                                                                                                                                                   
agcttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑aaaaagat                                                                                                                                              >  2:1113511/4‑145 (MQ=255)
  aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAATTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATaa                                                                                                                                             <  2:562500/149‑1 (MQ=255)
  aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATaa                                                                                                                                             <  1:1356169/149‑1 (MQ=255)
   ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAAc                                                                                                                                            <  2:373927/149‑1 (MQ=255)
   ttAACGGGAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAAc                                                                                                                                            <  2:374518/149‑1 (MQ=255)
        ggTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTa                                                                                                                                       >  2:1283060/1‑149 (MQ=255)
           aaCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAttt                                                                                                                                    <  1:981405/149‑1 (MQ=255)
             cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGACAGCATCTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAttttt                                                                                                                                  <  1:1421691/149‑1 (MQ=255)
                        aGCATCCTTACACAGAACGCACGACGTATACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:328296/149‑1 (MQ=255)
                          cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                     <  1:1163980/149‑1 (MQ=255)
                             ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  1:219582/149‑1 (MQ=255)
                               aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                >  2:230353/1‑149 (MQ=255)
                               aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                >  1:976194/1‑149 (MQ=255)
                                        aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                       <  1:1283060/149‑1 (MQ=255)
                                          gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  1:819234/149‑1 (MQ=255)
                                          gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  1:332213/149‑1 (MQ=255)
                                           cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCAAGAAGGCCGGATTTTCTGTTTTCCAGTCGC‑GAA‑AAAAGATAACGACTGTTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                    >  2:1420762/1‑149 (MQ=255)
                                             agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCTAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  <  1:853445/149‑1 (MQ=255)
                                                agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTTCCATATTGTTCCTCACATGca                                                                                               <  2:1273829/149‑1 (MQ=255)
                                                agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                               <  1:914230/149‑1 (MQ=255)
                                                  tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGGGCCCTCTTGTTCCTCACATGccca                                                                                             >  2:422622/1‑149 (MQ=255)
                                                    gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                           >  1:1038661/1‑149 (MQ=255)
                                                       aaCTTACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  1:1400308/149‑1 (MQ=255)
                                                       aaCTGACCGATACCCAGAGCTCTTTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  1:422622/149‑1 (MQ=255)
                                                        actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       >  1:947814/1‑149 (MQ=255)
                                                          tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  2:215720/1‑149 (MQ=255)
                                                          tgaatcatatCTACAGTTTTCTTCGATCAATATTTTTGCCTTTATTTTCCCATAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTTCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  2:263643/143‑1 (MQ=255)
                                                              tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  <  1:873244/132‑1 (MQ=255)
                                                              tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  >  2:873244/1‑132 (MQ=255)
                                                              tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                                 >  1:615111/1‑149 (MQ=255)
                                                               gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACCCATTGGTAATGaa                                                                                >  2:693791/1‑149 (MQ=255)
                                                                    cgagagATTTCTGTGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:609085/147‑1 (MQ=255)
                                                                    ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           >  2:1215309/1‑149 (MQ=255)
                                                                    ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGAAGACTTTTCTTGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:768271/149‑1 (MQ=255)
                                                                    ccAGAGATCTCGTCCAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:261082/149‑1 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         <  1:737053/149‑1 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  2:138356/1‑149 (MQ=255)
                                                                      agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGATTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACAAATTGGTAGTGAAAAAAAGa                                                                         >  2:1272124/1‑149 (MQ=255)
                                                                        agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAATGACa                                                                       >  1:834341/1‑149 (MQ=255)
                                                                           tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                       <  1:831616/114‑1 (MQ=255)
                                                                           tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                       >  2:831616/1‑114 (MQ=255)
                                                                                tctagaaaTATGTTTTCCTTTATTTACACAGATGGCCGGATGTTCTGGTTTCCAGTCGCGG‑A‑AAAAGATAACGACTACTTGTCTCATTTTGCAGTGCAATTATGTTCTTCACACACACACCTTGGTAATGAAAAAAAGACAAAACACGa                                                               <  2:14151/143‑1 (MQ=255)
                                                                                tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  1:52958/149‑1 (MQ=255)
                                                                                  gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAACCACGAgg                                                             >  1:1019290/1‑149 (MQ=255)
                                                                                  gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAACGAACAAAAACGAgg                                                             >  2:95320/1‑149 (MQ=255)
                                                                                         atGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                      <  1:107476/149‑1 (MQ=255)
                                                                                            ttttCCTTTAATTTCCCCGATGTCCTGTTTTTCTGGTTTCCAGTCGA‑GAA‑AAAAGATAACGACTATTTTTCTAATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                   <  1:553022/149‑1 (MQ=255)
                                                                                               tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAg                                                >  2:903955/1‑149 (MQ=255)
                                                                                                ccTTTATTTTGCCAGAAGGGCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGATTTTTTTTCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                               <  2:1320378/149‑1 (MQ=255)
                                                                                                        ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTGTTCTCATTTTGCAGTGCCATATTGTTCCTAACATGCACACATAGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       <  1:343339/149‑1 (MQ=255)
                                                                                                            gaGAAGGCCGGTTTTTCGGTTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                   <  2:834341/148‑1 (MQ=255)
                                                                                                            cAGAAGGCCGGAGTTTCCGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                   <  2:534015/149‑1 (MQ=255)
                                                                                                                aGGCCGGATTTACTGGAATCAAGACGC‑GAA‑AAAAGATAACGACTTTTTTTCACATTTTGCAGTGCTATATTGTTCCTCACATGCACACAATAAAAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:947814/149‑1 (MQ=255)
                                                                                                                 ggCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  2:280823/65‑1 (MQ=255)
                                                                                                                 ggCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  >  1:280823/1‑65 (MQ=255)
                                                                                                                  gCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTAATTAGAATTAAgg                             >  1:1355939/1‑149 (MQ=255)
                                                                                                                    cGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                           >  2:1102258/1‑149 (MQ=255)
                                                                                                                    cGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                           >  1:1370207/1‑149 (MQ=255)
                                                                                                                    cGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                           >  2:510380/1‑149 (MQ=255)
                                                                                                                     ggATTTTCTGGTTTGCAGGCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTTCTCACATGCATACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGTCAGTTATTAGAATTAAGGATg                          <  2:513837/149‑1 (MQ=255)
                                                                                                                     ggATTTTCTGGTTTCCAGTCGCAA‑A‑AAAAGAAAACGACTATTTTTCTCATTTTGCAGTGCCAAATTGATACTCACATGCACACATTGGTAATCAAAAAAAGACAAAACACCAGGTAAGGCGAAATAGCTCATAATCAGAATTAAgactg                          >  2:1415314/1‑145 (MQ=255)
                                                                                                                      gATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa                         <  1:307651/149‑1 (MQ=255)
                                                                                                                           tCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:230353/149‑1 (MQ=255)
                                                                                                                            cTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  1:510380/149‑1 (MQ=255)
                                                                                                                            cTGGATTCCCGTCGC‑GAA‑AAAAGTTAACGCCTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAATGGTAATGAAAAAAAGACAAAACTCGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  1:786051/149‑1 (MQ=255)
                                                                                                                                tttCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa               >  2:148385/1‑149 (MQ=255)
                                                                                                                                gttCCAGTCGC‑GAA‑AAAAGATCACTACTATTTTTATCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa               <  2:208703/148‑1 (MQ=255)
                                                                                                                                    cAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATGGCCAGTTATTAGAATTAACAATGCATTGGGTGAAGGGc           >  1:552703/1‑149 (MQ=255)
                                                                                                                                     aGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTACAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAAACAAAACACGACGTAAGGCGCACTAGCCAGTTATTAGAACTAAGGATGAATTGGGTGAACTGCt          >  2:851126/1‑149 (MQ=255)
                                                                                                                                          cgcgaaaaaaGATAACGACTATGTTCCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGCATTAAGGATGAATTGGGTGAAGTGTTGAt       <  2:1103150/143‑1 (MQ=255)
                                                                                                                                           gcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt      <  1:359624/144‑1 (MQ=255)
                                                                                                                                               aaaaaGATAACGACTATTGTTCTAATATTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:95320/148‑1 (MQ=255)
                                                                                                                                              |||                                                                                                                                                   
AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543882‑1544170

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