Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 988,555 A→G L342P (CTG→CCG)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913988,5550AG99.6% 880.5 / ‑4.1 275coding (1025/1203 nt)pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (139/135);  minor base . (0/1);  total (139/136)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 290 total aligned reads displayed.

CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644
                                                                                          |                                                                                         
cTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCggg                                                                                         <  1:5324874/93‑1 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      >  2:6661732/1‑93 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      <  1:2112229/93‑1 (MQ=255)
   tCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGa                                                                                      <  1:5913682/93‑1 (MQ=255)
    cAGAAAGTTTCGCCTCCGGTTGACCGTTACTCTCTACCAACTTAATGACAATTTTGAGGGTTGTTGCCTGGGGGATATCGCATGTCGGGCGAg                                                                                     <  2:5833460/93‑1 (MQ=255)
    cAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAg                                                                                     <  2:4063008/93‑1 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTc                                                                                   <  2:2801597/93‑1 (MQ=255)
      gAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTc                                                                                   >  2:1792600/1‑93 (MQ=255)
        aaGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTccc                                                                                 >  2:7593225/1‑93 (MQ=255)
          gTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCaa                                                                               <  1:4383400/93‑1 (MQ=255)
           tttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAAt                                                                              <  2:4433671/93‑1 (MQ=255)
            ttCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATa                                                                             >  1:9748080/1‑93 (MQ=255)
                 cACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaa                                                                        >  1:4977844/1‑93 (MQ=255)
                  cccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCaaaa                                                                       <  2:4147498/92‑1 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      >  2:6337541/1‑93 (MQ=255)
                   ccGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAAc                                                                      >  1:6319701/1‑93 (MQ=255)
                     ggTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGACCCAATACCAAAACtt                                                                    >  1:3329279/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:8203906/1‑93 (MQ=255)
                      gTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTa                                                                   >  1:9843560/1‑93 (MQ=255)
                       tttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTaa                                                                  <  2:7408171/93‑1 (MQ=255)
                        ttACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAt                                                                 <  2:1208268/93‑1 (MQ=255)
                         tACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAAtt                                                                >  2:5259612/1‑93 (MQ=255)
                          aCCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTg                                                               <  1:9304721/93‑1 (MQ=255)
                          aCCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTg                                                               >  2:8649862/1‑93 (MQ=255)
                           ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGc                                                              <  1:4308394/93‑1 (MQ=255)
                           ccGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGc                                                              <  1:6513673/93‑1 (MQ=255)
                            cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCa                                                             <  1:9641787/93‑1 (MQ=255)
                            cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCa                                                             <  2:7656418/93‑1 (MQ=255)
                            cGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCa                                                             <  2:8463946/93‑1 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            >  2:3587761/1‑93 (MQ=255)
                             gTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            >  2:2673776/1‑93 (MQ=255)
                             gTTACACTATACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAc                                                            >  2:7329240/1‑93 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           <  1:5042515/93‑1 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           <  1:7621397/93‑1 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           <  1:7505350/93‑1 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           <  1:4612149/93‑1 (MQ=255)
                              ttACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcg                                                           >  1:1068733/1‑93 (MQ=255)
                               tACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgc                                                          <  2:3121774/93‑1 (MQ=255)
                                acacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCAcgcg                                                         >  1:2874300/1‑93 (MQ=255)
                                 cacTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCgg                                                        <  2:4545313/93‑1 (MQ=255)
                                  cctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCggg                                                       <  2:7182572/92‑1 (MQ=255)
                                   ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGa                                                      >  2:1351387/1‑93 (MQ=255)
                                   ctctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGa                                                      >  1:8333433/1‑93 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     <  2:7045371/93‑1 (MQ=255)
                                    tctACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGaa                                                     >  1:2595727/1‑93 (MQ=255)
                                       aCCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAagag                                                  >  1:9402856/1‑93 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                >  1:7326120/1‑93 (MQ=255)
                                         cAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGaa                                                <  2:2033291/93‑1 (MQ=255)
                                          aaCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAg                                               <  2:8877604/93‑1 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGACGt                                              >  2:8189272/1‑93 (MQ=255)
                                           aCTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGt                                              <  1:9137215/93‑1 (MQ=255)
                                            cTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTg                                             <  1:5818278/93‑1 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:2160827/1‑93 (MQ=255)
                                             ttAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTgg                                            >  1:8711772/1‑93 (MQ=255)
                                              tAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGc                                           <  1:6082124/93‑1 (MQ=255)
                                                aTGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGCAGTGGCgg                                         >  1:68797/1‑93 (MQ=255)
                                                  gACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGta                                       <  2:9903401/93‑1 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      <  1:3896003/93‑1 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      <  2:5516982/93‑1 (MQ=255)
                                                   aCAATATTCAGGGGTTTTACCTGGGCGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGtat                                      <  2:4486413/93‑1 (MQ=255)
                                                         ttCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCa                                >  2:7001183/1‑93 (MQ=255)
                                                          tCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCaa                               <  2:6565801/93‑1 (MQ=255)
                                                          tCAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCaa                               >  1:9214102/1‑93 (MQ=255)
                                                           cAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAc                              <  2:9662835/93‑1 (MQ=255)
                                                           cAGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAc                              <  1:4346181/93‑1 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             <  2:4687238/93‑1 (MQ=255)
                                                            aGGGGTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAAcc                             >  1:4020883/1‑93 (MQ=255)
                                                             ggggTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTGAATTGCACGCGGGAAGAGAAGTGGCGGTAGAGCTCAACCg                            >  2:8921666/1‑93 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         >  2:753467/1‑93 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         <  1:8179322/93‑1 (MQ=255)
                                                                gTTTTACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCtt                         >  2:4938472/1‑93 (MQ=255)
                                                                   ttACCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGc                      >  1:8110699/1‑93 (MQ=255)
                                                                     aCCTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGt                    >  2:1473059/1‑93 (MQ=255)
                                                                      ccTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTa                   <  2:8707332/93‑1 (MQ=255)
                                                                      ccTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTa                   <  1:1838438/93‑1 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  <  2:4858021/93‑1 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  >  1:5528451/1‑93 (MQ=255)
                                                                       cTGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaa                  >  2:3170981/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  1:1797066/1‑93 (MQ=255)
                                                                        tGGGGGATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTaaa                 >  1:10021509/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAAt                >  1:7641118/1‑93 (MQ=255)
                                                                         gggggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAAt                <  2:6760781/93‑1 (MQ=255)
                                                                           gggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATcc              <  2:2509090/93‑1 (MQ=255)
                                                                           gggATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATcc              >  1:1343027/1‑93 (MQ=255)
                                                                            ggATATCGCAGGGCGGGCGAGTCCCAATACCAAACCTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCa             <  1:2923191/93‑1 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  1:7384220/1‑93 (MQ=255)
                                                                             gATATCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAg            >  2:8324041/1‑93 (MQ=255)
                                                                              atatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGa           >  2:426718/1‑93 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGAt          <  1:9025643/93‑1 (MQ=255)
                                                                               tatCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGAGGCGGTATAGCTCAACCGGTTTGCGTAAATGCCGAt          >  2:6905361/1‑93 (MQ=255)
                                                                                 tCGCAGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTg        >  1:2332646/1‑93 (MQ=255)
                                                                                     aGGTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCag    <  1:5713751/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   <  2:1522514/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   <  2:9318525/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  2:2374597/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   <  2:7952303/93‑1 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   >  1:1811761/1‑93 (MQ=255)
                                                                                      ggTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCaga   <  2:1899803/93‑1 (MQ=255)
                                                                                       gTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCagag  >  1:365270/1‑93 (MQ=255)
                                                                                       gTCGGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCagag  >  2:3867108/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGTCAGAAAGTTTCGCCACCGGTTTACCGTTACACTCTACCAACTTAATGACAATATTCAGGGGTTTTACCTGGGGGATATCGCAGGTCAGGCGAGTCCCAATACCAAAACTTAATTGCACGCGGGAAGAGAAGTGGCGGTATAGCTCAACCGCTTTGCGTAAATCCAGATTGTCAGAG  >  NC_000913/988465‑988644

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