Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 503,440 T→C S161G (AGT→GGT)  fsr ← fosmidomycin resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913503,4400TC87.6% 634.5 / 69.9 234S161G (AGT→GGT) fsrfosmidomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (102/103);  total (118/116)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

Only 100 of 247 total aligned reads displayed.

ATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGA  >  NC_000913/503349‑503531
                                                                                           |                                                                                           
 tCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGc                                                                                           >  2:1380499/1‑93 (MQ=255)
  cGCCAGCAGCGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGcc                                                                                          <  2:5609995/93‑1 (MQ=255)
      agcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAg                                                                                      <  1:4713632/93‑1 (MQ=255)
      agcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAg                                                                                      <  1:2897715/93‑1 (MQ=255)
      agcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAg                                                                                      <  2:3033244/93‑1 (MQ=255)
       gcagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGt                                                                                     <  2:2895981/93‑1 (MQ=255)
        cagTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGtt                                                                                    <  2:1481638/93‑1 (MQ=255)
         agTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTg                                                                                   <  2:377054/93‑1 (MQ=255)
          gTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgc                                                                                  <  1:3776479/93‑1 (MQ=255)
          gTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgc                                                                                  >  1:4080471/1‑93 (MQ=255)
          gTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgc                                                                                  >  1:6290900/1‑93 (MQ=255)
          gTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgc                                                                                  <  1:1466428/93‑1 (MQ=255)
           tGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgcc                                                                                 >  2:385369/1‑93 (MQ=255)
           tGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTgcc                                                                                 >  1:1084576/1‑93 (MQ=255)
               gCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGccgccg                                                                             <  1:2642599/93‑1 (MQ=255)
                cAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGa                                                                            <  2:3877171/93‑1 (MQ=255)
                  aGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGccc                                                                          <  1:1694145/93‑1 (MQ=255)
                  aGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGAcc                                                                          <  2:4611371/93‑1 (MQ=255)
                     acaAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGa                                                                       <  2:1495470/93‑1 (MQ=255)
                       aaaCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGaaa                                                                     <  2:4719059/93‑1 (MQ=255)
                       aaaCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGACGACCTGaaa                                                                     >  2:901181/1‑93 (MQ=255)
                        aaCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAag                                                                    <  2:491017/93‑1 (MQ=255)
                         aCCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAaga                                                                   <  1:1176677/93‑1 (MQ=255)
                          ccAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagat                                                                  >  2:1657998/1‑93 (MQ=255)
                          ccAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagat                                                                  <  1:5176935/93‑1 (MQ=255)
                           cAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAagata                                                                 <  1:4854503/93‑1 (MQ=255)
                           cAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagata                                                                 <  1:660380/93‑1 (MQ=255)
                            aGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAagatag                                                                <  1:715563/93‑1 (MQ=255)
                             ggCAGCGTTGCCTTTGCCATAAGGCGCGATAATCCCCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagataga                                                               <  2:3301631/93‑1 (MQ=255)
                             ggCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagataga                                                               <  1:3375192/93‑1 (MQ=255)
                             ggCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagataga                                                               >  1:4010877/1‑93 (MQ=255)
                              gCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAagatagat                                                              <  2:3699971/93‑1 (MQ=255)
                               gagcGTTGCCTTTGCCATAAGGCGCGATAAGCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGAtt                                                             <  2:240322/90‑1 (MQ=255)
                               cAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGAtt                                                             >  1:3391811/1‑93 (MQ=255)
                               cAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGAtt                                                             >  1:2165197/1‑93 (MQ=255)
                               cAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGAtt                                                             <  2:5249991/93‑1 (MQ=255)
                                 aCGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTgc                                                           <  2:2215664/93‑1 (MQ=255)
                                 aCGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTgc                                                           >  2:2664209/1‑93 (MQ=255)
                                   gTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTgcgc                                                         <  2:6504713/93‑1 (MQ=255)
                                    ttGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGcc                                                        <  2:4078243/93‑1 (MQ=255)
                                     tGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCa                                                       <  1:6116034/93‑1 (MQ=255)
                                      tccTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAg                                                      <  2:5749649/92‑1 (MQ=255)
                                         tttGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGcc                                                   >  1:1393061/1‑93 (MQ=255)
                                           tGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCAt                                                 <  2:1952739/93‑1 (MQ=255)
                                           tGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCAt                                                 <  1:4452913/93‑1 (MQ=255)
                                           tGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCAt                                                 <  2:624121/93‑1 (MQ=255)
                                             ccATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgc                                               <  1:3551646/93‑1 (MQ=255)
                                             ccATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgc                                               >  2:2876930/1‑93 (MQ=255)
                                              cATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgcc                                              <  2:6112426/93‑1 (MQ=255)
                                              cATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgcc                                              <  1:894581/93‑1 (MQ=255)
                                              cATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgcc                                              >  2:1897482/1‑93 (MQ=255)
                                               aTAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccg                                             >  1:2414851/1‑93 (MQ=255)
                                               aTAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccg                                             >  1:65168/1‑93 (MQ=255)
                                               aTAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccg                                             <  2:313052/93‑1 (MQ=255)
                                               aTAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccg                                             <  1:2435697/93‑1 (MQ=255)
                                                tAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccgc                                            >  2:4444256/1‑93 (MQ=255)
                                                tAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccgc                                            >  1:473913/1‑93 (MQ=255)
                                                tAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATgccgc                                            >  1:1108660/1‑93 (MQ=255)
                                                tAAGGCGCGATAATCACCGCAGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGCAAGATAGATTGCGCCAGGCCATgccgc                                            >  1:5573839/1‑93 (MQ=255)
                                                  aGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGccgccc                                          <  2:2563014/93‑1 (MQ=255)
                                                  aGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGccgccc                                          >  2:5326647/1‑93 (MQ=255)
                                                   ggCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGccgcccg                                         >  1:1996845/1‑93 (MQ=255)
                                                     cgcgATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCcgcccgcc                                       >  2:6318444/1‑93 (MQ=255)
                                                     cgcgATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCcgcccgcc                                       <  1:1868990/93‑1 (MQ=255)
                                                       cgATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCgg                                     >  2:4598438/1‑93 (MQ=255)
                                                        gATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGa                                    <  2:1684496/93‑1 (MQ=255)
                                                          tAATCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAg                                  >  2:5838597/1‑93 (MQ=255)
                                                           aaTCGCCGCCGCCAGCAAGGGTCCAAGGGAGCCTCCAAAGTTGCCGCCGACGTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGc                                 <  1:2744337/93‑1 (MQ=255)
                                                           aaTCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGc                                 >  2:5520492/1‑93 (MQ=255)
                                                           aaTCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGc                                 >  2:6407583/1‑93 (MQ=255)
                                                            aTCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGcc                                <  1:6122843/93‑1 (MQ=255)
                                                            aTCACCGCCGCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGca                                >  2:2089781/1‑92 (MQ=255)
                                                             tCACCGCCCCCACCAAGGGTCCCAGGGACCCTCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCa                               <  1:2587761/93‑1 (MQ=255)
                                                                ccgccgCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATAc                            <  2:3239168/93‑1 (MQ=255)
                                                                ccgccgCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATAc                            <  1:564026/93‑1 (MQ=255)
                                                                ccgccgCCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATAc                            >  2:864093/1‑93 (MQ=255)
                                                                 cgccgcCAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACg                           <  2:584448/93‑1 (MQ=255)
                                                                   ccgccAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACggg                         <  2:1628490/93‑1 (MQ=255)
                                                                    cgccAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGc                        <  1:1560885/93‑1 (MQ=255)
                                                                       cAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAc                     <  2:3386087/93‑1 (MQ=255)
                                                                       cAGCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAc                     >  1:4997373/1‑93 (MQ=255)
                                                                         gCAAGGGTCCCAGGGACCCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAcgc                   <  2:1295912/93‑1 (MQ=255)
                                                                         gCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAcgc                   <  2:5717796/93‑1 (MQ=255)
                                                                         gCAAGGGTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCAcgc                   <  2:474143/93‑1 (MQ=255)
                                                                               gTCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGagcag             >  1:2785964/1‑93 (MQ=255)
                                                                                tCCCAGGGAACCGCCAAAGTTGCTGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGagaaga            >  1:5339433/1‑93 (MQ=255)
                                                                                tCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGagaaga            <  1:4287384/93‑1 (MQ=255)
                                                                                tCCCAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGagaaga            >  2:4255353/1‑93 (MQ=255)
                                                                                 cccAGGGACCTGCCAAAGTTCCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGAt           <  1:3223362/93‑1 (MQ=255)
                                                                                 cccAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGTCCGCCGGAAGCCATACGGGCCACGCGAGAAGAt           <  2:6482017/93‑1 (MQ=255)
                                                                                  ccAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGAtt          >  1:5401012/1‑93 (MQ=255)
                                                                                  ccAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGAtt          >  1:3398972/1‑93 (MQ=255)
                                                                                   cAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTc         <  2:2338035/93‑1 (MQ=255)
                                                                                   cAGGGAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTc         >  2:6587865/1‑93 (MQ=255)
                                                                                     gggAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCg       <  2:4327941/93‑1 (MQ=255)
                                                                                     gggAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCg       <  1:1631511/93‑1 (MQ=255)
                                                                                      ggAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCgg      >  1:37596/1‑93 (MQ=255)
                                                                                       gAACCGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGa     <  1:6353390/93‑1 (MQ=255)
                                                                                          cTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGa  >  1:3103224/1‑93 (MQ=255)
                                                                                          cTGCCAAAGTGGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGa  <  2:808590/93‑1 (MQ=255)
                                                                                           |                                                                                           
ATCGCCAGCAGTGCCGCAAGCACAAACCAGGCAACGTTGCCTTTGCCATAAGGCGCGATAATCACCGCCGCCAGCAAGGGTCCCAGGGAACTGCCAAAGTTGCCGCCGACCTGAAAGATAGATTGCGCCAGGCCATGCCGCCCGCCGGAAGCCATACGGGCCACGCGAGAAGATTCCGGATGA  >  NC_000913/503349‑503531

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