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.A81.5% 158.2 / 29.8 54coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (4/6);  new base A (31/13);  total (35/19)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 185.7 / NA 54L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (35/19);  total (35/19)

TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543880‑1544170
                                                                                                                                               ||||                                                                                                                                                    
tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑‑a‑aaaag                                                                                                                                                >  1:607999/1‑147 (MQ=255)
    ggtAACGGTAATTGCGCTAAATAGCATCAAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCGTCGATCAATATGTTTTCCTTTAGTTTCCCAGAAGGCCGGATTTTCTGGTTTACAGACGAGAA‑‑A‑AAAGATaa                                                                                                                                             <  2:478759/147‑1 (MQ=255)
           gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGATTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTAt                                                                                                                                      >  1:516299/1‑149 (MQ=255)
                tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTtc                                                                                                                                 >  1:1137295/1‑149 (MQ=255)
                     tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATCTCCAGAAATTTCTTCGAGCAATATGTTTTCATTTATTTTCCCAGAAGGCCGGATTTTCTGATTTCCAGTCGCGAA‑‑A‑AAAGCTAACGACTATTTTTCTAAtt                                                                                                                            >  2:509300/1‑149 (MQ=255)
                          aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:967623/149‑1 (MQ=255)
                                         aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                        >  1:582119/1‑149 (MQ=255)
                                          aGGCAAGAAGTAGACAACTGACGGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑‑ACAAAGATAACGACTACTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                       <  1:1081040/149‑1 (MQ=255)
                                              aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCGGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTTTCATTTTACCGTGCCATATTGTTCTTcaca                                                                                                   >  1:245110/1‑149 (MQ=255)
                                              aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                   <  2:911920/149‑1 (MQ=255)
                                                      gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacaca                                                                                           >  2:1029173/1‑149 (MQ=255)
                                                       aCAACTGACTGATATTCGGAGATCTCTTCGATAAATATATTTCCCTTTATCGTATCAGTAAGCCGGATTTTC‑GGTTGAATATAAAA‑AAAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTTCTCACATGCACACAt                                                                                          <  2:937008/149‑1 (MQ=255)
                                                         aaCTGACTGATATCCAGAGATTTCTGCGAGCAATATGGTCTCCTTTATTTTCCCAGAAGGCCGGGTTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACCATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  1:1019828/149‑1 (MQ=255)
                                                          actgactgATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  1:1043607/149‑1 (MQ=255)
                                                            tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGTCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  2:84909/1‑149 (MQ=255)
                                                             gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    <  1:688460/149‑1 (MQ=255)
                                                              actgaTATCCAGAGATTTATTCGAGCAATATGTTTTCCTTTATGTTCTCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAATGGTAAt                                                                                   >  1:329979/1‑149 (MQ=255)
                                                                tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCAACTGCACACATTGTTAATGa                                                                                 >  1:466756/1‑149 (MQ=255)
                                                                  atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTCTGCAGTGCCATATTGTTCCTCACATGCACACCTTGGTAATGaaa                                                                               >  2:20511/1‑149 (MQ=255)
                                                                  atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGCTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTCTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATaaaa                                                                               >  2:595980/1‑149 (MQ=255)
                                                                     tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAATGCCATATTGTTCCTCACATGCACACATTGGTAGTGaaaaaa                                                                            >  2:38864/1‑149 (MQ=255)
                                                                       cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCAATATTGTTCCTCACATGCACACATTGTTAATGaaaaaaaa                                                                          >  2:735716/1‑148 (MQ=255)
                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  1:119501/1‑149 (MQ=255)
                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGAAATCGTCATTTTTTCTCATTTTGCAGTGCAATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         <  1:1071235/149‑1 (MQ=255)
                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  1:579364/1‑149 (MQ=255)
                                                                          agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                       >  2:430474/1‑149 (MQ=255)
                                                                             tttCTTCGAGCAATACGATTTCCTTTATTTTCCCAGAAGGCCGGATCTCCTGGTTTCCAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTCCCATATTCTTCCTCACAAGCCCACATTGTCAATCAAAAAACGCCaaaa                                                                    >  2:372001/1‑149 (MQ=255)
                                                                                cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCCCACATTGGTAATGAAAAAAAGACAAAacac                                                                 >  2:916515/1‑149 (MQ=255)
                                                                                cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                 >  1:629693/1‑149 (MQ=255)
                                                                                 ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCATTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACACAacccc                                                                >  2:115269/1‑146 (MQ=255)
                                                                                  tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTTTTTTTTTTTTTTTCCATATAGTTCTTCACATGTACACATTGGTTATGTAAAAATTACAAAaccatc                                                               >  2:529296/1‑145 (MQ=255)
                                                                                     aGCAATATGTTTTCCTTTACTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                            <  1:943552/149‑1 (MQ=255)
                                                                                          tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCACTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGCGGTAAGgc                                                       >  2:1049279/1‑149 (MQ=255)
                                                                                            gtttttCTTGTAGTTTCCCAGAAGGCGGGATGTTATGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACGATTTTTCTCATTTTGCAGTGCCATGTTGTTCCTCACATGCACACCTTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcgc                                                     <  2:543767/147‑1 (MQ=255)
                                                                                              ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAgagggaa                                                   >  2:518858/1‑143 (MQ=255)
                                                                                              tgtaccTTTATTTTCACAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAGTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                   <  1:754667/145‑1 (MQ=255)
                                                                                                     ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  2:152539/1‑149 (MQ=255)
                                                                                                             ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                    >  2:926834/1‑149 (MQ=255)
                                                                                                              cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                   >  1:95098/1‑149 (MQ=255)
                                                                                                                   ggCCGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  2:582119/149‑1 (MQ=255)
                                                                                                                      cGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAGGGCGCAATAGCCAGTTGTTAGAATTAAGGAt                           >  1:422110/1‑149 (MQ=255)
                                                                                                                      cGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAgg                                                        <  2:954708/120‑1 (MQ=255)
                                                                                                                      cGGATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAgg                                                        >  1:954708/1‑120 (MQ=255)
                                                                                                                       ggATTTTCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGACGCAATAGCCAGTTATTAGAATTAAGGATg                          >  1:190463/1‑149 (MQ=255)
                                                                                                                          ttttCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      >  2:352863/1‑70 (MQ=255)
                                                                                                                          ttttCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      <  1:352863/70‑1 (MQ=255)
                                                                                                                            ttCTGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                     >  1:630510/1‑149 (MQ=255)
                                                                                                                               tGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGGATTAGGGATGAATTGGGt                  >  1:657412/1‑149 (MQ=255)
                                                                                                                               tGGTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                  >  1:309856/1‑149 (MQ=255)
                                                                                                                               tGGTTTCCAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                  <  1:97065/149‑1 (MQ=255)
                                                                                                                                 gTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:91293/145‑1 (MQ=255)
                                                                                                                                 gTTTCCAGTCGCGAA‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    >  2:91293/1‑145 (MQ=255)
                                                                                                                                     ccAGTCGCG‑AA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGATGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTg            <  1:430474/149‑1 (MQ=255)
                                                                                                                                       aGTCGCG‑AA‑A‑AAAGACAATGACTTATTCTCTCATTTTGCAGGTCCATATTGTTCCTCACATGCACACATTGGTAATTAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAATGATGAATTTGGTGAAGTGCt          <  2:502832/149‑1 (MQ=255)
                                                                                                                                           tcgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        >  1:405874/8‑149 (MQ=255)
                                                                                                                                            cgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  1:152539/143‑1 (MQ=255)
                                                                                                                                             gcgaaa‑aaaGAGAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt      <  2:939082/144‑1 (MQ=255)
                                                                                                                                              cgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg     >  2:428328/5‑149 (MQ=255)
                                                                                                                                                 aa‑aaaTATAACGACTATTTTACTCATTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:416/148‑1 (MQ=255)
                                                                                                                                               ||||                                                                                                                                                    
TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543880‑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: ‑