Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,073,307 Δ1 bp 100% coding (705/1149 nt) rutA ← pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,073,3050A.60.3% 74.0 / 229.1 121coding (707/1149 nt)rutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base A (25/23);  new base . (30/43);  total (55/66)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 127 total aligned reads displayed.

GTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCA  >  NC_000913/1073208‑1073404
                                                                                                 |                                                                                                   
   ggTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaagccg                                                                                                >  1:1464154/1‑96 (MQ=255)
   ggTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGACCTGTATCATCCGCGCAGCGGTCGGGGCGaagccg                                                                                                >  1:3173622/1‑96 (MQ=255)
    gTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCg                                                                                               >  1:2905900/1‑100 (MQ=255)
       tCATCGGCAATCACCATACACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGCGCGAAGCCGTgg                                                                                            <  1:810162/99‑1 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTggg                                                                                           <  1:2355617/100‑1 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTggg                                                                                           <  2:2432065/100‑1 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTggg                                                                                           >  2:1163792/1‑99 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGgt                                                                                          >  2:1615410/1‑100 (MQ=255)
         aTCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGgtg                                                                                         >  2:3112605/1‑100 (MQ=255)
          tCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGgtg                                                                                         <  1:858565/99‑1 (MQ=255)
           cGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGgtgt                                                                                        <  2:439118/100‑1 (MQ=255)
            ggCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTa                                                                                       >  1:303704/1‑100 (MQ=255)
              cAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTAtt                                                                                     <  2:3022106/100‑1 (MQ=255)
                aTCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTa                                                                                   >  1:1805996/1‑100 (MQ=255)
                 tCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACg                                                                                 >  1:1363243/1‑100 (MQ=255)
                  cACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACg                                                                                 <  1:3278783/100‑1 (MQ=255)
                    ccATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCt                                                                              >  2:2102179/1‑100 (MQ=255)
                    ccATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCt                                                                              >  2:1685215/1‑100 (MQ=255)
                    ccATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCt                                                                              >  1:319917/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGTAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCt                                                                              >  2:3189968/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCtt                                                                             <  1:154944/100‑1 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCACAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCtt                                                                             <  1:197156/100‑1 (MQ=255)
                      aTAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCtt                                                                             >  1:1710271/1‑100 (MQ=255)
                       tAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCttt                                                                            >  2:704642/1‑100 (MQ=255)
                       tAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTg                                                                           <  2:1397929/100‑1 (MQ=255)
                        aaaCAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGc                                                                          <  2:1306477/100‑1 (MQ=255)
                         aacaacACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGc                                                                          <  2:1009225/100‑1 (MQ=255)
                          acaacaCATAAGAGCCACCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGcc                                                                         >  1:188229/1‑100 (MQ=255)
                           caacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGa                                                                       >  1:1041403/1‑100 (MQ=255)
                           caacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGa                                                                       >  1:3179208/1‑100 (MQ=255)
                            aacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGaa                                                                      <  1:511350/100‑1 (MQ=255)
                            aacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGaa                                                                      <  2:419352/100‑1 (MQ=255)
                              cacaTAAGAGCCAACGTCGCGTCCGGTTTTCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAAc                                                                    <  1:2043300/100‑1 (MQ=255)
                              cacaTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAAc                                                                    <  1:764429/100‑1 (MQ=255)
                                caTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACa                                                                   >  2:1533114/1‑100 (MQ=255)
                                caTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAg                                                                  >  1:312743/1‑100 (MQ=255)
                                    agagCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAg                                                               <  2:2912261/100‑1 (MQ=255)
                                     tagCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGt                                                              <  2:2573854/99‑1 (MQ=255)
                                      agCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGta                                                             >  1:2573834/1‑99 (MQ=255)
                                       gCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTg                                                            >  2:642052/1‑100 (MQ=255)
                                       gCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTg                                                            >  1:781943/1‑100 (MQ=255)
                                       gCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTg                                                            <  1:1211735/100‑1 (MQ=255)
                                       gCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGa                                                           >  1:2044562/1‑100 (MQ=255)
                                        ccAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCCGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGa                                                           <  1:1410102/99‑1 (MQ=255)
                                           aCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATc                                                       <  1:243155/100‑1 (MQ=255)
                                           aCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGATCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATc                                                       <  2:1872944/100‑1 (MQ=255)
                                             gTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCgg                                                     <  1:1435196/100‑1 (MQ=255)
                                             gTCGCGTCCGGTGTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAGAGCCGTGGGTGGCTTTACGCCTTTGCCGAAACAGAAGTTGACATCg                                                      <  1:1029996/100‑1 (MQ=255)
                                              tggcgTTCGGTTTGTGCTGCGGCCTGTTTCCTCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  1:2175885/98‑1 (MQ=255)
                                              tCGCGTCCGGTTTGCTCTGCGGCCTGTTTTATCCGCGCAGCTGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  1:472338/100‑1 (MQ=255)
                                              tCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  1:2295790/100‑1 (MQ=255)
                                              tCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  1:1280187/100‑1 (MQ=255)
                                              tCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATGTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  1:1327003/100‑1 (MQ=255)
                                               cgcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                    <  2:480103/100‑1 (MQ=255)
                                               cgcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCa                                                   <  1:3174450/100‑1 (MQ=255)
                                                gcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCAt                                                  >  2:2253377/1‑100 (MQ=255)
                                                  gTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATAC                                                <  1:1635977‑M1/100‑1 (MQ=255)
                                                     cGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGG                                             >  2:554076‑M1/1‑100 (MQ=255)
                                                      ggTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGG                                            <  1:2662370‑M1/100‑1 (MQ=255)
                                                      ggTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGG                                            <  1:2324298‑M1/100‑1 (MQ=255)
                                                           gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGga                                        >  2:2697438/1‑100 (MQ=255)
                                                           gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag                                       <  2:1198209/100‑1 (MQ=255)
                                                             tctGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgaga                                      >  2:13956/1‑100 (MQ=255)
                                                                gCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACg                                   >  1:1436861/1‑100 (MQ=255)
                                                                gCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGc                                  >  1:1921347/1‑100 (MQ=255)
                                                                 cGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGcc                                 <  1:2730035/100‑1 (MQ=255)
                                                                 cGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGcc                                 <  1:2730085/100‑1 (MQ=255)
                                                                  ggCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCa                                >  2:842147/1‑100 (MQ=255)
                                                                  ggCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCa                                >  2:839351/1‑100 (MQ=255)
                                                                   gCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAATGCCAt                               <  2:2266366/100‑1 (MQ=255)
                                                                   gCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCa                                <  1:3243964/99‑1 (MQ=255)
                                                                   gCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCAt                               >  1:2421550/1‑100 (MQ=255)
                                                                       gattCATCCGCGCGGCGGTCGGGGCGAAGCCTTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCa                           >  2:2285321/3‑100 (MQ=255)
                                                                       gTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGcc                            <  2:3359778/100‑1 (MQ=255)
                                                                       gTTTCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCa                           >  2:1995748/1‑100 (MQ=255)
                                                                         ttCATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGc                         >  2:2272680/1‑100 (MQ=255)
                                                                           cATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                        <  2:2879003/100‑1 (MQ=255)
                                                                           cATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                        >  2:570982/1‑100 (MQ=255)
                                                                           cATCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                        >  1:1244212/1‑99 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGt                       >  1:1605028/1‑100 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGt                       >  1:2347302/1‑100 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTc                      <  2:2592567/100‑1 (MQ=255)
                                                                             tCCGCGCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCg                     >  1:610155/1‑100 (MQ=255)
                                                                             tCCGCGCAGCGGTCGGGGCGAAGCCGTAGGTGTATTTACGCCTTGGCCGAAACAGAAGCTGAAAGCGGCATAACGGGCGGAGAACGCCATGCCAGCGTc                      <  1:209597/99‑1 (MQ=255)
                                                                               cgcgCAGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCt                   >  1:2003362/1‑100 (MQ=255)
                                                                                    aGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgc               <  1:2759879‑M2/100‑3 (MQ=255)
                                                                                    aGCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCaaagcgacgctggcatggcgttctccgcc              <  2:2463612‑M1/100‑30 (MQ=255)
                                                                                     gCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgcc              >  2:96481‑M2/1‑97 (MQ=255)
                                                                                     gCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgcc              >  2:1659633‑M2/1‑97 (MQ=255)
                                                                                     gCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgccc             >  2:2181905‑M2/1‑96 (MQ=255)
                                                                                     gCGGTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCaaagcgacgctggcatggcgttctccgccc             <  1:437402‑M1/100‑31 (MQ=255)
                                                                                       ggTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgcccg            <  1:748708‑M2/100‑6 (MQ=255)
                                                                                        gTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGACGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTgcccgg           <  2:1141333‑M2/100‑7 (MQ=255)
                                                                                        gTCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgc          <  1:2967275/100‑1 (MQ=255)
                                                                                         tCGGGGCGAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgcg         <  2:461566/100‑1 (MQ=255)
                                                                                          cGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCaaagcgacgctggcatggcgttctccgcccggta         <  2:2435240‑M1/100‑35 (MQ=255)
                                                                                          cGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgcg         >  2:1686219/1‑100 (MQ=255)
                                                                                            gggCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCaaagcgacgctggcatggcgttctccgcccggtatg       <  2:2877446‑M1/100‑37 (MQ=255)
                                                                                            gggCGAAGCCGTGGGTGTATTTACGCCGTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCaaagcgacgctggcatggcgtgctccgcccggtatgc      <  1:675098‑M1/100‑38 (MQ=255)
                                                                                                 aaaGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCa  >  2:2766761/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAGATCA  >  NC_000913/1073208‑1073404

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