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.A83.8% 230.5 / 41.7 74coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (5/7);  new base A (35/27);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
*NC_000913_3_pae_tpiA1,544,0220CA98.7% 264.7 / NA 74L13F (TTG→TTTnarUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (39/34);  minor base T (1/0);  total (40/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 = 1.00e+00

TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543877‑1544172
                                                                                                                                                   |||                                                                                                                                                     
tagtGAGATTAACGGTAGCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaa                                                                                                                                                    <  4:182042/147‑1 (MQ=255)
   tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAaaaaaga                                                                                                                                                 <  1:56442/150‑4 (MQ=255)
        ttAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACg                                                                                                                                             >  2:727349/1‑150 (MQ=255)
              gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAtt                                                                                                                                       >  3:568033/1‑150 (MQ=255)
                aaCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAtttt                                                                                                                                     >  1:75240/1‑150 (MQ=255)
                 aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTAttttt                                                                                                                                    >  3:448311/1‑150 (MQ=255)
                   tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATACAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGATAAACGATAAAGACTATTTTtct                                                                                                                                  >  3:241875/1‑150 (MQ=255)
                      gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCAt                                                                                                                               <  3:90751/150‑1 (MQ=255)
                        tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCAttt                                                                                                                             >  4:49294/1‑150 (MQ=255)
                                        acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCGTATTGt                                                                                                             <  1:105912/150‑1 (MQ=255)
                                        acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                             <  1:344163/150‑1 (MQ=255)
                                          aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           >  2:468672/1‑150 (MQ=255)
                                          aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           <  4:430823/150‑1 (MQ=255)
                                            aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTGCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                         <  2:192758/150‑1 (MQ=255)
                                              ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                       <  2:748418/150‑1 (MQ=255)
                                                cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                     >  1:575571/1‑150 (MQ=255)
                                                 aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                    >  4:441998/1‑150 (MQ=255)
                                                  agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                   >  1:58157/1‑150 (MQ=255)
                                                  agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                   <  1:411757/150‑1 (MQ=255)
                                                   gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                                  >  2:229205/1‑150 (MQ=255)
                                                          aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                                           >  4:287288/1‑150 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                                >  4:352661/1‑143 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                                <  3:352661/143‑1 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                         >  1:193894/1‑150 (MQ=255)
                                                            aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                         <  2:575571/150‑1 (MQ=255)
                                                              ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                       >  3:503807/1‑150 (MQ=255)
                                                               tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      <  4:580689/150‑1 (MQ=255)
                                                               tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      <  4:38211/150‑1 (MQ=255)
                                                                 actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATg                                                                                    >  4:119120/1‑150 (MQ=255)
                                                                       atCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                              >  3:497118/1‑150 (MQ=255)
                                                                        tCCAGAGATTTCTTCCAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCAAGTCGCGAAAAAAGACAACGACTATTTTTCTCATTTTGCAGTGACATATTGTTCCTCACCTGCACAAATAAATAAAGaaaaaaa                                                                             <  2:819374/150‑1 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                            <  4:568033/150‑1 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                            >  3:189885/1‑150 (MQ=255)
                                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCatat                                                                                                                <  3:462947/112‑1 (MQ=255)
                                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCatat                                                                                                                >  4:462947/1‑112 (MQ=255)
                                                                           agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                          >  1:115001/1‑150 (MQ=255)
                                                                            gagaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                         <  2:193894/150‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACCTATTTGTAATGAAAAAAAGACAAAACACGa                                                                 <  4:448311/150‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 >  1:448647/1‑150 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 <  1:727349/150‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 >  4:216742/1‑150 (MQ=255)
                                                                                         gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                            >  2:520504/1‑150 (MQ=255)
                                                                                             tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                        <  3:49294/150‑1 (MQ=255)
                                                                                                 ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                   >  3:252927/1‑150 (MQ=255)
                                                                                                   ttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAg                                                  >  1:820623/1‑150 (MQ=255)
                                                                                                    tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                                 >  2:872166/1‑150 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATGGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                               >  3:347393/1‑150 (MQ=255)
                                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                              <  3:287288/150‑1 (MQ=255)
                                                                                                        ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                             <  1:872166/150‑1 (MQ=255)
                                                                                                         tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                            >  3:487955/1‑150 (MQ=255)
                                                                                                           ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                          >  3:54807/1‑150 (MQ=255)
                                                                                                             accccAGAAGGCCCGATTTTCTGGTTCACAACCGCCAAAAAAGATAACGACTATTTTTCTCATTTTGAAGGGCCATTTTATTCATCACAAGCACACATTGTTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                        <  4:386695/148‑1 (MQ=255)
                                                                                                               cccAGATGGCCGGGTTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                      <  1:874961/150‑1 (MQ=255)
                                                                                                               cccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                      >  3:503137/1‑150 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                     <  3:130552/150‑1 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                     <  2:708481/150‑1 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                     >  4:330301/1‑150 (MQ=255)
                                                                                                                 cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTTTTTTTTTTTTGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                    >  2:916789/1‑150 (MQ=255)
                                                                                                                 cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                    >  1:570497/1‑150 (MQ=255)
                                                                                                                     gggCGGGATTTTCTGGTTTCCGGTCGCGAAAAAAGATAGCGAGTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                                <  2:638966/149‑1 (MQ=255)
                                                                                                                     aGGCCGGATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                                <  2:527448/150‑1 (MQ=255)
                                                                                                                          ggATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa                           >  1:869230/1‑150 (MQ=255)
                                                                                                                           gATTTTCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGaa                          >  3:249135/1‑150 (MQ=255)
                                                                                                                               ttCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                      <  1:504063/150‑1 (MQ=255)
                                                                                                                               ttCTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                      <  1:229205/150‑1 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                    <  4:539070/150‑1 (MQ=255)
                                                                                                                                  tGGTTTCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTg                   >  4:78332/1‑150 (MQ=255)
                                                                                                                                      ttCCAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGt               >  4:47692/1‑150 (MQ=255)
                                                                                                                                        ccAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc             <  3:330301/150‑1 (MQ=255)
                                                                                                                                         cAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                 >  1:635675/1‑145 (MQ=255)
                                                                                                                                         cAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                 <  2:635675/145‑1 (MQ=255)
                                                                                                                                         cAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGcg            >  1:917295/1‑149 (MQ=255)
                                                                                                                                         cAGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt            <  4:487955/150‑1 (MQ=255)
                                                                                                                                          aGTCGCGAAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTg           >  1:937206/1‑150 (MQ=255)
                                                                                                                                              tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt         >  1:814796/8‑150 (MQ=255)
                                                                                                                                               cgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt        <  4:503807/144‑1 (MQ=255)
                                                                                                                                                gcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg       <  2:224566/145‑1 (MQ=255)
                                                                                                                                                   aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  4:497118/148‑1 (MQ=255)
                                                                                                                                                   aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  3:119120/148‑1 (MQ=255)
                                                                                                                                                    aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg   >  3:488925/2‑150 (MQ=255)
                                                                                                                                                     aaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGa  >  3:76285/1‑150 (MQ=255)
                                                                                                                                                   |||                                                                                                                                                     
TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGCAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543877‑1544172

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