Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,620,042 T→C W276R (TGG→CGG)  purM → phosphoribosylaminoimidazole synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,620,0420TC86.4% 718.8 / 97.4 295W276R (TGG→CGG) purMphosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base T (23/17);  new base C (144/111);  total (167/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 295 total aligned reads displayed.

GCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGAT  >  NC_000913/2619952‑2620130
                                                                                          |                                                                                        
gcgcACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTgg                                                                                        >  2:4304126/1‑93 (MQ=255)
     ccTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGa                                                                                   >  1:5466985/1‑93 (MQ=255)
     ccTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGa                                                                                   >  2:62690/1‑93 (MQ=255)
     ccTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGa                                                                                   >  1:4750652/1‑93 (MQ=255)
      cTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGaa                                                                                  >  1:692071/1‑93 (MQ=255)
      cTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGaa                                                                                  >  2:5721600/1‑93 (MQ=255)
       tGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAg                                                                                 <  2:2070626/93‑1 (MQ=255)
       tGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAg                                                                                 >  2:1098785/1‑93 (MQ=255)
       tGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAg                                                                                 >  1:6261358/1‑93 (MQ=255)
            ggcggcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTc                                                                            >  2:3252590/1‑93 (MQ=255)
              cggcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCaa                                                                          <  2:1905461/93‑1 (MQ=255)
                gcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACt                                                                        >  2:4078820/1‑93 (MQ=255)
                gcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACt                                                                        >  2:6586512/1‑93 (MQ=255)
                gcggcTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACt                                                                        >  1:1286998/1‑93 (MQ=255)
                    cTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCt                                                                    >  1:1340052/1‑93 (MQ=255)
                     ttCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTg                                                                   >  2:3124920/1‑93 (MQ=255)
                     ttCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTg                                                                   >  2:711592/1‑93 (MQ=255)
                      tCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGc                                                                  >  2:3014135/1‑93 (MQ=255)
                        tGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCaa                                                                <  1:6414887/93‑1 (MQ=255)
                        tGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCaa                                                                >  2:1277012/1‑93 (MQ=255)
                         gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCaaa                                                               >  1:5398961/1‑93 (MQ=255)
                         gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCaaa                                                               <  1:4105587/93‑1 (MQ=255)
                         gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCaaa                                                               <  1:6663903/93‑1 (MQ=255)
                         gggAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCaaa                                                               <  2:4789201/93‑1 (MQ=255)
                           gAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACg                                                             >  1:6558819/1‑93 (MQ=255)
                            aaaaCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACgg                                                            >  2:3181326/1‑93 (MQ=255)
                            aaaaCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACgg                                                            >  2:7524945/1‑93 (MQ=255)
                             aaaCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGc                                                           <  1:696797/93‑1 (MQ=255)
                             aaaCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGc                                                           >  2:5900980/1‑93 (MQ=255)
                              aaCATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCa                                                          <  2:2283532/93‑1 (MQ=255)
                                 aTTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGt                                                       >  1:4832550/1‑93 (MQ=255)
                                 aTTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGt                                                       >  1:2520373/1‑93 (MQ=255)
                                  ttCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTa                                                      >  2:5282148/1‑93 (MQ=255)
                                  ttCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTa                                                      >  1:911443/1‑93 (MQ=255)
                                  ttCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTa                                                      >  1:1883460/1‑93 (MQ=255)
                                  ttCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTa                                                      <  1:5912486/93‑1 (MQ=255)
                                    ccGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAAc                                                    <  1:6143391/93‑1 (MQ=255)
                                     cgcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACg                                                   <  2:852521/93‑1 (MQ=255)
                                     cgcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACg                                                   <  2:2401265/93‑1 (MQ=255)
                                     cgcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACg                                                   >  2:7309049/1‑93 (MQ=255)
                                      gcgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGt                                                  <  2:389216/93‑1 (MQ=255)
                                       cgcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGtt                                                 <  1:4459668/93‑1 (MQ=255)
                                        gcgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTg                                                <  1:1482344/93‑1 (MQ=255)
                                         cgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                               <  1:142206/93‑1 (MQ=255)
                                         cgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                               <  1:4234266/93‑1 (MQ=255)
                                         cgTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGa                                               <  1:3828533/93‑1 (MQ=255)
                                          gTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAg                                              >  2:70525/1‑93 (MQ=255)
                                           tATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGc                                             >  1:7440987/1‑93 (MQ=255)
                                           tATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGc                                             <  1:3774342/93‑1 (MQ=255)
                                            aTTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCa                                            >  2:1515153/1‑93 (MQ=255)
                                             ttGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCAc                                           >  1:502644/1‑93 (MQ=255)
                                             ttGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCAc                                           >  1:5508697/1‑93 (MQ=255)
                                              tGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCAcc                                          >  1:2458348/1‑93 (MQ=255)
                                              tGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCAcc                                          <  1:2357979/93‑1 (MQ=255)
                                                ccAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCGGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCAt                                        <  1:6030109/93‑1 (MQ=255)
                                                ccAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCAt                                        >  2:265602/1‑93 (MQ=255)
                                                 cAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATg                                       >  2:1019416/1‑93 (MQ=255)
                                                 cAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATg                                       >  2:4134511/1‑93 (MQ=255)
                                                 cAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATg                                       <  2:6120256/93‑1 (MQ=255)
                                                  aGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGa                                      >  2:5120756/1‑93 (MQ=255)
                                                  aGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGa                                      >  1:5241764/1‑93 (MQ=255)
                                                  aGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGa                                      >  2:2838844/1‑93 (MQ=255)
                                                  aGATAACACTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGa                                      <  2:4534332/93‑1 (MQ=255)
                                                   gATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGaa                                     >  1:7185418/1‑93 (MQ=255)
                                                    ataataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGaaa                                    <  1:3644518/93‑1 (MQ=255)
                                                     taataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                   >  2:6217922/1‑93 (MQ=255)
                                                     taataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                   >  1:2489838/1‑93 (MQ=255)
                                                     taataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAAt                                   <  2:1385698/93‑1 (MQ=255)
                                                     taataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGGAACGTTGAGCACCATGAAAt                                   >  1:6543682/1‑93 (MQ=255)
                                                      aataCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATg                                  >  2:2106441/1‑93 (MQ=255)
                                                      aataCTAAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATg                                  >  2:4722364/1‑93 (MQ=255)
                                                        taCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTa                                >  2:2885980/1‑93 (MQ=255)
                                                         aCTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTAt                               <  1:6356427/93‑1 (MQ=255)
                                                          cTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATc                              >  2:6495629/1‑93 (MQ=255)
                                                          cTCAGGCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATc                              <  1:2032364/93‑1 (MQ=255)
                                                              ggCAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGGGTTCAACTGGCTGCAAACGGCAGGGAACGTTGAGCACCATGAAATGTATCGCAc                          >  1:2620474/1‑93 (MQ=255)
                                                                cAGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCt                        <  2:5418854/93‑1 (MQ=255)
                                                                 aGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                       <  2:5556827/93‑1 (MQ=255)
                                                                 aGTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCtt                       >  2:6889308/1‑93 (MQ=255)
                                                                  gTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTc                      <  2:3639828/93‑1 (MQ=255)
                                                                  gTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTc                      >  2:4881079/1‑93 (MQ=255)
                                                                  gTGATTGATGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGATGAGCACCATGAAATGTATCGCACCTTc                      >  2:1671215/1‑93 (MQ=255)
                                                                   tgattgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCa                     <  2:1832353/93‑1 (MQ=255)
                                                                     attgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAAc                   <  1:2059012/93‑1 (MQ=255)
                                                                      ttgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACt                  >  1:2908521/1‑93 (MQ=255)
                                                                      ttgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACt                  >  1:3211829/1‑93 (MQ=255)
                                                                      ttgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACt                  <  1:3045712/93‑1 (MQ=255)
                                                                      ttgatGAATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACt                  >  1:275974/1‑93 (MQ=255)
                                                                       tgatgaATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTg                 >  2:2324626/1‑93 (MQ=255)
                                                                        gatgaATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgc                >  2:924168/1‑93 (MQ=255)
                                                                        gatgaATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgc                <  1:577385/93‑1 (MQ=255)
                                                                           gaATCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTgcggc             >  1:3391064/1‑93 (MQ=255)
                                                                             aTCTTCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGt           <  1:3537631/93‑1 (MQ=255)
                                                                                ttCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCgg        >  2:4587904/1‑93 (MQ=255)
                                                                                ttCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCgg        >  2:3731831/1‑93 (MQ=255)
                                                                                ttCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCgg        >  1:3956673/1‑93 (MQ=255)
                                                                                ttCCTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCgg        >  1:1083758/1‑93 (MQ=255)
                                                                                   cTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGgat     >  2:4469631/1‑93 (MQ=255)
                                                                                   cTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGgat     <  1:1792971/93‑1 (MQ=255)
                                                                                   cTGGCAGCGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGgat     >  1:1009304/1‑93 (MQ=255)
                                                                                          |                                                                                        
GCGCACCTGACCGGCGGCGGCTTCTGGGAAAACATTCCGCGCGTATTGCCAGATAATACTCAGGCAGTGATTGATGAATCTTCCTGGCAGTGGCCGGAAGTGTTCAACTGGCTGCAAACGGCAGGTAACGTTGAGCACCATGAAATGTATCGCACCTTCAACTGCGGCGTCGGGATGAT  >  NC_000913/2619952‑2620130

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