Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 350,318 A→G Q361Q (CAA→CAG)  prpC → putative citrate synthase; propionate metabolism?

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913350,3180AG99.6% 864.8 / ‑5.3 277Q361Q (CAA→CAG) prpCputative citrate synthase; propionate metabolism?
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (131/145);  total (132/145)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 296 total aligned reads displayed.

TACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAAC  >  NC_000913/350229‑350407
                                                                                         |                                                                                         
  caacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCa                                                                                      <  1:3869343/93‑1 (MQ=255)
  caacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCa                                                                                      <  1:3047721/93‑1 (MQ=255)
  caacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCa                                                                                      <  2:601560/93‑1 (MQ=255)
     caTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGa                                                                                   <  2:3117973/93‑1 (MQ=255)
       tgatgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGaca                                                                                 >  1:3561329/1‑93 (MQ=255)
       tgatgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGaca                                                                                 >  2:2750319/1‑93 (MQ=255)
          tgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaaca                                                                              >  2:5119158/1‑93 (MQ=255)
          tgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaaca                                                                              >  1:4539514/1‑93 (MQ=255)
            ggCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACaaa                                                                            <  1:2188906/93‑1 (MQ=255)
            ggCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACaaa                                                                            >  1:6002870/1‑93 (MQ=255)
             gCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACacaa                                                                           >  2:2712622/1‑93 (MQ=255)
             gCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACaaaa                                                                           <  2:7860903/93‑1 (MQ=255)
             gCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACaaaa                                                                           <  1:5279103/93‑1 (MQ=255)
                ttCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTa                                                                        >  2:6613484/1‑93 (MQ=255)
                ttCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTa                                                                        >  1:3353302/1‑93 (MQ=255)
                ttCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCCCACATTATCGAACAGCGCCAGAACAACAAAATTa                                                                        >  1:5085781/1‑93 (MQ=255)
                 tCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTAt                                                                       >  1:6988564/1‑93 (MQ=255)
                  cccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATc                                                                      <  1:6717961/93‑1 (MQ=255)
                  cccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATc                                                                      <  2:4171809/93‑1 (MQ=255)
                    cACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCg                                                                    >  1:6527276/1‑93 (MQ=255)
                       cGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTcc                                                                 >  2:1916887/1‑93 (MQ=255)
                        gagaTGTTCACCCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCt                                                                <  1:2651024/93‑1 (MQ=255)
                        gagaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTTTCCGTCCt                                                                >  2:6452368/1‑93 (MQ=255)
                        gagaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCt                                                                >  2:1538861/1‑93 (MQ=255)
                        gagaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCt                                                                <  1:208231/93‑1 (MQ=255)
                          gaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTc                                                              >  2:4771377/1‑93 (MQ=255)
                          gaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTc                                                              >  1:1393812/1‑93 (MQ=255)
                           aTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGTGCGGCGCACATTATCGAACAGCGTCAAGACAACAAAATTATCCGTCCTTcc                                                             >  1:3324135/1‑93 (MQ=255)
                            tGTTCACACCACAGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCACACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCg                                                            >  1:7701560/1‑93 (MQ=255)
                             gTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACCAAATTATCCGTCCTTCCGc                                                           >  1:5471736/1‑93 (MQ=255)
                             gTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGc                                                           >  2:4665577/1‑93 (MQ=255)
                             gTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGc                                                           <  2:6711246/93‑1 (MQ=255)
                              ttCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGcc                                                          >  1:4225502/1‑93 (MQ=255)
                                 acacCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAAt                                                       >  2:2921190/1‑93 (MQ=255)
                                   accacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTa                                                     <  2:3049468/93‑1 (MQ=255)
                                   accacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTa                                                     <  1:2845038/93‑1 (MQ=255)
                                     cacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATg                                                   <  1:5127626/93‑1 (MQ=255)
                                     cacTGTTTGTTATCGCCCGCGTCAACGGCTGGGCGGCGCACATTATAGAACCGCGACAGGACAACAAAATTAAACGTCCTTCACCCAACTATg                                                   >  2:4078828/1‑93 (MQ=255)
                                      acTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGt                                                  >  1:2334935/1‑93 (MQ=255)
                                       cTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGtt                                                 <  2:355709/93‑1 (MQ=255)
                                       cTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGtt                                                 <  1:3990101/93‑1 (MQ=255)
                                        tgtttgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAGATTATCCGTCCTTCCGCCAATTATGTTg                                                <  2:7204871/93‑1 (MQ=255)
                                        tgtttgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTg                                                >  2:447040/1‑93 (MQ=255)
                                          tttgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGa                                              >  2:3845057/1‑93 (MQ=255)
                                          tttgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGa                                              >  2:792756/1‑93 (MQ=255)
                                            tgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGAcc                                            <  1:1048328/93‑1 (MQ=255)
                                             gttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCg                                           <  2:4249372/93‑1 (MQ=255)
                                             gttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCg                                           <  1:6384169/93‑1 (MQ=255)
                                              ttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCgg                                          >  2:7415758/1‑93 (MQ=255)
                                               tATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGa                                         <  2:4975668/93‑1 (MQ=255)
                                                ttCGCCCGGGGCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGaa                                        <  1:7818300/92‑1 (MQ=255)
                                                aTCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGaa                                        <  1:53754/93‑1 (MQ=255)
                                                 tCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAg                                       >  2:6352440/1‑93 (MQ=255)
                                                    cccGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                    >  2:6581650/1‑93 (MQ=255)
                                                    cccGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                    >  1:5623127/1‑93 (MQ=255)
                                                    cccGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                    >  1:5526408/1‑93 (MQ=255)
                                                    cccGCGTCACCGGCTGGGCAGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                    >  1:3153265/1‑93 (MQ=255)
                                                       gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                 <  1:3253956/93‑1 (MQ=255)
                                                       gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                 <  1:4984696/93‑1 (MQ=255)
                                                       gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                 >  2:1379588/1‑93 (MQ=255)
                                                        cgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGccc                                <  2:1402609/93‑1 (MQ=255)
                                                        cgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGccc                                <  1:4673228/93‑1 (MQ=255)
                                                         gTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCg                               >  2:1966313/1‑93 (MQ=255)
                                                          tCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGt                              >  1:867067/1‑93 (MQ=255)
                                                          tCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGt                              >  2:5350233/1‑93 (MQ=255)
                                                           cACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGtt                             <  1:4727512/93‑1 (MQ=255)
                                                           cACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGtt                             <  1:7427712/93‑1 (MQ=255)
                                                            cccGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                            <  2:5813718/92‑1 (MQ=255)
                                                            aCCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                            <  2:2974106/93‑1 (MQ=255)
                                                            aCCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                            >  1:6728016/1‑93 (MQ=255)
                                                             ccGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTg                           <  2:5884338/93‑1 (MQ=255)
                                                             ccGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTg                           >  2:3931140/1‑93 (MQ=255)
                                                              cGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGt                          >  2:1758987/1‑93 (MQ=255)
                                                              cGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGt                          >  1:1687229/1‑93 (MQ=255)
                                                               ggCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTc                         >  2:1267170/1‑93 (MQ=255)
                                                                gCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcg                        >  1:5830118/1‑93 (MQ=255)
                                                                gCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcg                        >  1:973155/1‑93 (MQ=255)
                                                                  tGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgcg                      >  1:1900928/1‑93 (MQ=255)
                                                                    ggcggcGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCt                    >  2:3212742/1‑93 (MQ=255)
                                                                        gcgcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGAt                <  1:2749299/93‑1 (MQ=255)
                                                                        gcgcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGAt                <  2:2497292/93‑1 (MQ=255)
                                                                         cgcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATa               <  2:381514/93‑1 (MQ=255)
                                                                          gcACATTATCGACCAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCAGGCTaa              >  2:2264870/1‑93 (MQ=255)
                                                                           cacaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAg             <  2:54644/93‑1 (MQ=255)
                                                                           cacaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAg             <  1:4354353/93‑1 (MQ=255)
                                                                            acaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTCTGCTGGACTGGAAGACCGCCCGTTTGTCGCGATGGATAAgc            >  1:4099528/1‑93 (MQ=255)
                                                                            acaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgc            <  2:6915856/93‑1 (MQ=255)
                                                                              aTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgcgc          <  2:6772607/93‑1 (MQ=255)
                                                                               ttATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGcc         >  1:1076271/1‑93 (MQ=255)
                                                                               ttATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCCCGCTGGATAAGCGcc         >  2:1844432/1‑93 (MQ=255)
                                                                                tATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCa        <  1:3733291/93‑1 (MQ=255)
                                                                                tATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCa        <  2:773587/93‑1 (MQ=255)
                                                                                 aTCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAg       <  2:3330969/93‑1 (MQ=255)
                                                                                 aTCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAg       >  2:7506978/1‑93 (MQ=255)
                                                                                     aacagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaaa   <  1:4122806/93‑1 (MQ=255)
                                                                                     aacagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaaa   >  1:7329521/1‑93 (MQ=255)
                                                                                     aacagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaaa   >  2:7626113/1‑93 (MQ=255)
                                                                                     aacagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaaa   <  2:4220541/93‑1 (MQ=255)
                                                                                      acagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAAc  <  1:6161146/93‑1 (MQ=255)
                                                                                      acagcGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAAc  <  1:7263617/93‑1 (MQ=255)
                                                                                         |                                                                                         
TACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAAC  >  NC_000913/350229‑350407

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